Machine Learning Applied to 3D Seismic Data from the Denver-Julesburg Basin Improves Stratigraphic Resolution in the Niobrara

Machine Learning Applied to 3D Seismic Data from the Denver-Julesburg Basin Improves Stratigraphic Resolution in the Niobrara

By Carolan Laudon, Sarah Stanley, Patricia Santogrossi 
Published with permission: Unconventional Resources Technology Conference (URTeC 2019)
July 2019

Abstract

Seismic attributes can be both powerful and challenging to incorporate into interpretation and analysis. Recent developments with machine learning have added new capabilities to multi-attribute seismic analysis. In 2018, Geophysical Insights conducted a proof of concept on 100 square miles of multi-client 3D data jointly owned by Geophysical Pursuit, Inc. (GPI) and Fairfield Geotechnologies (FFG) in the Denver-Julesburg Basin (DJ). The purpose of the study was to evaluate the effectiveness of a machine learning workflow to improve resolution within the reservoir intervals of the Niobrara and Codell formations, the primary targets for development in this portion of the basin.

The seismic data are from Phase 5 of the GPI/Fairfield Niobrara program in northern Colorado. A preliminary workflow which included synthetics, horizon picking and correlation of 28 wells was completed. The seismic volume was re-sampled from 2 ms to 1 ms. Detailed well time-depth charts were created for the Top Niobrara, Niobrara A, B and C benches, Fort Hays and Codell intervals. The interpretations, along with the seismic volume, were loaded into the Paradise® machine learning application, and two suites of attributes were generated, instantaneous and geometric. The first step in the machine learning workflow is Principal Component Analysis (PCA). PCA is a method of identifying attributes that have the greatest contribution to the data and that quantifies the relative contribution of each. PCA aids in the selection of which attributes are appropriate to use in a Self-Organizing Map (SOM). In this case, 15 instantaneous attribute volumes, plus the parent amplitude volume, were used in the PCA and eight were selected to use in SOMs. The SOM is a neural network-based machine learning process that is applied to multiple attribute volumes simultaneously. The SOM produces a non-linear classification of the data in a designated time or depth window.

For this study, a 60-ms interval that encompasses the Niobrara and Codell formations was evaluated using several SOM topologies. One of the main drilling targets, the B chalk, is approximately 30 feet thick; making horizontal well planning and execution a challenge for operators. An 8 X 8 SOM applied to 1 ms seismic data improves the stratigraphic resolution of the B bench. The neuron classification also images small but significant structural variations within the chalk bench. These variations correlate visually with the geometric curvature attributes. This improved resolution allows for precise well planning for horizontals within the bench. The 25 foot thick C bench and the 17 to 25 foot thick Codell are also seismically resolved via SOM analysis. Petrophysical analyses from wireline logs run in seven wells within the survey by Digital Formation; together with additional results from SOMs show the capability to differentiate a high TOC upper unit within the A marl which presents an additional exploration target. Utilizing 2D color maps and geobodies extracted from the SOMs combined with petrophysical results allows calculation of reserves for the individual reservoir units as well as the recently identified high TOC target within the A marl.

The results show that a multi-attribute machine learning workflow improves the seismic resolution within the Niobrara reservoirs of the DJ Basin and results can be utilized in both exploration and development.

Introduction and preliminary work

The Denver-Julesburg Basin is an asymmetrical foreland basin that covers approximately 70,000 square miles over parts of Colorado, Wyoming, Kansas and Nebraska. The basin has over 47,000 oil and gas wells with a production history that dates back to 1881 (Higley, 2015). In 2009, operators in the Wattenberg field began to drill and complete horizontal wells in the chalk benches of the Niobrara formation and within the Codell sandstone. As of October 2018, approximately 9500 horizontal wells have been drilled and completed within Colorado and Wyoming in the Niobrara and Codell formations (shaleprofile.com/2019/01/29/niobrara-co-wy-update-through-october-2018).

