Responding to the latest developments in rock physics research, this popular reference book has been thoroughly updated while retaining its comprehensive coverage of the fundamental theory, concepts, and laboratory results. It brings together the vast literature from the field to address the relationships between geophysical observations and the underlying physical properties of Earth materials - including water, hydrocarbons, gases, minerals, rocks, ice, magma and methane hydrates. This third edition includes expanded coverage of topics such as effective medium models, viscoelasticity, attenuation, anisotropy, electrical-elastic cross relations, and highlights applications in unconventional reservoirs. Appendices have been enhanced with new materials and properties, while worked examples (supplemented by online datasets and MATLAB® codes) enable readers to implement the workflows and models in practice. This significantly revised edition will continue to be the go-to reference for students and researchers interested in rock physics, near-surface geophysics, seismology, and professionals in the oil and gas industries.

The Rock Physics Handbook

From experimental studies in digital processing of seismic reflection data, geophysicists know that a seismic signal does vary in amplitude, shape, frequency and phase, versus propagation time To enhance the resolution of the seismic reflection method, we must investigate these variations in more detail. We present quantitative results of theoretical studies on propagation of plane waves for normal incidence, through perfectly elastic multilayered media. As wavelet shapes, we use zero-phase cosine wavelets modulated by a Gaussian envelope and the corresponding complex wavelets. A finite set of such wavelets, for an appropriate sampling of the frequency domain, may be taken as the basic wavelets for a Gabor expansion of any signal or trace in a two-dimensional (2-D) domain (time and frequency). We can then compute the wave propagation using complex functions and thereby obtain quantitative results including energy and phase of the propagating signals. These results appear as complex 2-D functions of time and frequency, i.e., as “instantaneous frequency spectra. ’ ’ Choosing a constant sampling rate on the logarithmic scale in the frequency domain leads to an appropriate sampling method for phase preservation of the complex signals or traces. For this purpose, we developed a Gabor expansion involving basic wavelets with a constant time duration/mean period ratio. For layered media, as found in sedimentary basins,

Wave propagation and sampling theory—Part I: Complex signal and scattering in multilayered media

The authors introduce a directionally oriented 2-D filter bank with the property that the individual channels may be critically sampled without loss of information. The passband regions of the component filters are wedge-shaped and thus provide directional information. It is shown that these filter bank outputs may be maximally decimated to achieve a minimum sample representation in a way that permits the original signal to be exactly reconstructed. The authors discuss the theory for directional decomposition and reconstruction. In addition, implementation issues are addressed where realizations based on both recursive and nonrecursive filters are considered. >

A filter bank for the directional decomposition of images: theory and design

Multifold ground coverage by seismic techniques such as the common reflection point method provides a multiplicity of wave travel path information which allows direct determination of root‐mean‐square velocities associated with such paths. Hyperbolic searches for semblance among appropriately gathered arrays of traces form the basis upon which velocities are estimated. Measured semblances are presented as a velocity spectral display. Interpretation of this information can give velocities with meaningful accuracy for primary as well as multiple events. In addition, the velocity data can help correctly label events. This paper outlines the fundamental principles for calculating velocity spectra displays. Examples are included which demonstrate the depth and detail of geological information which may be obtained from the interpretation of such displays.

Velocity spectra-digital computer derivation and applications of velocity functions

Preface 1. Introduction to rock physics 2. Rock physics interpretation of texture, lithology and compaction 3. Statistical rock physics: combining rock physics, information theory, and statistics to reduce uncertainty 4. Common techniques for quantitative seismic interpretation 5. Case studies: lithology and pore-fluid prediction from seismic data 6. Workflows and guide lines 7. Hands-on References Index.

Quantitative Seismic Interpretation: Applying Rock Physics Tools to Reduce Interpretation Risk

Offset-Dependent Reflectivity–Theory and Practice of AVO Analysis

