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: Bulletin Article

The origin of thermogenic natural gas in the shallow stratigraphy of northeastern Pennsylvania is associated, in part, with interbedded coal identified in numerous outcrops of the Upper Devonian Catskill and Lock Haven Formations. Historically documented and newly identified locations of Upper Devonian coal stringers are shown to be widespread, both laterally across the region and vertically throughout the stratigraphic section of the Catskill and Lock Haven Formations. Coal samples exhibited considerable gas source potential with total organic carbon as high as 44.40% by weight, with a mean of 13.66% for 23 sample locations analyzed. Upper Devonian coal is thermogenically mature; calculated vitrinite reflectances range from 1.25% to 2.89%, with most samples falling within the dry-gas window. Source potential is further supported by gas shows observed while drilling through shallow, identifiable coal horizons, which are at times located within fresh groundwater aquifers. Thermogenic gas detected in area water wells during predrill baseline sampling is determined not only to be naturally occurring, but also common in the region.

Desktop /Portals/0/PackFlashItemImages/WebReady/geologic-and-baseline-groundwater-evidence-for.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5776 Bulletin Article

See Also: CD DVD

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

See Also: DL Abstract

The results of regional deep seismic acquisition in the South Atlantic continental margins have shed new lights on the birth and development of sedimentary basins formed during the Gondwana breakup. Recent models of mantle exhumation as observed in the deep water Iberian margin have been applied extensively to the interpretation of several basins in the Eastern Brazilian and West African conjugate margins. However, the tectonic development of these basins is markedly different from the magma-poor margins, and in this lecture we emphasize the contrasts from the tectono-sedimentary features imaged in deep-penetrating seismic profiles that extend from the platform towards the oceanic crust, which indicate that the Red Sea constitutes a better analogue for the birth of divergent continental margins.

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

See Also: Online e Symposium

Upon successful completion of this course, you will be able to describe geomechanics in shale reservoirs and discuss differences between plays.

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-all-shale-gas-reservoirs-are-not-the-same.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1461 Online e-Symposium