The transition to horizontal drilling necessitated the acquisition of modern, 3D seismic data (long offset, wide azimuth) to properly image the complex faulting and fracturing within the basin. In 2011, Geophysical Pursuit, Inc., in partnership with the former Geokinetics Inc., embarked on a multi-year, multi-client seismic program that ultimately resulted in the acquisition of 1580 square miles of contiguous 3D seismic data. In 2018, Geophysical Pursuit, Inc. (GPI) and joint-venture partner Fairfield Geotechnologies (FFG) provided Geophysical Insights with seismic data in the Denver-Julesburg Basin to conduct a proof of concept evaluation of the effectiveness of a machine learning workflow to improve resolution within the reservoir intervals of the Niobrara and Codell formations, currently the primary targets for development in this portion of the basin. The GPI/FFG seismic data analyzed are 100 square miles from the Niobrara Phase 5 multi-client 3D program in northern Colorado (Figure 1). Prior to the machine learning workflow, a preliminary interpretation workflow was carried out, that included synthetics, horizon picking and well correlation on 28 public wells with digital data. The seismic volume was resampled from 2 ms to 1 ms. Time depth charts were made with detailed well ties for the Top Niobrara, Niobrara A, B, and C benches, Fort Hays and Codell. The interpretations, along with the re-sampled seismic amplitude volume, were loaded into the Paradise® machine learning application. The machine learning software has several options for computing seismic attributes, and two suites were selected for the study: standard instantaneous attributes and geometric attributes from the AASPI (Attribute Assisted Seismic Processing and Interpretation) consortium (http://mcee.ou.edu/aaspi/).

Figure 1: Map of GPI FFG multi-client program and study area outline

Geologic Setting of the Niobrara and Surrounding Formations

The Niobrara formation is late Cretaceous in age and was deposited in the Western Interior Seaway (Kaufmann, 1977). The Niobrara is subdivided into the basal Fort Hays limestone and the Smoky Hill member. The Smoky Hill member is further subdivided into three subunits informally termed Niobrara A, B, and C. These units consist of fractured chalk benches which are primary reservoirs with marls and shales between the benches which comprise source rocks and secondary reservoir targets. (Figure 2). The Niobrara unconformably overlies the Codell sandstone and is overlain by the Sharon Springs member of the Pierre shale.

The Codell is also late Cretaceous in age, and unconformably underlies the Fort Hays member of the Niobrara formation. In general, the Codell thins from north to south due to erosional truncation (Sterling, Bottjer and Smith, 2016). In the study area, the thickness of the Codell ranges from 18 to 25 feet. Lewis (2013) inferred an eastern provenance for the Codell with a limited area of deposition or subsequent erosion through much of the DJ Basin. Based upon geochemical analyses, Sterling and others (2016) state that hydrocarbons produced from the Codell are sourced from the Niobrara, primarily the C marl, and the thermal maturity provides evidence of migration into the Codell. The same study found that oil produced from the Niobrara C chalk was generated in-situ.

Figure 2 (Sonnenberg, 2015) shows a generalized stratigraphic column and a structure map for the Niobrara in the DJ Basin along with an outline of the DJ basin and the location of the Wattenberg Field within which the study area is contained.

Figure 2: Outline of the DJ Basin with Niobrara structure contours and generalized stratigraphic column that shows the source rock and reservoir intervals for late Cretaceous units in the basin (from Sonnenberg, 2015).

Figure 3 shows the structural setting of the Niobrara in the study area, as well as types of fractures which can be expected to provide storage capacity and permeability for reservoirs within the chalk benches (Friedman and others, 1992). The study area covers approximately 100 square miles and shows large antiforms on the western edge. The area is normally faulted with most faults trending northeast to southwest. The Top Niobrara time structure also shows extensive small-scale structural relief which is visualized in a curvature attribute volume as shown in Figure 4. This implies that a significant amount of fracturing is present within the Niobrara.

Figure 3: Gross structure of the Niobrara in the study area in seismic two-way travel time. Insets from Friedman and others, 1992, showing predicted fracture types from structural elements. Area shown is approximately 100 square miles.

Figure 4: Most positive curvature, K1 on top Niobrara. The faulting and fractures are complex with both NE-SW and NW-SE trends apparent. Area shown is approximately 100 square miles. Seismic data provided courtesy of GPI and FFG.

Meissner and others (1984) and Landon and others (2001) have stated that the Niobrara formation kerogen is Type-II and oil-prone. Landon and others, and Finn and Johnson (2005) have also stated that the DJ basin contains the richest Niobrara source rocks with TOC contents reaching eight weight percent. Niobrara petroleum production is dependent on fractures in the hard, brittle, carbonate-rich zones. These zones are overlain and/or interbedded with soft, ductile marine shales that inhibit migration and seal the hydrocarbons in the fractured zones.

Why Utilize Machine Learning?

In the study area, the Niobrara to Greenhorn section is represented in approximately 60 milliseconds of two-way travel time in the seismic data. Figure 5 shows an amplitude section through a well within the study area. Figure 6 is an index map of wells used in the study with the Anderson 11-2 well highlighted in red. It is apparent that the top Niobrara is a well resolved positive amplitude or peak which can be picked on either a normal amplitude section or an instantaneous phase display. The individual units within the Niobrara A bench, A marl, B bench, B marl, C bench, C marl, Fort Hays and Codell present a significant challenge for an interpreter to resolve using only one or two attributes. The use of simultaneous multiple seismic attributes holds promise to resolve thin beds and a machine learning approach is one methodology which has been documented to successfully resolve stratigraphy below tuning (Roden and others, 2015, Santogrossi, 2017).

Figure 5: Amplitude section shows the approximately 60 milliseconds between marked horizons which contain the Niobrara and Codell reservoirs. Trace spacing is 110 feet, vertical scale is two-way time in seconds. Seismic data are shown courtesy of GPI and FFG.

Figure 6: Index map of vertical wells used in study. The dashed lines connect well names to well locations. Wells were obtained from the Colorado Oil and Gas Conservation Commission public database.

Machine Learning Data Preparation

The Niobrara Phase 5 3D data used for this study consisted of a 32-bit seismic amplitude volume that covers approximately 100 square miles. The survey contained 5.118 seconds of data with a bin spacing of 110 feet. Machine learning classifications benefit from sharper natural clusters of information through one level of finer trace sampling. Machine learned seismic resolution also benefits from sample-by-sample classification when compared to conventional wavelet analysis. Therefore, the data were upsampled to 1 ms from its original 2 ms interval by Geophysical Insights. The 1 ms amplitude data were used for seismic attribute generation.

Focus should be placed on the time interval that encompasses the geologic units of interest. The time interval selected for this study was 0.5 seconds to 2.2 seconds.

A total of 44 digital wells were obtained, 40 of which were within the seismic survey.

Classification by Principal Component Analysis (PCA)

Multi-dimensional analysis and multi-attribute analysis go hand in hand. Because individuals are grounded in three-dimensional space, it is difficult to visualize what data in a higher number dimensional space looks like. Fortunately, mathematics doesn’t have this limitation and the results can be easily understood with conventional 2D and 3D viewers.

Working with multiple instantaneous or geometric seismic attributes generates tremendous volumes of data. These volumes contain huge numbers of data points which may be highly continuous, greatly redundant, and/or noisy. (Coleou et al., 2003). Principal Component Analysis (PCA) is a linear technique for data reduction which maintains the variation associated with the larger data sets (Guo and others, 2009; Haykin, 2009; Roden and others, 2015). PCA has the ability to separate attribute types by frequency, distribution, and even character. PCA technology is used to determine which attributes to use and which may be ignored due to their very low impact on neural network solutions.

Figure 7 illustrates the analysis of a data cluster in two directions offset by 90 degrees. The first principal component (eigenvector 1) analyses the data cluster along the longest axis. The second principal component (eigenvector 2) analyses the data cluster variations perpendicular to the first principal component. As stated in the diagram, each eigenvector is associated with an eigenvalue which shows how much variance is in the data.

Figure 7: 2 attribute data set demonstrating the concept of PCA

Eigenvectors and eigenvalues from inline 1683 were consistently used for Principal Component Analysis because line 1683 bisected the deepest well in the study area. The entire pre-Niobrara, Niobrara, Codell, and post-Niobrara depositional events were present in the borehole.

PCA results for the first two eigenvectors for the interval Top Niobrara to Top Greenhorn are shown in Figure 8. Results show the most significant attributes in the first eigenvector are Sweetness, Envelope, and Relative Acoustic Impedance; each contributes approximately 60% of the maximum value for the eigenvector. PCA results for the second eigenvector show Thin Bed and Instantaneous Frequency are the most significant attributes. Figure 9 shows instantaneous attributes from the first eigenvector (sweetness) and second eigenvector (thin bed indicator) extracted near the B chalk of the Niobrara. The table shown in Figure 9 lists the instantaneous attributes that PCA indicated contain the most significance in the survey and the eigenvector associated with the attribute. This selection of attributes comprises a ‘recipe’ for input to the Self-Organizing Maps for the interval Niobrara to Greenhorn.

Figure 8: Eigenvalue charts for Eigenvectors 1 and 2 from PCA for Top Niobrara to Top Greenhorn. Attributes that contribute more than 50% of the maximum were selected for input to SOM

Figure 9: Instantaneous attributes near the Niobrara B chalk. These are prominent attributes in Eigenvectors 1 and 2. On the right of the figure is a list of eight selected attributes for SOM analysis. Seismic data is shown courtesy of GPI and FFG.

Self-Organzing Maps

Teuvo Kohonen, a Finnish mathematician, invented the concepts of Self-organizing Maps (SOM) in 1982 (Kohonen, T., 2001). Self-Organizing Maps employ the use of unsupervised neural networks to reduce very high dimensions of data to a scale that can be easily visualized (Roden and others, 2015). Another important aspect of SOMs is that every seismic sample is used as input to classification as opposed to wavelet-based classification.

Figures 10 and 11 illustrate classification by SOM. Within the 3D seismic survey, samples are first organized into attribute points with similar properties called natural clusters in attribute space. Within each cluster new, empty, multi-attribute samples, named neurons, are introduced. The SOM neurons will seek out natural clusters of like characteristics in the seismic data and produce a 2D mesh that can be illustrated with a two- dimensional color map.

Figure 10: Example SOM classification of two attributes into 4 clusters (neurons)

In other words, the neurons “learn” the characteristics of a data cluster through an iterative process (epochs) of cooperative then competitive training. When the learning is completed each unique cluster is assigned to a neuron number and each seismic sample is now classified (Smith, 2016).

Figure 11: Illustration of how SOM works with 3D seismic volumes

Note that the two-dimensional color map in Figure 11 shows an 8X8 topology. Topology is important. The finer the topology of the two-dimensional color map the finer the data clusters associated with each neuron become. For example: an 8X8 topology distributes 64 neurons throughout an attribute set, while a 12X12 topology distributes 144 neurons. Finer topologies help to refine variations in lithologies, porosity, and other reservoir characteristics. Although there is no theoretical limit to a two-dimensional map topology, experience has shown that there is a practical limit to the number of neuron topologies for geological resolution. Conversely, a coarser neuron topology is associated with much larger data clusters and helps to define structural features. For the Niobrara project an 8X8 topology appeared to give the best stratigraphic resolution for instantaneous attributes and a 5X5 topology resolved the geometric attributes most effectively.

SOM Results for the Survey and their Interpretation

The SOM topology selected to best resolve the sub-Niobrara stratigraphy from the eight instantaneous attributes is an 8X8 hexagonal which yields 64 individual neurons. The SOM interval selected was Top Niobrara to Top Greenhorn. The next sequence of figures highlights the improved resolution provided by the SOM when compared to the original amplitude data. Figure 12 shows a north-south inline through the survey and through the Rotharmel 11-33 well which was one of the wells selected for petrophysical analysis. The original amplitude data is shown along with the SOM result for the interval.

Figure 12: North-South inline showing the original amplitude data (upper) and the 8X8 SOM result (lower) from Top Niobrara through Greenhorn horizons. Seismic data is shown courtesy of GPI and FFG.

The next image, Figure 13, zooms into the SOM and highlights the correlation with lithology from petrophysical analysis. The B chalk is noted by a stacked pattern of yellow-red-yellow neurons, with the red representing the maximum carbonate content within the middle of the chalk bench.

Figure 13: 8X8 Instantaneous SOM through Rotharmel 11-33 with well log composite. The B bench, highlighted in green on the wellbore, ties the yellow-red-yellow sequence of neurons. Seismic data is shown courtesy of GPI and FFG.

One can see on the SOM the sweet spot within the B chalk and that there is a fair amount of small-scale structural relief present. These results aid in the resolution of structural offset within the reservoir away from well control which is critical for staying in a 20 to 30 foot zone when drilling horizontally. Each classified sample is 1 ms in thickness which converted to depth equates to roughly 7 feet.

Figure 14 shows the K2 curvature attribute co-rendered with the SOM results in vertical sections. The Rotharmel 11-33 is at the intersection of the vertical sections. The curvature is extracted at the middle of the B chalk and shows good agreement with the SOM. The entire B bench is represented by only 5-6 ms of seismic data.

Figure 14: Most negative curvature, K2 rendered at the middle of the B chalk. Vertical sections are an 8X8 instantaneous SOM Top Niobrara to Top Greenhorn. Seismic data is shown courtesy of GPI and FFG.

A Marl Results

Seven wells within the survey were sent to a third party for petrophysical analysis (Figure 15). The analysis identified zones of interest within the Niobrara marls which are typically considered source rocks. The calculations show a high TOC zone in the upper A marl which the analysis identifies as shale pay (Figure 16). A seismic cross-section of the 8X8 instantaneous SOM (Figure 16) through the three wells depicted shows that this zone is well imaged. The neurons can be isolated and volumetric calculations derived from the representative neurons.

Figure 15: Index map for wells used in petrophysical analysis (in red)

Figure 16: Petrophysical results and SOM for three wells in the study area. The TOC curve (Track 12) and Shale pay curve (Track 10), highlighted in yellow, indicate the Upper A marl is both a rich source rock and a potential shale reservoir. Seismic data is shown courtesy of GPI and FFG.

Codell Results

The Codell sandstone in general and within the study area shows more heterogeneity in reservoir properties than the Niobrara chalk benches. The petrophysical analysis on 7 wells shows net pay ranging from zero feet to three feet. The gross thickness ranges from 17 feet to 25 feet. The SOM results reflect this heterogeneity, resolve the Codell gross interval throughout most of the study area, and thus, can be useful for horizontal well planning.

Figures 17 and 18 shows inline 60 through a well with the Top Niobrara to Greenhorn 8X8 SOM results. The 2D color map has been manipulated to emphasize the lower interval from approximately base Niobrara through the Codell. Figure 18 zooms into the well and shows the specific neurons associated with the Codell interval. Figures 19 shows a N-S traverse through four wells again with the Codell interval highlighted through use of a 2D color map. The western and southwest areas of the survey show a much more continuous character to the classification with only two neurons representing the Codell interval (6 and 48). Figure 20 shows both the N-S traverse and a crossline through the anomaly.

Figure 17: Instantaneous 8X8 SOM, Top Niobrara to Greenhorn. Seismic data is shown courtesy of GPI and FFG.

Figure 18: Detailed look at the Codell portion of the SOM at the Haythorn 4-12 with GR in background. The 2D color map shows how neurons can be isolated to show a specific stratigraphic interval. Seismic data is shown courtesy of GPI and FFG.

Figure 19: Traverse through 4 wells in the western part of the study area showing the isolation of the Codell sandstone within the SOM. The south west part of the line shows the Codell being represented by only 2 neurons (6 and 48). The colormap can be interrogated to determine which attributes contribute to any given neuron. Seismic data is shown courtesy of GPI and FFG.

Figure 20: View of the SW Codell anomaly where the neuron stacking pattern changes to two neurons only (6 and 47). Seismic data is shown courtesy of GPI and FFG.

Figure 21: 3D view of neurons isolated from the SOM in the Codell interval. The areas where red is prominent and continuous show the extent of Codell represented by neurons 6 and 47 only. Also, an area in the eastern part of the study is outlined. The Codell is not represented in this area by the six neurons highlighted in the 2D color map. Seismic data is shown courtesy of GPI and FFG.

Unfortunately, vertical well control was not available through this southwestern anomaly. To examine the extent of individual neurons within the SOM at Codell level, the next image, Figure 21, shows a 3D view of the isolated Codell neurons. The southwest anomaly is apparent as well as similar anomalies in the northern portion of the survey. What is also immediately apparent is that in the east-central portion of the survey, the Codell is not represented by the six neurons (6,7,47, 48, 55, 56) previously used to isolate it within the volume. Figure 22 takes a closer look at the SOM results through this area and also utilizes the original amplitude data. Both the SOM and the amplitude data show a change in character throughout the entire section, but the SOM results only change significantly in the lower Niobrara to Greenhorn portion of the interval.

The machine learning application has a feature in which individual neurons can be queried for statistics on how individual seismic attributes contribute to the cluster which makes up the neuron. Queries were done on all of the neurons within the Codell and shown are the results for neuron 6 which is one of 2 neurons characteristic of the southwestern Codell anomaly and on neuron 61in the area where the SOM changes significantly in Figure 23. Neuron 6 has equal contributions from Instantaneous Frequency, Hilbert, Thin Bed, and Relative Acoustic Impedance. Neuron 61 shows Instantaneous Q as the top attribute which is consistent with the interpretation of the section being structurally disturbed or highly fractured.

Figure 22: West-East crossline through two wells showing the SOM and amplitude data through the blank area from Figure 23. The seismic character and classification results differ significantly in this portion of the survey for the lower Niobrara, Fort Hays and Codell. This area is interpreted to be highly fractured. Seismic data is shown courtesy of GPI and FFG.

Figure 23: Example of attribute details for individual neurons (6 and 61). This shows the contribution of individual attributes to the neuron.

Structural Attributes

The machine learning workflow can be applied to geometric attributes. PCA and SOM need to be run separately from the instantaneous attributes since PCA assumes a Gaussian distribution of the attributes. This assumption doesn’t hold for geometric attributes but the SOM process does not assume any distribution and thus still finds patterns in the data. To produce a structural SOM, four attributes were selected from PCA: Curvature_K1, Similarity, Energy Ratio, Texture Entropy, and Texture Homogeneity. These were combined with the original amplitude data to generate SOMs from the Top Niobrara to Top Greenhorn interval. Several SOM topologies were generated with geometric attributes and a 5X5 yielded good results. Figure 24 shows the geometrical SOM results at the Top Niobrara, B bench, and Codell. The Top Niobrara level shows major faults, but not nearly as much structural disturbance as the mid-Niobrara B bench or the Codell level. The eastern part of the survey where the instantaneous classification changed also shows significant differences between the B bench and Codell and agrees with the interpretation that this is a highly fractured area for the lower Niobrara and Codell. The B bench appears more structurally disrupted than the Top Niobrara but shows fewer areal changes compared to Codell. Pressure and production data could help confirm how these features relate to reservoir quality.

Figure 24: 5X5 Structural SOM at 3 levels. There are significant changes both vertically and areally

Conclusions

Seismic multi-attribute analysis has always held the promise of improving interpretations via the integration of attributes which respond to subsurface conditions such as stratigraphy, lithology, faulting, fracturing, fluids, pressure, etc. Machine learning augments traditional interpretation and attribute analysis by utilizing attribute space to simultaneously classify suites of attributes into sample based, high dimension clusters that are subsequently visualized and further interpreted in the 3D seismic survey. 2D colormaps aid in their interpretation and visualization.

In the DJ Basin, we have resolved the primary reservoir targets, the Niobrara chalk benches and the Codell formation, represented within approximately 60 ms of data in two-way time, to the level of one to five neurons which is approximately 7 to 35 feet in thickness. Structural SOM classifications with a suite of geometric attributes better image the complex faulting and fracturing and its variations throughout the reservoir interval. The classification volumes are designed to aid in drilling target identification, reserves calculations and horizontal well planning.

Acknowledgements

The authors would like to thank their colleagues at Geophysical Insights for their valuable insight and suggestions and Digital Formation for the petrophysical analysis. We also thank Geophysical Pursuit, Inc. and Fairfield Geotechnologies for use of their data and permission to publish this paper.

References

Coleou, T., M. Poupon, and A. Kostia, 2003, Unsupervised seismic facies classification: A review and comparison
of techniques and implementation: The Leading Edge, 22, 942–953, doi: 10.1190/1.1623635.

Finn, T. M. and Johnson, R. C., 2005, Niobrara Total Petroleum System in the Southwestern Wyoming Province, Chapter 6 of Petroleum Systems and Geologic Assessment of Oil and Gas in the Southwestern Wyoming Province, Wyoming, Colorado, and Utah, USGS Southwestern Wyoming Province Assessment Team, U.S. Geological Survey Digital Data Series DDS–69–D.

Guo, H., K. J. Marfurt, and J. Liu, 2009, Principal component spectral analysis: Geophysics, 74, no. 4, p. 35–43.

Haykin, S., 2009, Neural networks and learning machines, 3rd ed.: Pearson.

Kauffman, E.G., 1977, Geological and biological overview— Western Interior Cretaceous Basin, in Kauffman, E.G., ed., Cretaceous facies, faunas, and paleoenvironments across the Western Interior Basin: The Mountain Geologist, v. 14, nos. 3 and 4, p. 75–99.

Kohonen, T., 2001, Self-organizing maps: Third extended addition, Springer, Series in Information Services, Vol. 30.

Landon, S.M., Longman, M.W., and Luneau, B.A., 2001, Hydrocarbon source rock potential of the Upper Cretaceous Niobrara Formation, Western Interior Seaway of the Rocky Mountain region: The Mountain Geologist, v. 38, no. 1, p. 1–18.

Lewis, R.K., 2013, Stratigraphy and Depositional Environments of the Late Cretaceous (Late Turonian) Codell Sandstone and Juana Lopez Member of the Carlile Shale, Southeast Colorado: Colorado School of Mines MS Thesis, 190 p.

Longman, M.W., Luneau, B.A., and Landon, S.M., 1998, Nature and distribution of Niobrara lithologies in the Cretaceous Western Interior Seaway of the Rocky Mountain Region: The Mountain Geologist, v. 35, no. 4, p. 137–170.

Luneau, B., Longman, M., Kaufman, P., and Landon, S., 2011, Stratigraphy and Petrophysical Characteristics of the Niobrara Formation in the Denver Basin, Colorado and Wyoming, AAPG Search and Discovery Article #50469.

Meissner, F.F., Woodward, J., and Clayton, J.L., 1984, Stratigraphic relationships and distribution of source rocks in the greater Rocky Mountain region, in Woodward, J., Meissner, F.F., and Clayton, J.L., eds., Hydrocarbon source rocks of the greater Rocky Mountain region: Rocky Mountain Association of Geologists Guidebook, p. 1–34.

Molenaar, C.M., and Rice, D.D., 1988, Cretaceous rocks of the Western Interior Basin, in Sloss, L.L., ed., Sedimentary cover-North American craton, U.S.: Geological Society of America, The Geology of North America, v. D–2, p. 77–82.

Roden, R., Smith, T., and Sacrey, D., 2015, Geologic pattern recognition from seismic attributes: Principal component analysis and self-organizing maps, Interpretation, Vol. 3, No. 4, p. SAE59-SAE83.

Santogrossi, P., 2017, Classification/Corroboration of Facies Architecture in the Eagle Ford Group: A Case Study in Thin Bed Resolution, URTeC 2696775, doi 10.15530-urtec-2017-<2696775>.

Smith, T., 2016, Why SOM is an Appealing Learning Machine, Internal Geophysical Insights Paper.

Sonnenberg, S.A., 2015. Geologic Factors Controlling Production in the Codell Sandstone, Wattenberg Field, Colorado. URTeC Paper 2145312 presented at the Unconventional Resources Technology Conference, San Antonio, TX, July 20-22.

Sonnenberg, S.A., 2015. New reserves in an old field, the Niobrara/Codell resource plays in the Wattenberg Field, Denver Basin, Colorado. EAGE First Break, v. 33, p. 55-62.

Sterling, R., Bottjer, R. and Smith, K., 2016, Codell SS, A review of the Northern DJ oil resource play Laramie County, WY and Weld, County, CO, AAPG Search and Discovery Article #10754.

Applications of Machine Learning for Geoscientists – Permian Basin

Applications of Machine Learning for Geoscientists – Permian Basin

By Carrie Laudon
Published with permission: Permian Basin Geophysical Society 60th Annual Exploration Meeting
May 2019

Abstract

Over the last few years, because of the increase in low-cost computer power, individuals and companies have stepped up investigations into the use of machine learning in many areas of E&P. For the geosciences, the emphasis has been in reservoir characterization, seismic data processing, and to a lesser extent interpretation. The benefits of using machine learning (whether supervised or unsupervised) have been demonstrated throughout the literature, and yet the technology is still not a standard workflow for most seismic interpreters. This lack of uptake can be attributed to several factors, including a lack of software tools, clear and well-defined case histories and training. Fortunately, all these factors are being mitigated as the technology matures. Rather than looking at machine learning as an adjunct to the traditional interpretation methodology, machine learning techniques should be considered the first step in the interpretation workflow.

By using statistical tools such as Principal Component Analysis (PCA) and Self Organizing Maps (SOM) a multi-attribute 3D seismic volume can be “classified”. The PCA reduces a large set of seismic attributes both instantaneous and geometric, to those that are the most meaningful. The output of the PCA serves as the input to the SOM, a form of unsupervised neural network, which, when combined with a 2D color map facilitates the identification of clustering within the data volume. When the correct “recipe” is selected, the clustered or classified volume allows the interpreter to view and separate geological and geophysical features that are not observable in traditional seismic amplitude volumes. Seismic facies, detailed stratigraphy, direct hydrocarbon indicators, faulting trends, and thin beds are all features that can be enhanced by using a classified volume.

The tuning-bed thickness or vertical resolution of seismic data traditionally is based on the frequency content of the data and the associated wavelet. Seismic interpretation of thin beds routinely involves estimation of tuning thickness and the subsequent scaling of amplitude or inversion information below tuning. These traditional below-tuning-thickness estimation approaches have limitations and require assumptions that limit accuracy. The below tuning effects are a result of the interference of wavelets, which are a function of the geology as it changes vertically and laterally. However, numerous instantaneous attributes exhibit effects at and below tuning, but these are seldom incorporated in thin-bed analyses. A seismic multi-attribute approach employs self-organizing maps to identify natural clusters from combinations of attributes that exhibit below-tuning effects. These results may exhibit changes as thin as a single sample interval in thickness. Self-organizing maps employed in this fashion analyze associated seismic attributes on a sample-by-sample basis and identify the natural patterns or clusters produced by thin beds. Examples of this approach to improve stratigraphic resolution in both the Eagle Ford play, and the Niobrara reservoir of the Denver-Julesburg Basin will be used to illustrate the workflow.

Introduction

Seismic multi-attribute analysis has always held the promise of improving interpretations via the integration of attributes which respond to subsurface conditions such as stratigraphy, lithology, faulting, fracturing, fluids, pressure, etc. The benefits of using machine learning (whether supervised or unsupervised) has been demonstrated throughout the literature and yet the technology is still not a standard workflow for most seismic interpreters. This lack of uptake can be attributed to several factors, including a lack of software tools, clear and well-defined case histories, and training. This paper focuses on an unsupervised machine learning workflow utilizing Self-Organizing Maps (Kohonen, 2001) in combination with Principal Component Analysis to produce classified seismic volumes from multiple instantaneous attribute volumes. The workflow addresses several significant issues in seismic interpretation: it analyzes large amounts of data simultaneously; it determines relationships between different types of data; it is sample based and produces high-resolution results and, reveals geologic features that are difficult to see in conventional approaches.

Principal Component Analysis (PCA)

Multi-dimensional analysis and multi-attribute analysis go hand in hand. Because individuals are grounded in three-dimensional space, it is difficult to visualize what data in a higher number dimensional space looks like. Fortunately, mathematics doesn’t have this limitation and the results can be easily understood with conventional 2D and 3D viewers.

Working with multiple instantaneous or geometric seismic attributes generates tremendous volumes of data. These volumes contain huge numbers of data points which may be highly continuous, greatly redundant, and/or noisy. (Coleou et al., 2003). Principal Component Analysis (PCA) is a linear technique for data reduction which maintains the variation associated with the larger data sets (Guo and others, 2009; Haykin, 2009; Roden and others, 2015). PCA can separate attribute types by frequency, distribution, and even character. PCA technology is used to determine which attributes may be ignored due to their very low impact on neural network solutions and which attributes are most prominent in the data. Figure 1 illustrates the analysis of a data cluster in two directions, offset by 90 degrees. The first principal component (eigenvector 1) analyses the data cluster along the longest axis. The second principal component (eigenvector 2) analyses the data cluster variations perpendicular to the first principal component. As stated in the diagram, each eigenvector is associated with an eigenvalue which shows how much variance there is in the data.

two attribute data set

Figure 1. Two attribute data set illustrating the concept of PCA

The next step in PCA analysis is to review the eigen spectrum to select the most prominent attributes in a data set. The following example is taken from a suite of instantaneous attributes over the Niobrara formation within the Denver­ Julesburg Basin. Results for eigenvectors 1 are shown with three attributes: sweetness, envelope and relative acoustic impedance being the most prominent.

two attribute data set

Figure 2. Results from PCA for first eigenvector in a seismic attribute data set

Utilizing a cutoff of 60% in this example, attributes were selected from PCA for input to the neural network classification. For the Niobrara, eight instantaneous attributes from the four of the first six eigenvectors were chosen and are shown in Table 1. The PCA allowed identification of the most significant attributes from an initial group of 19 attributes.

Results from PCA for Niobrara Interval

Table 1: Results from PCA for Niobrara Interval shows which instantaneous attributes will be used in a Self-Organizing Map (SOM).

Self-Organizing Maps

Teuvo Kohonen, a Finnish mathematician, invented the concepts of Self-Organizing Maps (SOM) in 1982 (Kohonen, T., 2001). Self-Organizing Maps employ the use of unsupervised neural networks to reduce very high dimensions of data to a classification volume that can be easily visualized (Roden and others, 2015). Another important aspect of SOMs is that every seismic sample is used as input to classification as opposed to wavelet-based classification.

Figure 3 diagrams the SOM concept for 10 attributes derived from a 3D seismic amplitude volume. Within the 3D seismic survey, samples are first organized into attribute points with similar properties called natural clusters in attribute space. Within each cluster new, empty, multi-attribute samples, named neurons, are introduced. The SOM neurons will seek out natural clusters of like characteristics in the seismic data and produce a 2D mesh that can be illustrated with a two- dimensional color map. In other words, the neurons “learn” the characteristics of a data cluster through an iterative process (epochs) of cooperative than competitive training. When the learning is completed each unique cluster is assigned to a neuron number and each seismic sample is now classified (Smith, 2016).

two attribute data set

Figure 3. Illustration of the concept of a Self-Organizing Map

Figures 4 and 5 show a simple example using 2 attributes, amplitude, and Hilbert transform on a synthetic example. Synthetic reflection coefficients are convolved with a simple wavelet, 100 traces created, and noise added. When the attributes are cross plotted, clusters of points can be seen in the cross plot. The colored cross plot shows the attributes after SOM classification into 64 neurons with random colors assigned. In Figure 5, the individual clusters are identified and mapped back to the events on the synthetic. The SOM has correctly distinguished each event in the synthetic.

Two attribute synthetic example of a Self-Organizing Map

Figure 4. Two attribute synthetic example of a Self-Organizing Map. The amplitude and Hilbert transform are cross plotted. The colored cross plot shows the attributes after classification into 64 neurons by SOM.

Synthetic SOM example

Figure 5. Synthetic SOM example with neurons identified by number and mapped back to the original synthetic data

Results for Niobrara and Eagle Ford

In 2018, Geophysical Insights conducted a proof of concept on 100 square miles of multi-client 3D data jointly owned by Geophysical Pursuit, Inc. (GPI) and Fairfield Geotechnologies (FFG) in the Denver¬ Julesburg Basin (DJ). The purpose of the study is to evaluate the effectiveness of a machine learning workflow to improve resolution within the reservoir intervals of the Niobrara and Codell formations, the primary targets for development in this portion of the basin. An amplitude volume was resampled from 2 ms to 1 ms and along with horizons, loaded into the Paradise® machine learning application and attributes generated. PCA was used to identify which attributes were most significant in the data, and these were used in a SOM to evaluate the interval Top Niobrara to Greenhorn (Laudon and others, 2019).

Figure 6 shows results of an 8X8 SOM classification of 8 instantaneous attributes over the Niobrara interval along with the original amplitude data. Figure 7 is the same results with a well composite focused on the B chalk, the best section of the reservoir, which is difficult to resolve with individual seismic attributes. The SOM classification has resolved the chalk bench as well as other stratigraphic features within the interval.

North-South Inline showing the original amplitude data (upper) and the 8X8 SOM result (lower) from Top Niobrara

Figure 6. North-South Inline showing the original amplitude data (upper) and the 8X8 SOM result (lower) from Top Niobrara through Greenhorn horizons. Seismic data is shown courtesy of GPI and FFG.

8X8 Instantaneous SOM through Rotharmel 11-33 with well log composite

Figure 7. 8X8 Instantaneous SOM through Rotharmel 11-33 with well log composite. The B bench, highlighted in green on the wellbore, ties the yellow-red-yellow sequence of neurons. Seismic data is shown courtesy of GPI and FFG

 

8X8 SOM results through the Eagle Ford

Figure 8. 8X8 SOM results through the Eagle Ford. The primary target, the Lower Eagle Ford shale had 16 neuron classes over 14-29 milliseconds of data. Seismic data shown courtesy of Seitel.

The results shown in Figure 9 reveal non-layer cake facies bands that include details in the Eagle )RUG,v basal clay-rich shale, high resistivity and low resistivity Eagle Ford shale objectives, the Eagle Ford ash, and the upper Eagle Ford marl, which are overlain disconformably by the Austin Chalk.

