It’s Just a Phase (Residue)

Contributors: Oswaldo Davogustto

In the March issue of Geophysical Corner, my colleagues Marcilio Matos and AAPG member Kurt Marfurt discussed the concept of phase unwrapping and the computation of phase residues.

This month, we go deeper: I elaborate on the application of phase residues to seismic data – and the resulting subsequent interpretation of geological features.

Some geologically induced spatial discontinuities, such as channels and faults, easily can be identified as seismic phase shifts or amplitude anomalies when they are above seismic resolution – but phase shifts from condensed sections and erosional unconformities can be subtle and not as easily detected.

Spectral decomposition is a proven, powerful means of identifying such discontinuities at specific frequencies that are otherwise buried in the seismic broadband response.

Although seismic acquisition and processing preserve seismic phase very well, little has been published about interpreting the phase components resulting from spectral decomposition.

Morlet complex wavelet transform phase residues can improve seismic spectral decomposition interpretation by detecting the phase discontinuities in the joint time-frequency spectral phase component – by evaluating the phase shifts that are derived from thickness changes in a wedge model.

We unwrap phase the phase traversing a rectangular contour about each time-frequency sample. In almost all incidents, the contour closes. However, in some cases we have a +180 or -180 degree phase anomaly. We display the location of these phase residue anomalies and correlate them to stratigraphic discontinuities and inconsistencies in seismic data quality.

We are able to map interference patterns between the wavelets that occur below seismic resolution.

For example, here we apply the phase residues to a seismic dataset that served as one of the first published applications of spectral decomposition.

The geology consists of an incised valley system in the Red Fork formation of the Anadarko Basin that had at least five stages of incision and fill. Previous works described that the fill of the mentioned incised valley as comprising lag deposits, shales, coals, muddy sands and sands.

This dataset also is known for the occurrence of invisible channels – channels that can be detected and correlated through logs, but that are below the resolution of the seismic data.

Based on the response of the phase residues and the tops interpreted from the logs, we identified the incision stages in the phase residues attribute. We interpreted each stage as a seismic horizon and constructed a geological model that honors the well, seismic and attribute data.

This geological model can be further used in reservoir modeling for reservoir properties, such as net-to-gross, porosity and permeability, with the versatility that these properties can change as the depositional environment changed from stage to stage.

Figure 1 shows our well-to-seismic calibration for two wells, A and B. Locations of wells A and B are shown in figure 2c.

The correlation coefficient is 75 percent for both wells.

We identify very distinctive patterns in the log response for each well – well A is located in the regional Red Fork facies and shows faster P-wave velocity, whereas well B is located in the incised valley system facies and displays lower P-wave velocities from the sonic log.

In figure 2 we display the seismic data (a), the phase residues response of a representative time section (b) and a time slice through seismic amplitude at 1.8 s (c).

Seismic data (a) are able to resolve only one of the incision stages identified on the time slice of the data (c). This is a common problem in midcontinent datasets, where channels are identified in time slices but not on vertical sections.

Using phase residues (b) we identify three anomaly responses that correlate with the channel-like features interpreted as incised valleys in the time slice data (c).

Using the well and the phase residues data, we interpreted the incision stages.

In figure 3 we show a chair diagram of the seismic data with a time slice at 1.8 s (a) and a 3-D view of the geo-cellular grid constructed from the combined well log and phase residues interpretation (b). This detailed geo-cellular grid allows us to model the properties of each incision stage and the regional Red Fork as a separate event with their on-reservoir properties and modeling technique.

In conclusion, we demonstrated how the use of phase residues can be effectively applied to reveal and enhance important stratigraphic features not otherwise revealed by conventional seismic amplitude.

We have developed a workflow that combines well data with phase attributes in order to produce a well-to-seismic consistent stratigraphic model.

I would like to thank Mark Falk and Al Warner for their support and advice in this project. I also would like to thank Chesapeake Energy Corporation and CGG-Veritas for donating the data for this project, and to Schlumberger and CGG-Veritas for facilitating the software used in these displays.

Comments (0)


Geophysical Corner

Geophysical Corner - Satinder Chopra
Satinder Chopra, award-winning chief geophysicist (reservoir), at Arcis Seismic Solutions, Calgary, Canada, and a past AAPG-SEG Joint Distinguished Lecturer began serving as the editor of the Geophysical Corner column in 2012.

Oswaldo Davogustto is an AAPG Student member at 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

Desktop /Portals/0/images/_site/AAPG-newlogo-vertical-morepadding.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 4046 Book
Desktop /Portals/0/PackFlashItemImages/WebReady/book-s64-Heavy-oil-and-Oil-sand-Petroleum-Systems.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 4062 Book

See Also: Bulletin Article

Understanding the distribution and geometry of reservoir geobodies is crucial for net-to-gross estimates and to model subsurface flow. This article focuses on the process of dolomitization and resulting geometry of diagenetic geobodies in an outcrop of Jurassic host rocks from northern Oman. Field and petrographic data show that a first phase of stratabound dolomite is crosscut by a second phase of fault-related dolomite. The stratabound dolomite geobodies are laterally continuous for at least several hundreds of meters (sim1000 ft) and probably regionally and are one-half meter (1.6 ft) thick. Based on petrography and geochemistry, a process of seepage reflux of mesosaline or hypersaline fluids during the early stages of burial diagenesis is proposed for the formation of the stratabound dolomite. In contrast, the fault-related dolomite geobodies are trending along a fault that can be followed for at least 100 m (328 ft) and vary in width from a few tens of centimeters to as much as 10 m (sim1–33 ft). Petrography, geochemistry, and high homogenization temperature of fluid inclusions all point to the formation of the dolomite along a normal fault under deep burial conditions during the Middle to Late Cretaceous. The high 87Sr/86Sr ratio in the dolomite and the high salinity measured in fluid inclusions indicate that the dolomitizing fluids are deep basinal brines that interacted with crystalline basement. The dolomitization styles have an impact on the dimension, texture, and geochemistry of the different dolomite geobodies, and a modified classification scheme (compared to the one from Jung and Aigner, 2012) is proposed to incorporate diagenetic geobodies in future reservoir modeling.
Desktop /Portals/0/PackFlashItemImages/WebReady/linking-process-dimension-texture-and-geochemistry.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3772 Bulletin Article

See Also: Energy Policy Blog

As Arctic ice thins and retreats, growing commercial access–from cruise ships to oil and gas drilling–will increase the risk of an oil spill. A recent National Research Council report found that resources, technology, research, manpower, funding and logistics are inadequate to respond to an Arctic oil spill.

Desktop /Portals/0/PackFlashItemImages/WebReady/National-Academies-US-Inadequately-Prepared-for-an-Arctic-Oil-Spill-2014-04apr-30-hero.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12896 Energy Policy Blog

Witnesses stated that BLM policies create unnecessary obstacles to oil and gas activities on federal lands at a recent House Natural Resources, Energy and Mineral Resources subcommittee hearing: “Energy in America: BLM’s Red-Tape Run Around and its Impact on American Energy Production”

Desktop /Portals/0/PackFlashItemImages/WebReady/oil-and-gas-drilling-on-federal-lands-2014-02feb-24.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12376 Energy Policy Blog