Finding the Value of 'Impossible'

Two constraints can be imposed on the frequency behavior of a digital seismic wiggle trace:

  • Each frequency component of the trace is a positive quantity.

Frequencies can have a negative algebraic sign in some mathematical operations, such as Fourier transforms, but when the oscillatory behavior of an actual wiggle trace is visually examined to estimate seismic frequency, the frequency estimate is always a positive number.

The situation is much like trying to dial a radio to a negative frequency to listen to music. Negative frequencies are just not a part of the real world.

  • No frequency component can exceed the Nyquist limit fN, which is defined as fN = 1/(2Δt), where Δt is the time sample interval of the digital data in units of seconds.

Because of these constraints, when an instantaneous frequency attribute is calculated along the time extent of a seismic trace, the frequency value should never be negative at any time coordinate, nor should the value ever exceed the Nyquist limit.

Instantaneous frequencies, however, always exhibit negative values and values greater than Nyquist at numerous locations throughout a 3-D seismic data volume.

Please log in to read the full article

Two constraints can be imposed on the frequency behavior of a digital seismic wiggle trace:

  • Each frequency component of the trace is a positive quantity.

Frequencies can have a negative algebraic sign in some mathematical operations, such as Fourier transforms, but when the oscillatory behavior of an actual wiggle trace is visually examined to estimate seismic frequency, the frequency estimate is always a positive number.

The situation is much like trying to dial a radio to a negative frequency to listen to music. Negative frequencies are just not a part of the real world.

  • No frequency component can exceed the Nyquist limit fN, which is defined as fN = 1/(2Δt), where Δt is the time sample interval of the digital data in units of seconds.

Because of these constraints, when an instantaneous frequency attribute is calculated along the time extent of a seismic trace, the frequency value should never be negative at any time coordinate, nor should the value ever exceed the Nyquist limit.

Instantaneous frequencies, however, always exhibit negative values and values greater than Nyquist at numerous locations throughout a 3-D seismic data volume.

The term “impossible frequencies” is used in this article to refer to these anomalous frequency values that violate the constraints just stated.


Figure 1 shows an example of impossible frequencies occurring in a 3-D seismic volume; the surface being analyzed in this instance is a horizon that passes through the Brushy Canyon formation in an area of the Delaware Basin in New Mexico.

At this prospect area, the Brushy Canyon is a sequence of fine-grained siliciclastic units deposited in a deepwater turbidite environment. The color bar used to display the instantaneous frequency behavior was chosen to emphasize impossible frequencies, with all negative frequencies shown as bright red and all frequencies that equal or exceed the Nyquist limit (125 Hz in this case) shown as darkest black.

The arrows positioned along arbitrary lines A, XL91 and XL57 mark where each profile crosses a reasonably continuous trend of these impossible frequencies.

Vertical profiles along lines A and XL57 are displayed as figure 2 to illustrate the wavelet character that is associated with these impossible frequencies. The position of the attribute surface in figure 1 is 10 ms below the K horizon along these profiles.

The location of each anomalous-frequency arrow shown along the map horizon is indicated by a circle.

At each of these circled coordinates there is a reflection termination or some type of abrupt change in the waveform of the reflection wavelet. Interpreters working with these particular attribute displays concluded that most of these anomalous frequencies were marking stratigraphic pinch-outs or reservoir compartment boundaries.

If true, then “impossible frequencies” have great interpretational value.


Applying the assumption that, in this instance, impossible frequencies define boundaries of reservoir compartments, a tentative reservoir-compartment map such as that shown in figure 3 can be constructed.

The compartments on this map were created by drawing boundaries along trends suggested by the alignments of the impossible frequencies. Each boundary defined by this “connect-the-dots” method is arbitrary; other interpreters would connect the dots in different ways and cause some compartments to have sizes and shapes that differ from those shown on this map.

Is the assumption that impossible frequencies define compartment boundaries valid in this case?

At this prospect, there is circumstantial evidence that the answer is “yes.”

  • Drilling has confirmed that wells in area A penetrate small reservoir compartments at this stratigraphic level – the same principle implied by the frequency-based map.
  • A horizontal well drilled across area B has been the most prolific Brushy Canyon producer across a broad area of the basin, confirming that the well penetrated a large compartment, as indicated by the compartment map.

Will impossible frequencies be valuable indicators of stratigraphic terminations and/or compartment boundaries in other areas and in other depositional environments?

Always consider the possibility that they will.

Incidentally, some commercial interpretation software arbitrarily deletes the negative algebraic sign of any negative frequencies created by their algorithms, which may not be a desirable thing to do. Those negative frequencies can be important information.

You may also be interested in ...