Eagle Ford SOM classification shown with well results

Figure 9. Eagle Ford SOM classification shown with well results. The SOM resolves a high resistivity interval, overlain by a thin ash layer and finally a low resistivity layer. The SOM also resolves complex 3-dimensional relationships between these facies

Convolutional Neural Networks (CNN)

A promising development in machine learning is supervised classification via the applications of convolutional neural networks (CNNs). Supervised methods have, in the past, not been efficient due to the laborious task of training the neural network. CNN is a deep learning seismic classification. We apply CNN to fault detection on seismic data. The examples that follow show CNN fault detection results which did not require any interpreter picked faults for training, rather the network was trained using synthetic data. Two results are shown, one from the North Sea, Figure 10, and one from the Great South Basin, New Zealand, Figure 11.

Side by side comparison of coherence attribute to CNN fault probability attribute, North Sea

Figure 10. Side by side comparison of coherence attribute to CNN fault probability attribute, North Sea

Side by side comparison of coherence attribute to CNN fault probability attribute, North Sea

Figure 11. Comparison of Coherence to CNN fault probability attribute, New Zealand

Conclusions

Advances in compute power and algorithms are making the use of machine learning available on the desktop to seismic interpreters to augment their interpretation workflow. Taking advantage of today’s computing technology, visualization techniques, and an understanding of machine learning as applied to seismic data, PCA combined with SOMs efficiently distill multiple seismic attributes into classification volumes. When applied on a multi-attribute seismic sample basis, SOM is a powerful nonlinear cluster analysis and pattern recognition machine learning approach that helps interpreters identify geologic patterns in the data and has been able to reveal stratigraphy well below conventional tuning thickness.

