Quantifying Confidence in Horizon-Picking

Risk analysis is a crucial task in making drilling decisions and involves many factors, such as well logs, modeling results, production maps and interpretation quality.

In his book on 3-D seismic interpretation, AAPG award-winning member Alistair Brown presents a workflow for the quantification of interpretation confidence. In this workflow, picks at 0, 1, and 2s indicated low, medium and high reflector quality. The interpreter then generates a confidence map from a coarse grid of picked lines.

In practice, such interpretation confidence maps are commonly excluded from risk analysis, simply because such quantification is not easy.

In this article we demonstrate the quantification of horizon-picking confidence, using two seismic attributes that are sensitive to chaotic features - namely the Sobel-filter and disorder attributes.

Our study area is located within the Halten Terrace, Norwegian North Sea. The area involves rift-related geologic structure, particularly a system of listric faults with a weak, soft layer of salt between basement and the upper sedimentary rocks.

Figure 1a shows the time structure map of an interpreted horizon in the study area.

Figure 2 shows representative vertical slices through the seismic amplitude data.

While the horizon is relatively easy to pick in many areas, there are other areas where it is contaminated by steeply dipping migration alias artifacts. Autopickers work poorly on this horizon.

In order to quantify the confidence of the horizon picking task, we calculate attributes that are sensitive to chaotic features, such as salt, karst and seismic noise. The general idea is that the noisier the data, the less confidence the interpreter will have in picking a horizon.

The Sobel-filter implementation of coherence (the same Sobel filter as in your digital camera software) independently computes first derivatives of the seismic amplitudes between neighboring traces along the X and Y directions and combines them to form a coherence-like image. Disorder, on the other hand, cascades second derivatives in the X, Y and time directions.

Coherence algorithms are designed to emphasize continuous reflectors disrupted by incoherent structural and stratigraphic edges. In contrast, the disorder algorithm is design to emphasize noise and considers edges to be signal.

Both noise estimates are computed along local reflector dip and are normalized by the energy of the data within the analysis window.

Figures 1b and 1c show the results of the Sobel filter and disorder attributes extracted and smoothed along the same horizon in figure 1a. Most of the horizon corresponds to relatively low coherence and high disorder, suggesting that seismic data quality is generally low.

Such data quality impacts the continuity of time-structure maps.

In line AA' shown in figure 2a, the right part of the image corresponds to a smooth time-structure map and high values of coherence and low values of disorder (appearing as green in figures 1b and c) corresponding to a smoother part of the map in figure 1a.

In contrast, line CC' in figure 1c exhibits poor data quality at the target horizon that gives rise to lower coherence and higher disorder displayed as yellow and red in figures 1b and c, and also results in a less smooth time-structure map in figure 1a.

Interestingly, the horizon on the west side of line CC' (figure 1c), shows high coherence (in green) but medium disorder (in yellow). Note that while the horizon is picked as a (white) peak it is overlain by a higher coherence event that appears as a (black) trough. The coherence algorithm appears to measure the continuity of this higher amplitude neighboring reflector.

In this example, the disorder attribute represents data quality more accurately.

In summary, seismic attributes that are sensitive to chaotic features and noisy data, such as coherence and disorder, can be used to quantify horizon-picking confidence. Of the two attributes, disorder is relatively insensitive to faults and provides the more accurate result.

While both attributes are a measure of data quality along a picked reflector, they are not a measure of erroneously picking a more coherent neighboring reflector. Such interpreter error may be the biggest risk of all in the final map.

Authors' note: Thanks to Debapriya Paul for providing geologic information and seismic interpretation data of the study area. AASPI and Petrel were used in this project. Seismic data were provided courtesy of CGG.

(Editor's note: AAPG member Thang Ha is a master's student in geophysics at the University of Oklahoma; AAPG member Kurt Marfurt is his adviser there.)

