Some Seismic Instantaneous Bandwidth Attribute Applications

The bandwidth of a signal is a measure of the breadth of its spectrum, often expressed as the distance measured in Hertz (or better still, in octaves) between the half-power points of a smooth spectrum.

Bandwidth is the best measure of resolution. Narrow-band, low-frequency spectra cannot represent thin beds, while narrow-band, high-frequency spectra provide ringy, ambiguous images of thin beds. Ideally, the bandwidth is estimated by directly measuring the data spectrum using either Fourier transforms or spectral decomposition. In contrast, most interpretation packages include an instantaneous bandwidth attribute defined as the ratio of the derivative of the instantaneous envelope to the instantaneous envelope itself. An alternative measure is to estimate two times the standard deviation of the spectrum about its mean. The mathematics of both approximations assume the seismic data exhibit an approximately Gaussian spectrum – a poor assumption if the processing shop has worked hard to make the spectrum as flat as possible.

Given this limitation, instantaneous bandwidth still can be quite useful.

Please log in to read the full article

The bandwidth of a signal is a measure of the breadth of its spectrum, often expressed as the distance measured in Hertz (or better still, in octaves) between the half-power points of a smooth spectrum.

Bandwidth is the best measure of resolution. Narrow-band, low-frequency spectra cannot represent thin beds, while narrow-band, high-frequency spectra provide ringy, ambiguous images of thin beds. Ideally, the bandwidth is estimated by directly measuring the data spectrum using either Fourier transforms or spectral decomposition. In contrast, most interpretation packages include an instantaneous bandwidth attribute defined as the ratio of the derivative of the instantaneous envelope to the instantaneous envelope itself. An alternative measure is to estimate two times the standard deviation of the spectrum about its mean. The mathematics of both approximations assume the seismic data exhibit an approximately Gaussian spectrum – a poor assumption if the processing shop has worked hard to make the spectrum as flat as possible.

Given this limitation, instantaneous bandwidth still can be quite useful.

In general, temporally short seismic events exhibit high bandwidth, while temporally broad and/or ringing events exhibit narrow bandwidth. Instantaneous bandwidth is a sample-by-sample computation that is commonly used to measure vertical and lateral changes in attenuation. Seismic data undergo two types of attenuation: intrinsic attenuation, in which mechanical energy is converted to heat energy and is only a function the rate of compression/dilation (the frequency); and geometric attenuation, in which interference phenomena give rise to constructive and destructive interference. For random heterogeneities, geometric attenuation decreases the higher frequency components of the spectrum. For more ordered and less commonly encountered heterogeneities such as sabkhas, coal cyclothems, and laterally periodic fractures, some frequencies may be attenuated, and others enhanced. If the lower components of the spectrum are unchanged and the upper part attenuated, the bandwidth decreases.

Although first introduced to the seismic community in the 1970s, only a limited number of published papers show the value of this attribute. In this column, we show the value of bandwidth when applied to data from the Utica area of Eastern Ohio, United States and from the Gulf of Mexico.

Multiple versus Merged Surveys

Many mature areas for hydrocarbon exploration are covered by adjacent or overlapping 3-D seismic surveys acquired by different operators at different times using different technology. These surveys very commonly exhibit different amplitudes frequencies and phase responses in which the differences in frequency and phase might be simple – such as the use of different geophones with different resonant frequencies, or more complex, depending on the details of statics, filtering, velocity analysis and imaging. Ideally, these surveys are reprocessed from the original common shot gathers to provide a merged survey with a consistent datum, amplitude and phase. Although the amplitudes of the measured spectral components are consistent, because of differences in vibrator sweeps and geophone responses, some surveys will have only recorded a narrower spectrum, giving rise to changes in bandwidth across the survey.

There are several methods used during processing to enhance the frequency content of the seismic data. Starting with the application of the traditional time-variant spectral whitening, which might distort the seismic amplitudes, one might use inverse Q-filtering, spectral decomposition-based methods and thin-bed reflectivity inversion.

We have demonstrated the application of some of these methods in previous installments of Geophysical Corner (for example, check out the July 2009 and May 2014 columns). In this article, we demonstrate the application of two such methods of frequency enhancement applied to a merge of two-component seismic surveys. We refer to the two methods as 1 and 2; the performance comparison of these methods is not the focus of this analysis. In figure 1a, we show a segment of a seismic section from the merged surveys that exhibits more high frequency content to the left and less high frequency content on the right, resulting in lower vertical resolution. The overall frequency content of the displayed time window is shown in the inset as well as the average half-power bandwidth for the entire survey. Equivalent sections after the application the methods 1 and 2 are shown in figure 1b and c. Note that the application of method 2 provides a superior image to that of method 1.

Figure 2 shows the instantaneous bandwidth attribute computed for the input and two frequency-enhanced volumes shown in figure 1. Notice the higher instantaneous bandwidth on the left and lower instantaneous bandwidth on the right halves of figures 2a and 2b. In contrast, the instantaneous bandwidth in figure 2c is much better balanced. This application demonstrates that the instantaneous bandwidth attribute can be used as an aid to assess the applications of different frequency-enhancement procedures.

Imaging Salt Bodies

In many regions of the world, oil and gas accumulations are found trapped against salt domes. Consequently, the flanks of salt domes are considered potential exploration targets. For this reason, the boundaries of salt bodies need to be imaged and interpreted well. More recently, object extraction machine learning applications have aimed to pick or interpret salt bodies in an automated fashion, employing neural network applications that rely on differences between the signature of the seismic signal outside and inside the salt domes. In figure 3, the seismic signature within the salt body is lower amplitude and lower frequency content than the sedimentary layers around it. The time slices in figure 3a and 3b show low values of instantaneous envelope and bandwidth within the salt dome and higher values in the sedimentary layers. These observations indicated that both attributes are good candidates to use in training a neural network for salt body extraction.

Quantitative estimates of seismic data quality include measures of continuity, signal-to-noise ratio, temporal resolution. The instantaneous bandwidth attribute is a good first-order measure of temporal resolution and can be used to map lateral and vertical changes in seismic data quality and to discriminate between different seismic facies. The instantaneous bandwidth can also be useful in ascertaining not only the frequency content of the signal, but also to recognize noise that may contaminate the data.

Acknowledgement: The first author would like to thank Gary Margrave for drawing his attention to the instantaneous bandwidth attribute.

Comments (1)

Instantaneous Bandwidth
Strictly speaking, the definition given for instantaneous bandwidth might better be termed "relative amplitude change"; instantaneous bandwidth is defined as the absolute value of the derivative of the trace envelope divided by the trace envelope, and further divided by 2 PI so as to have units of Hz. This definition is independent of any assumption of a Gaussian distributed spectrum. The earliest appearance of this measure is in a 1988 SPIE paper by Cohen and Lee; it does not appear to have become a seismic attribute until 1993. A measure of local bandwidth derived from short-time window spectra as described here may have appeared earlier, but I do not know when it was first published. These two seemingly disparate measures are reconciled in an average sense through the so-called bandwidth formula (described in Cohen's 1995 book on time-frequency analysis). It is this formula that justifies treating relative amplitude change as a measure of bandwidth.
3/26/2020 3:16:38 PM

You may also be interested in ...