In the fault interpretation domain, recent development of a Convolutional Neural Network that works directly on amplitude data shows promise to efficiently create fault probability volumes without the requirement of a labor-intensive training effort.

References

Coleou, T., M. Poupon, and A. Kostia, 2003, Unsupervised seismic facies classification: A review and comparison of techniques and implementation: The Leading Edge, 22, 942–953, doi: 10.1190/1.1623635.

Guo, H., K. J. Marfurt, and J. Liu, 2009, Principal component spectral analysis: Geophysics, 74, no. 4, 35–43.

Haykin, S., 2009. Neural networks and learning machines, 3rd ed.: Pearson

Kohonen, T., 2001,Self organizing maps: Third extended addition, Springer, Series in Information Services, Vol. 30.

Laudon, C., Stanley, S., and Santogrossi, P., 2019, Machine Leaming Applied to 3D Seismic Data from the Denver-Julesburg Basin Improves Stratigraphic Resolution in the Niobrara, URTeC 337, in press

Roden, R., and Santogrossi, P., 2017, Significant Advancements in Seismic Reservoir Characterization with Machine Learning, The First, v. 3, p. 14-19

Roden, R., Smith, T., and Sacrey, D., 2015, Geologic pattern recognition from seismic attributes: Principal component analysis and self-organizing maps, Interpretation, Vol. 3, No. 4, p. SAE59-SAE83.

