‘Hilbert Transform’ Remains a Valuable Tool

Geological interpretation of seismic data is commonly done by analyzing patterns of seismic amplitude, phase and frequency in map and section views across a prospect area. Although many seismic attributes have been utilized to emphasize geologic targets and to define critical rock and fluid properties, these three simple attributes – amplitude, phase and frequency – remain the mainstay of geological interpretation of seismic data.

Any procedure that extracts and displays any of these seismic parameters in a convenient and understandable manner is an invaluable interpretation tool.

A little more than 30 years ago, M.T. Taner and Robert E. Sheriff introduced the concept of using the Hilbert transform to calculate seismic amplitude, phase and frequency instantaneously – meaning a value for each parameter is calculated at each time sample of a seismic trace.

That Hilbert transform approach now forms the basis by which almost all amplitude, phase and frequency attributes are calculated by today’s seismic interpretation software

The Complex Seismic Trace

The action of the Hilbert transform is to convert a seismic trace x(t) into what first appears to be a mysterious complex seismic trace z(t) as shown on figure 1.

In this context, the term “complex” is used in its mathematical sense, meaning it refers to a number that has a real part and an imaginary part. The term does not imply that the data are difficult to understand.

This complex trace consists of the real seismic trace x(t) and an imaginary seismic trace y(t) that is the Hilbert transform of x(t).

On figure 1 these two traces are shown in a three-dimensional data space (x, y, t), where t is seismic time, x is the real-data plane, and y is the imaginary-data plane. The actual seismic trace is confined to the real-data plane; the Hilbert transform trace is restricted to the imaginary-data plane.

These two traces combine to form a complex trace z(t), which appears as a helix that spirals around the time axis.

The projection of complex trace z(t) onto the real plane is the actual seismic trace x(t); the projection of z(t) onto the imaginary plane is the Hilbert transform trace y(t).

At any coordinate on the time axis, a vector a(t) can be calculated that extends perpendicularly away from the time axis to intercept the helical complex trace z(t) as shown on figure 2. The length of this vector is the amplitude of the complex trace at that particular instant in time – hence the term “instantaneous amplitude.”

The amplitude value is calculated using the equation for a(t) shown on the figure.

The orientation angle Φ(t) that defines where vector a(t) is pointing (figure 2) is defined as the seismic phase at time coordinate t – hence the term “instantaneous phase.” Numerically, the phase angle is calculated using the middle equation listed on figure 2.

As time progresses, vector a(t) moves down the time axis, constantly rotating about the time axis as it maintains contact with the spiraling helical trace z(t).

Mathematically, frequency can be defined as the rate of change of phase. This fundamental definition allows instantaneous frequency ω(t) to be calculated from the time derivative of the phase function as shown by the bottom equation on figure 2.

The calculation of these three interpretation attributes – amplitude, phase and frequency – are illustrated on figures 3 and 4. Application of the three equations listed on figure 2 yields first the instantaneous amplitude for one seismic trace x1(t) (figure 3), and then instantaneous phase and frequency are shown on figure 4 for a different seismic trace x2(t).

Note that the instantaneous frequency function is occasionally negative – a concept that has great interpretation value, as has been discussed in a previous article (April 2008 Geophysical Corner).

For those of you who click on a menu choice to create a seismic attribute as you interpret seismic data, you now see what goes on behind the screen to create that attribute.

Comments (0)


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

Next month – Part two

In the two-part series on the role of reference surfaces in seismic interpretation: Defining reflection events and their polarities.

See Also: Bulletin Article

The continuity of clay smears evolving in sealed direct shear experiments of initially intact sandstone-mudrock sequences was quantified to large displacements up to more than ten times the thickness of the sealing layer. The sample blocks consisted of a preconsolidated clay-rich seal layer, which was embedded and synthetically cemented in quartz sand. The mineralogy and mechanical properties of the clay layer and the reservoir sandstones were varied systematically to mimic a range of natural clastic rock sequences. The fluid-flow response across the fault zone was monitored continuously during deformation using a new type of direct shear cell. The displacement at which seals break down is closely linked to the amount of phyllosilicates in the seal layer. Contrary to expectations, softer seal layers do not seal better than stiff seal layers for a given clay content. In the testing range of normal effective stresses between 4 to 24 MPa (580–3481 psi) covering maximum burial depth conditions of approximately 800 m (2625 ft) to approximately 4 km (2 mi) (assuming normal fault tectonics), a systematic trend is also observed, indicating better smear continuity by increasing the effective normal stress. Predominantly brittle processes such as slicing and wear, and not ductile drag or plastic flow, appear to be responsible for the generation of clay smears. The test results offer the prospect of incorporating critical shale smear factors (i.e., normalized displacement at which seal breakdown occurs) into probabilistic fault seal algorithms that consider important properties that can be measured or estimated, namely, clay content and fault-normal effective stress.
Desktop /Portals/0/PackFlashItemImages/WebReady/fault-transmissibility-in-clastic-argillaceous-sequences.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 3720 Bulletin Article

See Also: CD DVD

These proceedings contain 283 manuscripts and extended abstracts presented during the three days of plenary, ePaper and oral sessions and are stored conveniently on a 4GB flash drive. The content is a product of the joint venture between AAPG, SEG and SPE and the collaboration of key disciplines and technologies engaged in the development of North American resource plays.

Desktop /Portals/0/images/_site/AAPG-newlogo-vertical-morepadding.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 12766 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 4052 CD-DVD

See Also: Online e Symposium

The geochemistry of formation fluids (water and hydrocarbon gases) in the Uinta Basin, Utah, is evaluated at the regional scale based on fluid sampling and compilation of past records.

Desktop /Portals/0/PackFlashItemImages/WebReady/oc-es-green-river-shales-geochemical-basin-study.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 1436 Online e-Symposium

See Also: Short Course

The course is a practical and applied introduction to geochemical techniques routinely employed in shale-gas condensate and tight-oil reservoir assessment with an emphasis on tools and techniques. Participants should have a solid background in petroleum geology.
Desktop /Portals/0/PackFlashItemImages/WebReady/sc-introductory-geochemistry-for-condensate-rich-shales-and-tight-oil.jpg?width=50&h=50&mode=crop&anchor=middlecenter&quality=90amp;encoder=freeimage&progressive=true 10167 Short Course