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

Alternative Resources, Structure, Geochemistry and Basin Modeling, Sedimentology and Stratigraphy, Geophysics, Business and Economics, Engineering, Petrophysics and Well Logs, Environmental, Geomechanics and Fracture Analysis, Compressional Systems, Salt Tectonics, Tectonics (General), Extensional Systems, Fold and Thrust Belts, Structural Analysis (Other), Basin Modeling, Source Rock, Migration, Petroleum Systems, Thermal History, Oil Seeps, Oil and Gas Analysis, Maturation, Sequence Stratigraphy, Clastics, Carbonates, Evaporites, Seismic, Gravity, Magnetic, Direct Hydrocarbon Indicators, Resource Estimates, Reserve Estimation, Risk Analysis, Economics, Reservoir Characterization, Development and Operations, Production, Structural Traps, Oil Sands, Oil Shale, Shale Gas, Coalbed Methane, Deep Basin Gas, Diagenetic Traps, Fractured Carbonate Reservoirs, Stratigraphic Traps, Subsalt Traps, Tight Gas Sands, Gas Hydrates, Coal, Uranium (Nuclear), Geothermal, Renewable Energy, Eolian Sandstones, Sheet Sand Deposits, Estuarine Deposits, Fluvial Deltaic Systems, Deep Sea / Deepwater, Lacustrine Deposits, Marine, Regressive Deposits, Transgressive Deposits, Shelf Sand Deposits, Slope, High Stand Deposits, Incised Valley Deposits, Low Stand Deposits, Conventional Sandstones, Deepwater Turbidites, Dolostones, Carbonate Reefs, (Carbonate) Shelf Sand Deposits, Carbonate Platforms, Sebkha, Lacustrine Deposits, Salt, Conventional Drilling, Directional Drilling, Infill Drilling, Coring, Hydraulic Fracturing, Primary Recovery, Secondary Recovery, Water Flooding, Gas Injection, Tertiary Recovery, Chemical Flooding Processes, Thermal Recovery Processes, Miscible Recovery, Microbial Recovery, Drive Mechanisms, Depletion Drive, Water Drive, Ground Water, Hydrology, Reclamation, Remediation, Remote Sensing, Water Resources, Monitoring, Pollution, Natural Resources, Wind Energy, Solar Energy, Hydroelectric Energy, Bioenergy, Hydrogen Energy
Desktop /Portals/0/PackFlashItemImages/WebReady/book-s59-Carbon-Dioxide-Sequestration-in-Geological-Media.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 5827 Book

See Also: Bulletin Article

A new hierarchical architectural classification for clastic marginal-marine depositional systems is presented and illustrated with examples. In ancient rocks, the architectural scheme effectively integrates the scales of sedimentology (core, outcrop) and sequence stratigraphy (wireline-log correlation, reflection seismic). The classification also applies to modern sediments, which allows for direct comparison of architectural units between modern and ancient settings. In marginal-marine systems, the parasequence typically defines reservoir flow units. This classification addresses subparasequence scales of stratigraphy that commonly control fluid flow in these reservoirs. The scheme consists of seven types of architectural units that are placed on five architectural hierarchy levels: hierarchy level I: element (E) and element set (ES); hierarchy level II: element complex (EC) and element complex set (ECS); hierarchy level III: element complex assemblage (ECA); hierarchy level IV: element complex assemblage set (ECAS); and hierarchy level V: transgressive-regressive sequence (T-R sequence). Architectural units in levels I to III are further classified relative to dominant depositional processes (wave, tide, and fluvial) acting at the time of deposition. All architectural units are three-dimensional and can also be expressed in terms of plan-view and cross-sectional geometries. Architectural units can be linked using tree data structures by a set of familial relationships (parent-child, siblings, and cousins), which provides a novel mechanism for managing uncertainty in marginal-marine systems. Using a hierarchical scheme permits classification of different data types at the most appropriate architectural scale. The use of the classification is illustrated in ancient settings by an outcrop and subsurface example from the Campanian Bearpaw–Horseshoe Canyon Formations transition, Alberta, Canada, and in modern settings, by the Mitchell River Delta, northern Australia. The case studies illustrate how the new classification can be used across both modern and ancient systems, in complicated, mixed-process depositional environments.
Desktop /Portals/0/PackFlashItemImages/WebReady/a-hierarchical-approach-to-architectural-classification.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3769 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 4315 CD-DVD

See Also: DL Abstract

When evaluating paleosystems, there will always be a shortage of data constraints and a surplus of plausible geological scenarios for a basin evaluation. Modelling paleosystems with constraints from the modern has been used as a successful approach to better understand petroleum systems.

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

See Also: Short Course

Here is an introduction to the tools and techniques that geologists and geophysicists use to locate gas and oil, that drillers use to drill the wells and that petroleum engineers use to test and complete the wells and produce the gas and oil. Exercises throughout the course provide practical experience in well log correlation, contouring, interpretation of surface and subsurface, contoured maps, seismic interpretation, well log interpretation, and decline curve analysis.

Desktop /Portals/0/PackFlashItemImages/WebReady/sc-basic-petroleum-geology-for-the-non-geologist.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 13584 Short Course