Santogrossi, P., 2017, Classification/Corroboration of Facies Architecture in the Eagle Ford Group: A Case Study in Thin Bed Resolution, URTeC 2696775, doi 10.15530-urtec-2017-<2696775>.

Approach Aids Multiattribute Analysis

Approach Aids Multiattribute Analysis

By: Rocky Roden, Geophysical Insights, and Deborah Sacrey, Auburn Energy
Published with permission: American Oil and Gas Reporter
September 2015

Seismic attributes, which are any measurable properties of seismic data, aid interpreters in identifying geologic features that are not understood clearly in the original data. However, the enormous amount of information generated from seismic attributes and the difficulty in understanding how these attributes when combined define geology, requires another approach in the interpretation workflow.

To address these issues, “machine learning” to evaluate seismic attributes has evolved over the last few years. Machine learning uses computer algorithms that learn iteratively from the data and adapt independently to produce reliable, repeatable results. Applying current computing technology and visualization techniques, machine learning addresses two significant issues in seismic interpretation:

• The big data problem of trying to interpret dozens, if not hundreds, of volumes of data; and

• The fact that humans cannot understand the relationship of several types of data all at once.

Principal component analysis (PCA) and self-organizing maps (SOMs) are machine learning approaches that when applied to seismic multiattribute analysis are producing results that reveal geologic features not previously identified or easily interpreted. Applying principal component analysis can help interpreters identify seismic attributes that show the most variance in the data for a given geologic setting, which helps determine which attributes to use in a multiattribute analysis using self-organizing maps. SOM analysis enables interpreters to identify the natural organizational patterns in the data from multiple seismic attributes.

