Complex Traces: There’s an ‘App’ for That

Last month we introduced the concept of a complex seismic trace; this month we’ll show how a complex trace provides a rigorous way to set the boundaries of data windows associated with distinct seismic reflections – and we’ll define the polarities of each of those reflection events.

This complex trace application is important because it is necessary to determine the polarity of every reflection event that spans a layered system in order to determine whether impedance increases or decreases from layer to layer – which in turn provides insight into the lithology, porosity and type of pore fluid in each rock layer.

The principal problem involved in determining the polarity of a seismic reflection event is the challenge of deciding what part of the seismic response represents the reflection event.

Questions that have to be answered include:

  • Where does the reflection event start and stop?
  • How many peaks and troughs are embedded in the reflection event?
  • Which peak or trough of a reflection event should be used to define reflection polarity?

The amplitude-envelope function determined from a complex seismic trace provides a way to define the start time, stop time, wavelet character and polarity of overlapping – but distinct – reflection events.

An example seismic trace, its complex-trace equivalent and the associated amplitude envelope are shown as figure 1.

Defining Reflection Events

As shown on figure 1, the amplitude envelope of a complex seismic trace is an oscillating function that has alternating maxima and minima. The data window between two successive minima of an amplitude-envelope function defines a distinct packet of seismic energy.

Terms that have been used to describe this interval between successive amplitude-envelope minima are energy packet, wavelet packet and reflection event. Once you equate the term “reflection event” with energy packet (or with wavelet packet), you can then ask the question:

“How many reflection events occur between time coordinates A and B on figure 1?”

You will get the definitive answer “13.”

A wavelet packet such as any of those defined on figure 1 may be a reflection from a single interface, or it may be a composite of several reflections from closely spaced interfaces. In either case, a wavelet packet represents the shortest-time concentration of reflection of energy that can be recognized in a seismic response.

Because amplitude-envelope minima can be determined numerically after an amplitude envelope is calculated, the start time, stop time and time extent of a reflection event can be defined with mathematical rigor, as shown by each of the labeled “events” on figure 1, and do not have to be left to interpreter judgment.

The basic seismic wavelet that is embedded in the seismic trace on the left of figure 1 is shown in the center part of the figure. A reader can compare this wavelet with its associated reflection trace on the left of the display to attempt to decide how many reflection events exist across the time interval A to B.

In classroom and workshop exercises, people have tended to conclude that the number of reflection events ranges from a low of five or six to a high of 17 or 18. Using the mathematical concept of amplitude-envelope minima to define the boundaries of a reflection event, the correct answer is 13 reflection events (right-hand panel) as already stated.

Defining Reflection Polarity

When a reflection event is defined by this energy packet concept, the polarity of the reflection event can be defined as the algebraic sign of the real-trace extremum (either peak or trough) that is closest to the maximum of the amplitude-envelope that encompasses the energy packet.

Using this concept, the polarity of reflection events 5 and 10 on figure 1 are positive, and the polarities of reflection events 7 and 12 are negative. Thus a complex-trace allows seismic reflection polarity to be defined with the same mathematical rigor that defines the time extent of each reflection event.

A second illustration of energy packets being used to define distinct reflection events and their polarities is provided as figure 2.

In this case, there are excellent examples of energy packets distinguishing overlapping reflection events (events 7 and 8 and events 11 and 12) and defining the data windows spanned by faint, low-amplitude reflections (events 2 and 3).

Comments (0)


Geophysical Corner

The Geophysical Corner is a regular column in the EXPLORER that features geophysical case studies, techniques and application to the petroleum industry.


Image Gallery

See Also: Book

This special issue honors the legacy of J. Fred Read, a pioneer in carbonate sedimentology and stratigraphy. He taught at Virginia Tech for 38 years and, along with his students, published more than 120 papers. Many of these students have become leaders in carbonate research.

Desktop /Portals/0/images/_site/AAPG-newlogo-vertical-morepadding.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12465 Book

See Also: Bulletin Article

We describe the structure, microstructure, and petrophysical properties of fault rocks from two normal fault zones formed in low-porosity turbiditic arkosic sandstones, in deep diagenesis conditions similar to those of deeply buried reservoirs. These fault rocks are characterized by a foliated fabric and quartz-calcite sealed veins, which formation resulted from the combination of the (1) pressure solution of quartz, (2) intense fracturing sealed by quartz and calcite cements, and (3) neoformation of synkinematic white micas derived from the alteration of feldspars and chlorite. Fluid inclusion microthermometry in quartz and calcite cements demonstrates fault activity at temperatures of 195degC to 268degC. Permeability measurements on plugs oriented parallel with the principal axes of the finite strain ellipsoid show that the Y axis (parallel with the foliation and veins) is the direction of highest permeability in the foliated sandstone (10–2 md for Y against 10–3 md for X, Z, and the protolith, measured at a confining pressure of 20 bars). Microstructural observations document the localization of the preferential fluid path between the phyllosilicate particles forming the foliation. Hence, the direction of highest permeability in these fault rocks would be parallel with the fault and subhorizontal, that is, perpendicular to the slickenlines representing the local slip direction on the fault surface. We suggest that a similar relationship between kinematic markers and fault rock permeability anisotropy may be found in other fault zone types (reverse or strike-slip) affecting feldspar-rich lithologies in deep diagenesis conditions.
Desktop /Portals/0/PackFlashItemImages/WebReady/Bulletin-hero-2013-07jul.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5694 Bulletin Article

See Also: Energy Policy Blog

18 percent of AAPG members are women, up from 10 percent in 2006.

Desktop /Portals/0/PackFlashItemImages/WebReady/future-oil-and-gas-opportunities-for-women-and-minorities-2014-04apr-16-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12894 Energy Policy Blog

See Also: Online e Symposium

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-eagle-ford-shale-prospecting-with-3d-seismic-data.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1482 Online e-Symposium

This presentation discusses one operator’s approach to fully integrate data captured in the Marcellus Shale in order to optimize horizontal well performance.

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-marcellus-shale-integrating-information-to-optimize.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1476 Online e-Symposium