Figure 2 - Three vertical slices through the seismic amplitude volume showing the yellow picks used to make the map in figure 1a. (a) In line AA' the horizon on the east (right) side is relatively continuous and easy to pick. (b) In line BB' the data quality is poor along the entire picked line. (c) In line CC' the left side of yellow horizon is also noisy, but corresponds to a high coherence (green) area in figure 1b. In this example, the coherence map is sensitive to the overlying, higher amplitude continuous (black) trough.

Figure 1- (a) Time-structure map of yellow horizon shown in figure 2 (below), and corresponding horizon slices through the (b) coherence and (c) Disorder volumes. Coherence is sensitive to structural and stratigraphic edges as well as noise. By design, disorder is insensitive to edges and only sensitive to chaotic noise.

Comments (0)


Geophysical Corner - Thang Ha

AAPG member Thang Ha is a master's student at the University of Oklahoma, Norman, Okla.

Geophysical Corner - Kurt Marfurt
AAPG member Kurt J. Marfurt is with the University of Oklahoma, Norman, Okla.

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 volume provides an overview of current research on fossil and modern dinoflagellates, as well as highlighting research areas for future collaboration, following the DINO9 International Conference in Liverpool. The volume is organized into four themes, with a review paper for each theme written by the key-note speaker. Each theme also includes a future research foci note following discussion during the conference.

Desktop /Portals/0/PackFlashItemImages/WebReady/book-Biological-and-Geological-Perspectives-of-Dinoflagellates-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 8358 Book

See Also: Bulletin Article

Organic-carbon–rich shales of the lower Marcellus Formation were deposited at the toe and basinward of a prograding clinothem associated with a Mahantango Formation delta complex centered near Harrisburg, Pennsylvania. Distribution of these organic-carbon–rich shales was influenced by shifts in the delta complex driven by changes in rates of accommodation creation and by a topographically high carbonate bank that formed along the Findlay-Algonquin arch during deposition of the Onondaga Formation. Specifically, we interpret the Union Springs member (Shamokin Member of the Marcellus Formation) and the Onondaga Formation as comprising a single third-order depositional sequence. The Onondaga Formation was deposited in the lowstand to transgressive systems tract, and the Union Springs member was deposited in the transgressive, highstand, and falling-stage systems tract. The regional extent of parasequences, systems tracts, and the interpreted depositional sequence suggest that base-level fluctuations were primarily caused by allogenic forcing—eustasy, climate, or regional thermal uplift or subsidence—instead of basement fault reactivation as argued by previous workers. Paleowater depths in the region of Marcellus Formation black mudrock accumulation were at least 330 ft (100 m) as estimated by differences in strata thickness between the northwestern carbonate bank and basinal facies to the southeast. Geochemical analysis indicates anoxic to euxinic bottom-water conditions. These conditions were supported by a deep, stratified basin with a lack of circulation.
Desktop /Portals/0/PackFlashItemImages/WebReady/sequence-stratigrapy-and-depositional-environments-of.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 7963 Bulletin Article

A seismic interpretation study of reflection profiles was undertaken to identify sedimentary and structural components of the tectonic development of the West Greenland continental margin. The new structural framework reveals a long and complex tectono-stratigraphy and evolution.

Desktop /portals/0/images/bulletin/2015/09sep/Tectonic-variation-and-structural-evolution-of-.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 22512 Bulletin Article

See Also: DL Abstract

Using examples from shale reservoirs worldwide, I demonstrate the diversity of shale-hosted fracture systems and present evidence for how and why various fractures systems form. Core and outcrop observations, strength tests on shale and on fractures in core, and geomechanical models allow prediction of fracture patterns and attributes that can be taken into account in well placement and hydraulic fracture treatment design. Both open and sealed fractures can interact with and modify hydraulic fracture size and shape. Open fractures can enhance reservoir permeability but may conduct treatment fluids great distances, in some instances possibly aseismically.

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

See Also: Map

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