Multiple-attribute analyses are beneficial when single attributes are indistinct. These natural patterns or clusters represent geologic information embedded in the data and can help identify geologic features, geobodies, and aspects of geology that often cannot be interpreted by any other means. SOM evaluations have proven to be beneficial in essentially all geologic settings, including unconventional resource plays, moderately compacted onshore regions, and offshore unconsolidated sediments.

This indicates the appropriate seismic attributes to employ in any SOM evaluation should be based on the interpretation problem to be solved and the associated geologic setting. Applying PCA and SOM can not only identify geologic patterns not seen previously in the seismic data, it also can increase or decrease confidence in features already interpreted. In other words, this multiattribute approach provides a methodology to produce a more accurate risk assessment of a geoscientist’s interpretation and may represent the next generation of advanced interpretation.

Seismic Attributes

A seismic attribute can be defined as any measure of the data that helps to visually enhance or quantify features of interpretation interest. There are hundreds of types of attributes, but Table 1 shows a composite list of seismic attributes and associated categories routinely employed in seismic interpretation. Interpreters wrestle continuously with evaluating the numerous seismic attribute volumes, including visually co-blending two or three attributes and even generating attributes from other attributes in an effort to better interpret their data.

This is where machine learning approaches such as PCA and SOM can help interpreters evaluate their data more efficiently, and help them understand the relationships between numerous seismic attributes to produce more accurate results.