Two-dimensional, fenced 2-D, and 3-D isosurface displays of some realistic 3-D seismic models built in the lower Miocene Powderhorn Field, Calhoun County, Texas, demonstrate that a seismic event does not necessarily follow an impedance boundary defined by a geological time surface. Instead, the position of a filtered impedance boundary relative to the geological time surface may vary with seismic frequency because of inadequate resolution of seismic data and to the en echelon or ramp arrangement of impedance anomalies of sandstone. Except for some relatively time-parallel seismic events, the correlation error of event picking is large enough to distort or even miss the majority of the target zone on stratal slices. In some cases, reflections from sandstone bodies in different depositional units interfere to form a single event and, in one instance, an event tying as many as six depositional units (interbedded sandy and shaly layers) over 50 m was observed. Frequency independence is a necessary condition for selecting time-parallel reference events. Instead of event picking, phantom mapping between such reference events is a better technique for picking stratal slices, making it possible to map detailed depositional facies within reservoir sequences routinely and reliably from 3-D seismic data.

Stratal slicing, Part I: Realistic 3-D seismic model

A new technique has been developed which makes it possible to process a seismic record‐section in such a way that all seismic events with dips in a given range are preserved with no alteration over a wide frequency band, while all seismic events with dips outside the specified range are uniformly and severely attenuated. By applying this process to a noisy record‐section, a record‐section may be obtained which has all events within a specified dip range perfectly preserved, and very high‐velocity noise essentially eliminated, a result which is impossible by simple wave‐number filtering or conventional array usage. In structurally complex areas where several steeply dipping events interfere, the technique may be applied to separate the events with different dips. In areas where a normal‐moveout contrast exists between primaries and multiples, the technique may be used for wide‐band multiple attenuation. By application of a “rotating Pie‐Slice” to micro‐spread noise data, seismic noise may be separated on t...

Wide-band velocity filtering; The Pie-Slice process

In offshore shooting the validity of previously recorded seismic data has been severely limited by multiple reflections within the water layer. The magnitude of this problem is dependent on the thickness and the nature of the boundaries of the water layer. The effect of the water layer is treated as a linear filtering mechanism, and it is suggested that most apparent water reverberation records probably contain some approximate subsurface structural information, even in their present form. The use of inverse filtering techniques for the removal or attenuation of the water reverberation effect is discussed. Examples show the application of the technique to conventional magnetically recorded offshore data. It has been found that the effectiveness of the method is strongly dependent on the instrumental parameters used in the recording of the original data.

Water reverberations—their nature and elimination

We discuss, in a two-part article, the benefits of 90°-phase wavelets in stratigraphic and lithologic interpretation of seismically thin beds. In Part 1, seismic models of Ricker wavelets with selected phases are constructed to assess interpretability of composite waveforms in increasingly complex geologic settings. Although superior for single surface and thick-layer interpretation, zero-phase seismic data are not optimal for interpreting beds thinner than a wavelength because their antisymmetric thin-bed responses tie to the reflectivity series rather than to impedance logs. Nonsymmetrical wavelets (e.g., minimum-phase wavelets) are generally not recommended for interpretation because their asymmetric composite waveforms have large side lobes. Integrated zero-phase traces are also less desirable because they lose high-frequency components in the integration process. However, the application of 90°-phase data consistently improves seismic interpretability. The unique symmetry of 90°-phase thin-bed response eliminates the dual polarity of thin-bed responses, resulting in better imagery of thin-bed geometry, impedance profiles, lithology, and stratigraphy. Less amplitude distortion and less stratigraphy-independent, thin-bed interference lead to more accurate acoustic impedance estimation from amplitude data and a better tie of seismic traces to lithology-indicative wireline logs. Field data applications are presented in part 2 of this article.

Milo M. Backus

Papers

Offset-Dependent Reflectivity–Theory and Practice of AVO Analysis

Stratal slicing, Part I: Realistic 3-D seismic model

Wide-band velocity filtering; The Pie-Slice process

Water reverberations—their nature and elimination

Interpretive advantages of 90°-phase wavelets: Part 1 — Modeling