Principal Component Analysis

Principal component analysis is a linear mathematical technique for reducing a large set of seismic attributes to a small set that still contains most of the variation in the large set. In other words, PCA is a good approach for identifying the combination of the most meaningful seismic attributes generated from an original volume.

principal component analysis for multiattribute analysis

Results from Principal Component Analysis in Paradise® utilizing 18 instantaneous seismic attributes are shown here. 1A shows histograms of the highest eigenvalues for in-lines in the seismic 3-D volume, with red histograms representing eigenvalues over the field. 1B shows the average of eigenvalues over the field (red), with the first principal component in orange and associated seismic attribute contributions to the right. 1C shows the second principal component over the field with the seismic attribute contributions to the right. The top five attributes in 1B were run in SOM A and the top four attributes in 1C were run in SOM B.

The first principal component accounts for as much of the variability in the data as possible, and each succeeding component (orthogonal to each preceding component) accounts for as much of the remaining variability. Given a set of seismic attributes generated from the same original volume, PCA can identify the attributes producing the largest variability in the data, suggesting these combinations of attributes will better identify specific geologic features of interest.

Even though the first principal component represents the largest linear attribute combinations best representing the variability of the bulk of the data, it may not identify specific features of interest. The interpreter should evaluate succeeding principal components also because they may be associated with other important aspects of the data and geologic features not identified with the first principal component.

In other words, PCA is a tool that, when employed in an interpretation workflow, can give direction to meaningful seismic attributes and improve interpretation results. It is logical, therefore, that a PCA evaluation may provide important information on appropriate seismic attributes to take into generating a self-organizing map.

Self-Organizing Maps

The next level of interpretation requires pattern recognition and classification of the often subtle information embedded in the seismic attributes. Taking advantage of today’s computing technology, visualization techniques and understanding of appropriate parameters, self-organizing maps distill multiple seismic attributes efficiently into classification and probability volumes. SOM is a powerful non- linear cluster analysis and pattern recognition approach that helps interpreters identify patterns in their data that can relate to desired geologic characteristics such as those listed in Table 1.

Seismic data contain huge amounts of data samples and are highly continuous, greatly redundant and significantly noisy. The tremendous amount of samples from numerous seismic attributes exhibit significant organizational structure in the midst of noise. SOM analysis identifies these natural organizational structures in the form of clusters. These clusters reveal significant information about the classification structure of natural groups that is difficult to view any other way. The natural groups and patterns in the data identified by clusters reveal the geology and aspects of the data that are difficult to interpret otherwise.

Offshore Case Study

Offshore Case Study 01

This shows SOM A results from Paradise on a north-south inline through the field. 1A shows the original stacked amplitude. 2B shows SOM results with the associated five-by-five color map displaying all 25 neurons. 2C shows SOM results with four neurons elected that isolate attenuation effects.

Offshore Case Study 02

SOM B results from Paradise are shown on the same in-line as Figure 2. 3A is the original stacked amplitude. 3B shows SOM results with the associated five-by-five color map. 3C is the SOM results with a color map showing two neurons that highlight flat spots in the data.

 

A case study is provided by a lease located in the Gulf of Mexico offshore Louisiana in 470 feet of water. This shallow field (approximately 3,900 feet) has two producing wells that were drilled on the upthrown side of an east-west trending normal fault and into an amplitude anomaly identified on the available 3-D seismic data. The normally pressured reservoir is approximately 100 feet thick and is located in a typical “bright spot” setting, i.e. a Class 3 AVO geologic setting (Rutherford and Williams, 1989).

The goal of this multiattribute analysis is to more clearly identify possible direct hydrocarbon indicator characteristics such as flat spots (hydrocarbon contacts) and attenuation effects and to better understand the reservoir and provide important approaches for decreasing the risk of future exploration in the area.

Initially, 18 instantaneous seismic attributes were generated from the 3-D data in the area. These were put into a PCA evaluation to determine which produced the largest variation in the data and the most meaningful attributes for SOM analysis.

The PCA was computed in a window 20 milliseconds above and 150 milliseconds below the mapped top of the reservoir over the entire survey, which encompassed approximately 10 square miles. Each bar in Figure 1A represents the highest eigenvalue on its associated in-line over the portion of the survey displayed.

An eigenvalue shows how much variance there is in its associated eigenvector, and an eigenvector is a direction showing the spread in the data. The red bars in Figure 1A specifically denote the in-lines that cover the areal extent of the amplitude feature, and the average of their eigenvalue results are displayed in Figures 1B and 1C.

Figure 1B displays the principal components from the selected in-lines over the anomalous feature with the highest eigenvalue (first principal component), indicating the percentage of seismic attributes contributing to this largest variation in the data. In this first principal component, the top seismic attributes include trace envelope, envelope modulated phase, envelope second derivative, sweetness and average energy, all of which account for more than 63 percent of the variance of all the instantaneous attributes in this PCA evaluation.

Figure 1C displays the PCA results, but this time the second highest eigenvalue was selected and produced a different set of seismic attributes. The top seismic attributes from the second principal component include instantaneous frequency, thin bed indicator, acceleration of phase, and dominant frequency, which total almost 70 percent of the variance of the 18 instantaneous seismic attributes analyzed. These results suggest that when applied to a SOM analysis, perhaps the two sets of seismic attributes for the first and second principal components will help define different types of anomalous features or different characteristics of the same feature.

The first SOM analysis (SOM A) incorporates the seismic attributes defined by the PCA with the highest variation in the data, i.e., the five highest percentage contributing attributes in Figure 1B.

Several neuron counts for SOM analyses were run on the data, and lower count matrices revealed broad, discrete features, while the higher counts displayed more detail and less variation. The SOM results from a five-by-five matrix of neurons (25) were selected for this article.

 

Detecting Attenuation

The north-south line through the field in Figures 2 and 3 show the original stacked amplitude data and classification results from the SOM analyses. In Figure 2B, the color map associated with the SOM classification results indicates all 25 neurons are displayed. Figure 2C shows results with four interpreted neurons highlighted.

Based on the location of the hydrocarbons determined from well control, it is interpreted from the SOM results that attenuation in the reservoir is very pronounced. As Figures 2B and 2C reveal, there is apparent absorption banding in the reservoir above the known hydrocarbon contacts defined by the wells in the field. This makes sense because the seismic attributes employed are sensitive to relatively low-frequency, broad variations in the seismic signal often associated with attenuation effects.

This combination of seismic attributes employed in the SOM analysis generates a more pronounced and clearer picture of attenuation in the reservoir than any of the seismic attributes or the original amplitude volume individually. Downdip of the field is another undrilled anomaly that also reveals apparent attenuation effects.

The second SOM evaluation (SOM B) includes the seismic attributes with the highest percentages from the second principal component, based on the PCA (see Figure 1). It is important to note that these attributes are different from the attributes determined from the first principal component. With a five-by-five neuron matrix, Figure 3 shows the classification results from this SOM evaluation on the same north-south line as Figure 2, and it identifies clearly several hydrocarbon contacts in the form of flat spots. These hydrocarbon contacts are confirmed by the well control.

Figure 3B defines three apparent flat spots that are further isolated in Figure 3C, which displays these features with two neurons. The gas/oil contact in the field was very difficult to see in the original seismic data, but is well defined and can be mapped from this SOM analysis.

The oil/water contact in the field is represented by a flat spot that defines the overall base of the hydrocarbon reservoir. Hints of this oil/water contact were interpreted from the original amplitude data, but the second SOM classification provides important information to clearly define the areal extent of reservoir.

Downdip of the field is another apparent flat spot event that is undrilled and is similar to the flat spots identified in the field. Based on SOM evaluations A and B in the field, which reveal similar known attenuation and flat spot results, respectively, there is a high probability this undrilled feature contains hydrocarbons.

West Texas Case Study

Unlike the Gulf of Mexico case study, attribute analyses on the Fasken Ranch in the Permian Basin involved using a “recipe” of seismic attributes, based on their ability to sort out fluid properties, porosity trends and hydrocarbon sensitivities. Rather than use principal component analysis to see which attributes had the greatest variation in the data, targeted use of specific attributes helped solve an issue regarding conventional porosity zones within an unconventional depositional environment in the Spraberry and Wolfcamp formations.

The Fasken Ranch is located in portions of Andrews, Ector, Martin and Midland counties, Tx. The approximately 165,000-acre property, which consists of surface and mineral rights, is held privately. This case study shows the SOM analysis results for one well, the Fasken Oil and Ranch No. 303 FEE BI, which was drilled as a straight hole to a depth of 11,195 feet. The well was drilled through the Spraberry and Wolfcamp formations and encountered a porosity zone from 8,245 to 8,270 feet measured depth.

This enabled the well to produce more than four times the normal cumulative production found in a typical vertical Spraberry well. The problem was being able to find that zone using conventional attribute analysis in the seismic data. Figure 4A depicts cross-line 516, which trends north-south and shows the intersection with well 303. The porosity zone is highlighted with a red circle.

water oil contact

4A is bandwidth extension amplitude volume, highlighting the No. 303 well and porosity zone. Wiggle trace overlay is from amplitude volume. 4B is SOM classification volume, highlighting the No. 303 well and porosity zone. Topology was 10-by-10 neurons with a 30-millisecond window above and below the zone of interest. Wiggle trace overlay is from amplitude volume.

Seven attributes were used in the neural analysis: attenuation, BE14-100 (amplitude volume), average energy, envelope time derivative, density (derived through prestack inversion), spectral decomposition envelop sub-band at 67.3 hertz, and sweetness.

Figure 4B is the same cross-line 516, showing the results of classifying the seven attributes referenced. The red ellipse shows the pattern in the data that best represents the actual porosity zone encountered in the well, but could not be identified readily by conventional attribute analysis.

Figure 5 is a 3-D view of the cluster of neurons that best represent porosity. The ability to isolate specific neurons enables one to more easily visualize specific stratigraphic events in the data.

neural cluster with colormap

This SOM classification volume in 3-D view shows the combination of a neural “cluster” that represents the porosity zone seen in the No. 303 well, but not seen in surrounding wells.

 

 

Conclusions

Seismic attributes help identify numerous geologic features in conventional seismic data. Applying principal component analysis can help interpreters identify seismic attributes that show the most variance in the data for a given geologic setting, and help them determine which attributes to use in a multiattribute analysis using self-organizing maps. Applying current computing technology, visualization techniques, and understanding of appropriate parameters for SOM enables interpreters to take multiple seismic attributes and identify the natural organizational patterns in the data.

Multiple-attribute analyses are beneficial when single attributes are indistinct. These natural patterns or clusters represent geologic information embedded in the data and can help identify geologic features that often cannot be interpreted by any other means. Applying SOM to bring out geologic features and anomalies of significance may indicate this approach represents the next generation of advanced interpretation.

 

Editor’s Note

The authors wish to thank the staff of Geophysical Insights for researching and developing the applications used in this article. The seismic data for the Gulf of Mexico case study is courtesy of Petroleum Geo-Services. Thanks to T. Englehart for insight into the Gulf of Mexico case study. The authors also would like to acknowledge Glenn Winters and Dexter Harmon of Fasken Ranch for the use of the Midland Merge 3-D seismic survey in the West Texas case study.

Rocky Roden ROCKY RODEN runs his own consulting company, Rocky Ridge Resources Inc., and works with oil companies around the world on interpretation technical issues, prospect generation, risk analysis evaluations, and reserve/resource calculations. He is a senior consulting geophysicist with Houston-based Geophysical
Insights, helping develop advanced geophysical technology for interpretation.
He also is a principal in the Rose and Associates DHI Risk Analysis Consortium,
which is developing a seismic amplitude risk analysis program and worldwide
prospect database. Roden also has worked with Seismic Microtechnology
and Rock Solid Images on integrating advanced geophysical software applications.
He holds a B.S. in oceanographic technology-geology from Lamar University
and a M.S. in geological and geophysical oceanography from Texas A&M University.
Deborah Sacrey DEBORAH KING SACREY is a geologist/geophysicist with 39 years of oil and gas exploration experience in the Texas and Louisiana Gulf Coast, and Mid-Continent areas. For the past three years, she has been part of a Geophysical Insights team working to bring the power of multiattribute neural analysis of seismic data to the geoscience public. Sacrey received a degree in geology from the University of Oklahoma in 1976, and immediately started working for Gulf Oil. She started her own company, Auburn Energy, in 1990, and built her first geophysical workstation using
Kingdom software in 1995. She specializes in 2-D and 3-D interpretation
for clients in the United States and internationally. Sacrey is a DPA certified
petroleum geologist and DPA certified petroleum geophysicist.