Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information
Petroleum Engineering Handbook
Larry W. Lake, Editor-in-Chief
Volume V – Reservoir Engineering and Petrophysics
Edward D. Holstein, Editor
Copyright 2007, Society of Petroleum Engineers
Chapter 3C - Acoustic Logging
Petroleum applications of acoustic-wave-propagation theory and physics include both surface- and borehole-geophysical methods. These data-acquisition methods cover a broad range of scales from millimeters to hundreds of meters (Fig. 3C.1). Acoustic logging is a subset of borehole-geophysical acoustic techniques. This chapter provides an overview of borehole acoustic-logging theory, modern tool design, data processing methods, and data applications. Table 3C.1 lists other common surface- and borehole-geophysical methods. The chapter on the Fundamentals of Geophysics, in this volume, and the chapter on Reservoir Geophysics, in the Emerging and Peripheral Technologies volume of this Handbook, also discuss these methods.
Fig. 3C.1 – Diagram showing the maximum and minimum ranges vs. the resolution for various acoustic methods. Acoustic logging includes the sonic and ultrasonic ranges (courtesy of SPE).
A virtual explosion in the volume of acoustic research conducted over the past 20 years has resulted in significant advances in the fundamental understanding of downhole acoustic measurements. These advances, in turn, have greatly influenced practical logging technology by allowing logging-tool designs to be optimized for specific applications.
Borehole acoustic-logging measurements are used in a wide variety of geophysical, geological, and engineering applications and play an important role in evaluating reservoirs, reducing exploration and production risks, selecting well locations, designing completions, and increasing hydrocarbon recovery (Table 3C.2).
Modern logging tools include conventional borehole-compensated (BHC) monopole devices as well as the newer array devices—both monopole and multipole (monopole/dipole)—and logging-while-drilling (LWD) acoustic services. These logging tools provide acoustic measurements in all borehole mud types (but not in air- or foam-filled boreholes) in vertical, deviated, and horizontal wells, in both open and cased hole. They are combinable with other logging devices and are available in a variety of sizes to accommodate a range of borehole and casing diameters. Specialized tool designs are used for cement and casing evaluation and borehole imaging.
Historically, the primary and the most routine uses of acoustic logs in reservoir engineering have been porosity determination, identification of gas-bearing intervals, and cement evaluation. Continuing developments in tool hardware and in interpretation techniques have expanded the utility of these logs in formation evaluation and completion (fracture) design and evaluation. This chapter discusses the potential applications of these logs to allow the reader to evaluate the appropriate applications for a particular well or field. However, site-specific assessments are required to determine whether acoustic logs, with the proper planning, can provide the desired results in a cost-effective manner.
- 1 Acoustic Theory and Wave Propagation
- 2 Acoustic-Logging Tools
- 3 Evolution of Acoustic-Logging Tools
- 3.1 Velocity/Porosity Logging
- 3.2 Long-Spaced Acoustic Tools
- 3.3 Monopole Array Devices
- 3.4 Through-Casing Acoustic Measurements
- 3.5 This page contains the following errors:
- 3.6 Below is a rendering of the page up to the first error.
- 3.7 Logging While Drilling
- 3.8 Ultrasonic Reflection (Pulse/Echo) Acoustic Devices
- 3.9 Cement-Bond Logging
- 4 Conventional Applications
- 5 Advanced Data Analysis And Applications
- 6 Summary
- 7 Nomenclature
- 8 Subscripts
- 9 References
- 10 SI Metric Conversion Factors
Acoustic Theory and Wave Propagation
The principles of borehole acoustic logging (and surface seismic methods) are based on the theory of wave propagation in an elastic medium, as detailed in several sources. The oscillating motion generated by a sound source (transducer) in an elastic medium (rock formation) is called an elastic wave or acoustic wave (also called head or body waves). Wave theory predicts how an acoustic signal propagates through the borehole and formation. Snell’s law explains how the acoustic signal behaves at the velocity boundary separating the borehole and the formation, that is, how it is transmitted into the formation and back to the receivers. Elasticity is the property of matter that causes it to resist deformation in volume or shape. It is the elastic nature of rock formations that permits wave propagation. Acoustic waves have four measurable properties: velocity, amplitude, amplitude attenuation, and frequency. Acoustic logging tools are designed to measure one or more of these properties, with velocity (slowness) being the most common.
The waveform recorded at the logging tool’s receivers is a composite signal containing different energy modes, each with a different frequency, velocity, and amplitude. For borehole logging, the modes of primary interest are (in order of arrival) compressional, shear, and Stoneley (tube) waves (Fig. 3C.2). The waveform is recorded as acoustic amplitude as a function of time.
These waves are transmitted through the medium some distance from the origin of displacement. The particles of the medium do not travel with the wave, but only vibrate around their mean central position. Acoustic waves are classified according to the direction of particle displacement with respect to the direction of wave propagation as either longitudinal (i.e., particle displacement is parallel to the direction of propagation) or transverse (i.e., direction of particle displacement is perpendicular to the direction of propagation). In acoustic logging, the longitudinal wave is known as the compressional wave and the transverse wave is known as the shear wave. The presence of the borehole excites two additional acoustic energy modes, called guided waves: normal (pseudo-Rayleigh) and tube (Stoneley) waves.
Acoustic-wave velocity is controlled by a number of factors: lithology (mineralogy), cementation, clay content, texture, porosity, pore-fluid composition and saturation, overburden-and pore-fluid pressure (stress), and temperature. The rock’s mechanical properties, elastic dynamics, and density are a constant for a particular homogeneous and isotropic material. Acoustic-wave velocity can be related to rock elastic properties through three constants of proportionality, elastic moduli (e.g., Young’s, shear, and bulk), and Poisson’s ratio. This serves as the basis for mechanical-property evaluation by acoustic logs (see the additional discussion under Geomechanical Applications—Rock Mechanical Properties). In reality, most petroleum reservoirs contain varying pore sizes, pore fill (e.g., clays), fractures, etc. and consequently, are neither truly isotropic nor homogeneous. Furthermore, in fluid-saturated rocks, these acoustic properties also depend on the type and volume of fluids present.
Compressional WavesCompressional (P, primary, or pressure) waves are longitudinal waves that are transmitted through an elastic formation by compression or pressure. Particle motion is parallel to the direction of wave propagation (Fig. 3C.3). They can travel through solids, liquids, and gases and are the fastest wave type—they represent the acoustic first arrival. Of all acoustic wave types, they are least affected by faults, unconsolidated formations, and borehole fluids, and are, therefore, the most reliable. The wave is transmitted by both the rock matrix (i.e., the framework) and the fluid present in the pore throats. A compression, together with an adjacent rarefaction preceding or following it, constitutes a complete cycle. The distance between complete cycles is called the wave length and the number of cycles propagating through a point in the medium per unit time is the frequency. The velocity of elastic-wave propagation in an isotropic homogeneous medium can be derived from a combination of the theory of elasticity with Newton’s law of motion. Compressional-wave velocity (or travel time) is a function of the density and elasticity of the medium and is a constant for a given material.
Shear and Borehole Flexural Waves
Shear (S, secondary) waves are transverse waves that are transmitted by lateral displacement of particles in a rigid elastic formation. Particle motion is perpendicular to the direction of motion (Fig. 3C.3). Normally, shear waves are the second arrival in an acoustic wave train. In most reservoir rocks, shear waves generally have higher amplitudes than compressional waves but lower velocities, by as much 40 to 50%. There are two types of borehole shear waves—direct and indirect, also known as refracted or induced. Indirect shear waves are induced in a formation through a process known as mode conversion in which some of the compressional energy is transferred from the borehole fluid into the rock formation. Monopole transmitters generate these indirect shear waves while dipole transmitters generate direct shear waves by inducing a flexural (asymmetric mode) in the borehole. Shear-wave propagation requires a medium that has shear strength (rigidity). Consequently, shear waves can only travel in solids, not in liquids or gas. In liquids and gas, the shear head-wave generated within the formation is converted into a compressional wave and propagated back across the borehole fluid to the acoustic receivers as a later-arriving compressional wave.
Unconsolidated or poorly consolidated sandstones ("soft" or "slow" rocks) are less rigid and more compressible than well-consolidated ("hard" or "fast") rocks. When the formation shear-wave velocity is less than the acoustic velocity of the borehole fluid (Vs < Vf), a rock formation is called "slow." There is no refracted shear-wave from monopole devices in slow formations and low-frequency dipole transmission and reception is required to adequately detect low-frequency flexural arrivals for the shear-wave slowness determination. However, if a monopole-array tool is used in these conditions, a shear-wave slowness can be estimated from Stoneley-wave velocity dispersion. In very slow formations, where Vc < Vf, special processing may be required to extract the formation compressional signal.
Flexural-wave velocity varies with frequency—a phenomenon called dispersion. In contrast, the compressional and shear headwaves generated by monopole sources are generally not dispersive. At very low frequencies, the flexural wave travels at the formation shear velocity. This dispersion effect diminishes as the wavelength of the flexural-wave increases and is generally minimal when the wavelength is at least three times the borehole diameter. Fast formations exhibit a center frequency slightly greater than 3 kHz, while slow formations exhibit a center frequency of ≈ 1 kHz, or less. The received frequency spectrum of a dipole array is a function of transmitter frequency, rock properties, and borehole size. Modern dipole transmitters are broadband transmitters, i.e., they operate over a range of frequencies, to account for dispersion and to accommodate different formation types.
Stoneley (tube) waves are high-amplitude guided waves that are generated by a radial (symmetric) flexing of the borehole as the acoustic energy passes from the borehole fluid into the rock formation. They propagate at low frequencies along the fluid-rock interface at the borehole wall; hence, they are sensitive to the rock properties adjacent to the borehole wall. They are the slowest acoustic mode. They can be measured in both open and cased boreholes, but in cased holes Stoneley-wave features are primarily controlled by the casing rigidity. Similarly to shear waves, Stoneley waves are also dispersive; i.e., wave velocity varies with frequency—the amount of dispersion is related to formation rock properties. However, Stoneley waves are notable for several special properties: there is no cut-off frequency; dispersion is very mild; for all frequencies, Stoneley-wave velocity is less than fluid velocity; and group velocity nearly equals phase velocity over the frequency range.
All acoustic waves undergo attenuation, a reduction in signal amplitude away from the source. For logging devices this means radially away from the borehole wall. Signal attenuation results from the geometric spread of energy through reflection, refraction, and scattering, and through absorption by the medium through which the acoustic energy travels. Attenuation, usually expressed in dB/ft, is characteristic of different materials and increases with frequency of the acoustic wave. Generally, attenuation is large in slow formations and very small to negligible in fast formations. Because of these features, Stoneley waves are used to identify acoustic leakage away from the borehole that may be caused by formation permeability or the presence of fractures.
Acoustic-logging devices are comprised of transmitters (sources), receiver arrays, and accompanying electronics. They have been designed to measure one or more acoustic-wave properties. Acoustic sources (transmitters) generally consist of piezoelectric transducers that generate the acoustic signal by converting electrical signals into a sonic vibration that travels through the borehole and adjacent rock formations. Monopole (axisymmetric) transducers generate omnidirectional acoustic waves around the tool circumference, while dipole (nonaxisymmetric) transducers generate azimuthally oriented acoustic waves (Fig. 3C.3).
Modern receivers are piezoelectric crystals that transform the received (measured) acoustic signal back into electrical signals. Different logging tools use different piezoelectric materials and operate at different frequencies to measure different energy modes (wave types). The pressure variation produced by an acoustic wave displaces the piezoelectric material causing it to ring or oscillate. When the receiver oscillates, it develops a small voltage. This voltage is amplified and the raw data are processed downhole and sent to the surface-acquisition system in the form of a waveform signal (wireline) or acoustic slowness (wireline and LWD). Digital data are stored, either at the surface (wireline) or downhole (LWD), for wellsite log presentations and post-acquisition processing and playback.
The energy (amplitude) of acoustic waves is lost (attenuated or dispersed) primarily by travel through the borehole fluid and rock matrix. Additional attenuation may result from a number of other factors that include internal particle friction within the propagation medium, changes in acoustic impedance [the product of density (ρ) and acoustic velocity (V)] at interfaces (boundaries) between different mediums, borehole rugosity, and signal cancellation resulting from tool eccentering. In general, the largest signal possible occurs when the instrument is in the center of the hole; a dipole tool in an extremely large borehole may be exception to this rule. For high frequencies, the monopole signal is reduced by as much 50% of the centralized value by displacing the instrument only inches from the center of the hole. Consequently, whenever possible, centralizers should be used with acoustic-logging tools.
The acoustic wave travels through the formation much faster than it does through the borehole fluid and reaches the receiver by the longer formation route first. This is true as long as the critical spacing is less than actual spacing. Critical spacing is the transmitter-to-receiver spacing at which the fluid-signal and formation-signal arrive at the receiver at the same time. The critical-spacing value depends on the diameter of the logging sonde, the diameter of the borehole, the time of travel through the fluid, and the time of travel through the rock.
In soft formations, and under certain conditions, the critical spacing can exceed actual spacing and the fluid arrival, which exists for a noncentralized tool, can interfere with the acoustic signal. When this occurs, data-processing methods can exclude this fluid arrival, thereby reducing or possibly eliminating interference from the fluid arrival. Longer transmitter-receiver spacings are used to minimize this occurrence.
The transmitter-to-receiver spacing in modern monopole tools is set to enable separation of the compressional- and shear-energy packets to allow for accurate measurement of both in fast formations. When logging with an array-acoustic device, the receivers nearest the transmitter see stronger acoustic waves than the more distant receivers.
Monopole ExcitationThe transmitter emits acoustic energy uniformly around the tool. In fast formations, this energy excites three waves that travel down the borehole wall: compressional, shear, and Stoneley. The compressional wave travels away from the transmitter with a velocity, Vf, in the mud. When these waves reach the borehole face, they are reflected, refracted, and converted according to Snell’s law (Fig. 3C.4).
For angles of incidence less than the compressional-wave critical angle (θc, Fig. 3C.4), part of the energy is transmitted into the formation in the form of compressional wave (Wave A, Fig. 3C.4), another part is converted as a (refracted) shear wave (Wave B, Fig. 3C.4), and the remainder is reflected back into the mud as a compressional wave (Wave C, Fig. 3C.4). The transmitted waves travel at velocity Vp and Vs in the formation, close and parallel to the borehole wall, while continuously radiating energy back into the mud as converted compressional-waves, at the same compressional-wave critical angle at which it entered. It is this radiated energy that is detected by the receivers.
If the formation shear-wave velocity is slower than borehole-fluid compressional velocity (Vs < Vf), shear waves cannot be refracted along the borehole wall, and no shear wave is measured. Beyond the shear-wave critical angle, all the incident energy is reflected back into the mud to form the guided waves. The Stoneley wave travels at approximately the velocity of compressional waves in the borehole fluid. Compressional and Stoneley arrivals are always present. In the absence of a refracted shear-wave arrival (i.e., in formations in which Vp < Vf), the Stoneley wave can be used to estimate formation shear-wave velocity when a formation bulk-density measurement is available, using certain assumptions. However, because of the uncertainty associated with these assumptions, dipole shear-wave measurements are recommended in slow formations. Stoneley-wave amplitude decreases (attenuates) significantly at high frequencies and modern tools use low-frequency transmitters (< 1 to 12 kHz) to ensure acquisition of the Stoneley arrival in slow formations.
The dipole transmitter exerts a differential pressure on one side of the transmitter element that creates a flexural wave in the borehole, much like the wave produced when a vertical rope is shaken from side to side. The flexural wave is dispersive, but at low frequencies this wave travels down the borehole at the formation shear velocity. Receivers, sensitive only to differential pressures, are used to detect this flexural wave. Because the receivers are not sensitive to axially symmetric pressure fields, both the compressional head wave and the Stoneley waves are suppressed. This is desirable because it simplifies data processing. The desired output is the velocity of the formation shear wave. If the wavelength of the flexural wave is at least three times the diameter of the borehole, the flexural wave travels at very nearly the formation bulk-shear velocity. However, because this is a dispersive mode, if the wavelength is shorter (because of higher frequency), this flexural mode will travel slower than the shear velocity and dispersion corrections are needed.
All types of well logs, both open- and cased-hole, should be accompanied by complete documentation to ensure good-quality logs and sound interpretation. It is important to remember that any acoustic analysis represents an interpretation of measured acoustic waves. The wellsite engineer or geologist must ensure that all data pertaining to a particular log run, including borehole information (bit and casing sizes and depths), tool configuration, borehole fluid, formation parameters, and tool centralization are recorded in the well-log header for future reference. Cased-hole logs should also contain information on cement composition, casing weight and thickness, if the log was run with pressure, and the amount of that pressure.
Evolution of Acoustic-Logging Tools
The most commonly used acoustic-wave property acquired in borehole logging is the compressional-wave velocity. Modern velocity-logging tools measure the time, Δt, required for a compressional or shear wave to travel through a fixed distance of formation; it is recorded as a function of depth. This parameter, Δt, referred to as the interval transit time, transit time, travel time, or slowness, is the reciprocal of the velocity of the compressional waves, Δt = 1/ Vp. For the formations typically encountered in acoustic logging, travel times range from 40 to 250 μsec/ft, corresponding to velocities ranging from 25,000 to 4,000 ft/sec (Table 3C.3).
The resolution of any acoustic method is a function of the signal wavelength; the lower limit is one-quarter of the propagating wavelength. Seismic-reflection exploration methods typically operate in the frequency range of 10 to 50 Hz. In typical petroleum reservoir rocks, the resolution of these methods is approximately 30 to 160 ft (10 to 50 m), depending on depth (signal strength). In contrast, conventional BHC logging tools operate within the frequency range of 10 to 40 kHz, while newer array devices operate at even lower frequencies, 1 to 12 kHz. Logging-tool resolution is also a function of both the array aperture and the methodology employed for the array processing. Consequently, acoustic-logging tools typically have resolutions on the order of 1.0 to 4.0 ft (0.3 to 1.2 m) (see the section on Resolution Enhancement in this chapter).
Velocity/Porosity LoggingAcoustic well logging developed out of the need for downhole velocity (time-depth) measurements to improve the accuracy (calibrate) of surface seismic measurements. Surface seismic maps out sub-surface structures referenced to time and borehole acoustic tools provide a bridge to understand how time is related to depth. Downhole geophones were introduced in the 1930s to provide acoustic travel times to the surface, and continuous-velocity-logging tools were introduced in the 1950s (Table 3C.4). Soon after the introduction of the continuous-velocity log, it was recognized that these data also provided an excellent means for stratigraphic correlation, lithologic identification, and for evaluation of formation porosity.
These simple devices were soon replaced by two receiver designs that had the advantage of eliminating the need to correct for travel time in the drilling mud (for more on the historical development of acoustic logging, consult the literature). Modern borehole-compensated designs, introduced in the 1960s, use two transmitters and two receivers to compensate for variations in borehole diameter and tool position within the borehole (Fig. 3C.5). Two separate values of two receiver (R1 and R2) and two transmitters (T1 and T2) interval-transit times are provided, and an average of the two effectively compensates for any problems (Eq. 3C.1).
where t1 = travel time between T2 and R1; t2 = travel time between T2 and R2; t3 = travel time between T1 and R1; t4 = travel time between T1 and R2 ; and X12 = distance between R1 and R2. Fig. 3C.5 illustrates travel paths that show that the averages of AA′, BB′, and CC′ are essentially equal.
Conventional BHC monopole-acoustic logs, with their short transmitter-to-receiver spacing, have shallow depths of investigation and they largely measure mud filtrate that fills the pore space in the invaded (flushed) zone around the borehole. Long-spaced and array devices can acquire measurements beyond the filtrate and altered zone.
Long-Spaced Acoustic ToolsTypical BHC devices have a transmitter-to-receiver spacing (TR) of 3 to 5 ft. These work well in many circumstances; however, in cases which borehole enlargement prevents acquisition of reliable data, due to the increase of the critical spacing, or in which the drilling process damages or alters the shales surrounding the borehole, long-spacing tools (TR = 8 to 15 ft) may be necessary, or advised, to obtain accurate measurements. In contrast to conventional BHC devices, in which the transmitters and receivers are arranged symmetrically, long-spacing tools use an asymmetric arrangement with the receivers at varying distances from the monopole transmitter (Fig. 3C.6). Consequently, these devices have deeper depths of investigation that make them less susceptible to borehole conditions such as enlargement and shale alteration. These tools operate in both open and cased holes.
Because long-spaced and array tools use an asymmetric configuration of transmitters and receivers, borehole compensation is achieved through a process called depth-derived borehole compensation (DDBHC) (Fig. 3C.6). The processing is accomplished downhole by either the tool’s electronics or in the surface recording system and uses a depth-based delay to create synthetic transmitter arrays from multiple tool positions. The compensated travel-time measurement, Δt, is determined through the following procedure. At position three (Fig. 3C.6), the transmitter (T) is at the depth where the far receiver (R2) will be when the tool is moved to position one. The interval-transit time (A) between the transmitter (T) and near receiver (R1), which includes mud and formation signals, is recorded and delayed (memorized). At position two (Fig. 3C.6), the transmitter (T) is at the same depth the near receiver (R1) will be when the tool moves to position one. The interval transit time or waveform (B) between T and R2 is recorded and delayed. When the tool reaches position one, the two interval transit times (A and B) are equal to the interval transit time that would result if a second transmitter were below the receivers. The correct compensated value of Δt is obtained by combining the two delayed values of transit time (recorded at positions two and three) with transit times C and D, recorded at position one. Compensated transit time (Δt) is then correctly represented by
where x = the distance between R1 and R2 .
Analog recording of full acoustic waveforms, called "amplitude" logging, was developed in the 1960s. However, it was not until digital technology, instrumentation, and signal-processing methods were introduced in the late 1970s that the recording of full-waveform data became routine. While these techniques enabled the extraction of shear-wave data from conventional BHC acoustic data, determination of shear arrivals using conventionally spaced BHC tools (TR = 3 to 5 ft) suffers from interference between late-arriving compressional waves and shear-wave arrivals. The use of long-spaced tools reduced this problem by allowing greater temporal separation between the different wave packets and provided accurate estimates of shear–wave slowness.
Monopole Array Devices
Modern array tools are a natural outgrowth of the long-spaced tool design. Additional receivers (4 to 13) were added to provide the statistical redundancy needed to enhance extraction of wave arrival times; some designs also include multiple transmitters. The monopole transmitters in these tools use lower frequencies (e.g., 1 to 12 kHz vs. 20 to 40 kHz) and have broader frequency ranges than earlier tools to permit acquisition of high-quality acoustic waveforms. These devices are typically comprised of several sections or subs that house the tool components: electronics, receiver array, acoustic isolator, and transmitters. The acoustic isolator, placed between the transmitter and receiver sections, prevents or minimizes and delays direct sound transmission between transmitters and receivers. The electronics section provides timing and control for the transmitter and receiver sections, digitizes the received acoustic waves, analyzes the acquired waveforms, and transmits the data to the surface-data-acquisition system, all in real time.
Array tools can record full waveforms: compressional, shear, and Stoneley arrivals. These tools operate in either a single or a variety of combination-acquisition modes that include full waveform, compression Δt, and cement-bond logging. The number of modes that can be activated during a single log run is a function of the logging speed. Acquisition of full waveforms permits the use of waveform-correlation techniques for waveform amplitude, coherent slowness from the coherent-wave moveout, arrival-time processing, and most importantly, allows for detailed post-acquisition processing that improves the interpreted results. Because these techniques are insensitive to cycle skipping, they are particularly effective in gas-saturated, rugose, and washed-out boreholes.
Through-Casing Acoustic Measurements
The extended transmitter-to-receiver offset (6 to 19 ft) provided by array instruments allows them to quantify formation-compressional and shear-wave energy through casing. Successful cased-hole operation requires a good cement bond to provide the necessary acoustic coupling to the formation and to minimize or eliminate the casing arrival. New processing techniques may allow valid acoustic evaluation even in cases of poor bonding.<html><parsererror style="display: block; white-space: pre; border: 2px solid #c77; padding: 0 1em 0 1em; margin: 1em; background-color: #fdd; color: black">
This page contains the following errors:
Below is a rendering of the page up to the first error.
Tools using dipole transmitters were conceived as early as the 1960s, but were not actually developed until the 1980s. In contrast to monopole logging tools, dipole acoustic devices can excite a low-frequency flexural wave in the borehole at shear velocity. Low-frequency (< 1 kHz) dipole sources allow for shear-velocity determination that is much closer to seismic shear waves and permits acquisition of direct-shear velocities in slow and fast formations. However, increased noise (i.e., a lower signal-to-noise ratio) is one limitation of low-frequency operation. Noise has been reduced through improved acquisition electronics, the use of semi-rigid tool designs, and by choosing the operational mode of the dipole source. A semi-rigid tool body not only reduces the influence of the tool body on the measurement but also permits operation in deviated wells.
At high frequencies, or when the borehole diameter is large, flexural-mode propagation is slower and a dispersion correction is needed to obtain the shear velocity from the measured flexural velocity. This dispersion correction is a function of mud compressional velocity, formation compressional and shear velocities, the ratio of formation and mud densities, and the product of borehole diameter and processing frequency. Few, if any corrections are required if the flexural wavelength (velocity/frequency) is at least three times the borehole diameter, which is why low frequencies (< 1 kHz) are used. Where a correction is necessary, it is typically only a few percent, but can be higher under certain conditions.
The latest commercial tool designs are multipole array devices that operate in open and cased hole. These tools typically integrate multiple transmitters (monopole and dipole) and one or more arrays of monopole and dipole receivers. Multiple monopole transmitters, or a single "programmable" transmitter, provide the preferred frequency-range for optimal acquisition of conventional BHC and long-spaced compressional, shear, low-frequency Stoneley modes along with cement-bond logs. The dipole transmitter, more recently a variable-frequency (wide bandwidth) transmitter, provides the crossed-dipole mode. The receivers are positioned axially along the length of the tool, for short- and long-spaced measurements, and may be interlaced or independent. At each axial position a group of receivers (e.g., 4 or 8) is positioned azimuthally around the circumference of the sonde. The receivers are oriented with the transmitter to allow for alignment of directional (multimode) source excitation and data acquisition. This allows radial imaging of acoustic parameters and measured properties. The dipole receivers may be aligned inline with the dipole transmitters or orthogonal (crossed) to them for crossed-dipole analysis.
Multipole tools provide enhanced Stoneley-wave and crossed-dipole shear-wave data for analysis of formation permeability and anisotropy (e.g., stress and fractures) (see later discussion).
Logging While DrillingAcoustic slowness measurements are a relatively recent addition to the suite of LWD measurements. Similar to wireline devices, LWD tools consist of transmitter, isolator, and receiver-array sections contained in either a single or separate drill-collars and use monopole-type (axisymmetric) transmitters. However, unlike wireline devices, which are small relative to the borehole size, the rigid collar, structural design, and large diameter of acoustic LWD drill collars may actually interfere with the physics of acoustic-energy propagation, making it more difficult to decouple the transmitter-to-receiver signal traveling along the tool body and complicating generation of borehole guided modes, such as Stoneley and flexural.
Downhole processing provides acoustic slowness data in real time and also allows for storage of raw waveform data in downhole memory that can be downloaded during bit trips. The receiver array provides data redundancy and is coupled with a much narrower recording bandwidth to enhance the signal-to-noise ratio in the presence of the drilling-related noise. This results in more accurate measurements of acoustic slowness. As with all types of LWD tools, the primary advantage of LWD acoustic measurements is the acquisition of data before significant fluid invasion or alteration of the formation can occur, and providing the data in sufficient time to influence drilling decisions for improved safety, well placement, and well productivity. Table 3C.6 summarizes the advantages and disadvantages of LWD acoustic measurements. LWD acoustic data recorded in casing during bit trips offers the potential for cased-hole evaluation, e.g. formation velocity and cement-evaluation, during bit trips.
Acoustic LWD can provide a limited real-time seismic-while-drilling capability in combination with a surface source. There is also great interest among drillers and geophysicists to develop a viable drilling-based system in which the drillbit acts as the acoustic source and the sensors (ruggized geophones) are located in the borehole, as part of an LWD bottomhole assembly. Both types of systems provide checkshot surveys for time-depth correlation with pre-drill surveys and look-ahead seismic capability.
Ultrasonic Reflection (Pulse/Echo) Acoustic DevicesReflection (pulse/echo) acoustic devices were introduced in 1967 with the borehole televiewer (BHTV). In contrast to conventional acoustic-logging devices, which record the transmission of acoustic waves through the formation, pulse/echo devices record the travel time, amplitude, and azimuth of ultrasonic acoustic pulses (echoes) that are reflected off the formation wall in openhole or off the casing or cement in cased hole (Fig. 3C.8). The difference between the acoustic impedance of the borehole fluid and formation determines the magnitude (amplitude) of the transmitted ultrasound pulse that is reflected off the formation and back to the transducer. A portion of the transmitted signal is not reflected and continues to travel into the formation (openhole) or casing (cased hole) (Fig. 3C.8).
The peak amplitude of the reflected signal is used to generate a 360° image of the borehole wall and the travel-time measurements are used as a caliper to provide a measurement of borehole geometry or casing corrosion. The primary factors that contribute to the measured pulse/echo amplitude are borehole-fluid ultrasonic attenuation, borehole-fluid/formation reflection coefficient, physical features of the formation, and transducer-beam angle of incidence on the formation. Loss in signal amplitude (image quality) results from conditions that either scatter, absorb, or spread the acoustic energy, such as tool eccentering, irregularities in the borehole shape and surface, high-density and some oil-based drilling muds, and contrasts in acoustic impedance between the borehole fluid and borehole wall or casing (Table 3C.7). A separate transducer, commonly in the "mud-sub," continually measures the borehole-fluid velocity for the caliper measurement.
Cement-Bond LoggingProper cement placement between the well casing and the formation is essential to support the casing (shear bond), to prevent fluid from leaking to the surface, and for isolating producing zones from water-bearing zones (hydraulic bond). Acoustic logs provide the primary means for evaluating the mechanical integrity and quality of the cement bond. Acoustic logs do not measure cement quality directly, rather, this value is inferred from the degree of acoustic coupling of the cement to the casing and to the formation. Properly run and interpreted, cement-bond logs (CBL) provide highly reliable estimates of well integrity and zone isolation. Just as filtrate invasion and formation alteration may produce changes in formation acoustic properties, and thus variation in acoustic logs over time, so too, cement-bond logs may vary over time as the cement cures and its properties change.
Modern acoustic cement-evaluation (bond) devices are comprised of monopole (axisymmetric) transmitters (one or more) and receivers (two or more). They operate on the principle that acoustic amplitude is rapidly attenuated in good cement bond but not in partial bond or free pipe. These cased-hole wireline tools measure compressional-wave travel time (transit time), amplitude (first pipe arrival), and attenuation per unit distance. Conventional CBL tools provide omnidirectional measurements, while the newer radial cement-evaluation tools provide azimuthally sensitive measurements for channel evaluation.
When the acoustic wave generated by the transmitter reaches the casing, part is refracted down the casing (amplitude and travel-time measurement), part travels through the mud (fluid arrival), and other parts are refracted into the annulus and the formation and received back (formation arrival). Amplitude, measured directly or as an attenuation ratio, is the primary bond measurement and is used to provide quantitative estimations of cement compressive strength, bond index, and qualitative interpretation of the cement-to-formation interface. Tool response depends on the acoustic impedance of the cement, which, in turn is function of density and velocity. On the basis of empirical data, the log can be calibrated directly in terms of cement compressive strength. However, in foamed cements or when exotic additives are used, these calibrations can be inaccurate. In these situations, users are advised to consult with the logging service company regarding the appropriate calibrations.
A typical cement-log presentation includes a correlation curve (gamma ray), travel time (μsec), amplitude (mV), and attenuation (dB/ft) curves, and a full-waveform display (μsec). Presentation of the full acoustic waveform assists in resolving bond ambiguities arising from use of an amplitude measurement alone and provides qualitative information about the cement-to-formation bond. Waveform displays may be in variable density (VDL) or intensity (also called microseismograms) formats, oscilloscope waves (also known as x-y or "signature"), or both (Fig. 3C.11).
TR spacing typically ranges from 3 to 5 ft. The shorter spacing (e.g., 3 ft) provides optimum signal level and resolution at high attenuation rates and is normally used for amplitude and travel-time (TT) measurements. A longer spacing, commonly ≥ 5 ft, is used for the full-waveform recording because longer TR spacing provides greater separation of the casing and formation-signal arrival times. This separation allows for easier analysis of the formation-signal strength and is used to monitor cement-to-formation bonding.
These tools typically operate at higher frequencies than conventional openhole tools—between 20 and 30 kHz. As with openhole tools, cement-bond tools require centralization to ensure accurate measurements. Centering in the cased hole is more critical because the higher-operating frequencies (i.e., shorter wavelengths) and the tool measurement are based on signal amplitude. Tool eccentering reduces signal amplitude and travel time (Fig. 3C.14).
A casing cement job can result in one or more of the following situations: free pipe, good bond, bond to casing only, and partial bond. In the first scenario, free pipe, there is no cement bond between the casing and cement. Consequently, there is no acoustic coupling with the formation and most of the transmitted acoustic energy is confined to the casing and the borehole fluid. As a result, a free-pipe acoustic signal is long-lived, high-amplitude, and of uniform frequency.
In the second scenario, good bond, cement is properly bonded to casing and to the formation. This provides good acoustic coupling and most of the acoustic energy is transmitted to the formation, resulting in little (weak) to no casing signals and little amplitude until the arrival of the strong formation signal.
The third scenario, bond to casing only, is a common condition in which cement is bonded to the casing but not to the formation. This can occur because the mudcake dries and shrinks away from cement, or because the cement did not bond with mudcake in poorly consolidated formations. In this situation, energy traveling through the casing is attenuated drastically because of the highly attenuating cement sheath. At the same time, the annulus outside the cement sheath provides poor acoustic coupling. The result is that little energy is transferred to the annular fluid and virtually none is transferred to the formation. This condition is indicated by the lack of later-arriving formation energy. A similar response can be caused by the presence of formation gas in shallow, high-porosity zones.
In the last scenario, partial bond, a space exists within an otherwise well-bonded casing. This may occur with the presence of a microannulus or channels within the cement. The resulting waveform is comprised of a casing signal and a formation signal; the casing signal arrives first, followed by the formation signal.
When channeling occurs, it is generally localized and nonuniform; that is, it occurs over relatively short intervals and can frequently be identified by variations in the amplitude response. Channeling is significant because it prevents a hydraulic seal. In contrast, a microannulus (a small gap between the casing and cement sheath) may extend over long sections of casing but may not prevent a hydraulic seal. Microannulus may result from thermal expansion or contraction of the pipe during cementing or to the presence of contaminants, such as grease or mill varnish, on the casing’s exterior surface. A common practice is to run cement-bond logs with the casing under pressure to expand the casing against the cement, thereby decreasing any microannulus that might exist. If the initial log run was not under pressure and the log indicates poor bond, the presence of a microannulus can be evaluated by running a second bond log under pressure to see if there is a difference. Pressuring the casing improves the acoustic coupling to the formation and the casing signal will decrease and the formation signal will become more obvious (Fig. 3C.15). However, if only channeling exists, pressuring the casing will not significantly change the log. When conducting a cement evaluation, information on the type of cement used is essential. For example, foam cements, which intentionally create void spaces in the cured cement, can be misinterpreted as partial bond if normal cement is assumed. Fig. 3C.16 summarizes this discussion, and Table 3C.8 lists additional factors that may affect interpretation of bond quality from the amplitude response.
The bond index (BI) is a qualitative measure of cement bond based on signal amplitude. This dimensionless quantity is the ratio of measured attenuation to maximum attenuation:
A BI value of 1.0 represents a perfect cement bond. A value of less than 1.0 indicates an incomplete bond. This technique requires attenuation measurements in zones with 100% bond and in free pipe.
Radial-Cement Evaluation. Radial-cement-evaluation devices were developed to overcome some limitations of conventional cement-bond tools and to permit more accurate evaluation of cement distribution behind casing by providing the precise location of partial bond and channeling. These tools use one or more azimuthally sensitive transducers to evaluate cement quality around the circumference of the casing. Data from these tools are presented as individual log curves or as azimuthal images ("maps") of cement quality generated by interpolating between the individual azimuthal measurements (Fig. 3C.17). In addition, each tool design also provides a conventional 5-ft VDL waveform measurement to provide information about the cement-to-formation bond.
The ultrasonic tools compute the acoustic impedance of the material beyond the casing. To do this, repeated acoustic pulses are directed at the casing to make it resonate in its thickness mode and the energy level (attenuation) of the decaying reflected wave is measured. Good cement bond to casing produces a rapid damping (higher impedance) of this resonance; poor cement bond results in longer resonance decay (lower impedance). Measurements from these devices are influenced by the same factors as openhole televiewer devices.
The pad device makes multiple short-spaced, compensated, azimuthal-attenuation measurements. Because the pads are in direct contact with the casing—in contrast to ultrasonic measurements—measurements are unaffected by gas in the borehole, fast formations, heavy-mud conditions, or minor tool eccentricity. The attenuation in each segment is measured in two directions using a pair of acoustic receivers and two transmitters. The two measurements are combined to form a result that compensates for surface roughness and the effects of minor residual cement on the inside of the casing. Transmitting elements and the firing sequence are controlled to direct (steer) and enhance the acoustic-energy output of both the pad transmitters and the VDL transmitter. This has the effect of improving the signal strength of both the casing and cement-to-formation arrivals, respectively. This technique improves VDL interpretation, particularly in soft formations in which the standard VDL may wash out.
The use of new high-performance low-density, foam, and complex cements is increasing. However, the presence of gas in cement slurries, as an inert component or as contamination, may seriously affect ultrasonic-tool interpretation. New interpretation methods integrate ultrasonic and attenuation measurements from conventional tools to provide improved cement evaluation in these conditions. The latest ultrasonic tool has a conventional pulse-echo transducer plus a flexural transmitter and two flexural receivers that provide greater depth of investigation. Interpretation techniques combining these different measurements provide improved evaluation in lightweight cements, especially in the annulus, beyond the casing-cement bond. Table 3C.10 summarizes the capabilities and guidelines for running the different types of cement-bond evaluation tools.
Fig. 3C.18 – Casing evaluation log display. Holes in the casing are visible in the series of ultrasonic images that are based on amplitude (left) and corrected travel time (right). The center 3D images show the pipe in 90° quadrants. The image shading is generated from the amplitude data (courtesy of SPE).
The acoustic caliper generated from the pulse/echo travel time provides the casing inside diameter (an average of all transducers or a single circumferential scan). An estimate of casing ovality is obtained using only the maximum and minimum measurements. Then, if the nominal value of the outside casing diameter is assumed, changes in thickness can be calculated and internal defects identified. Frequency analysis determines the casing resonant frequency from the acoustic waveform; casing thickness is inversely related to the resonant frequency. By combining travel time and resonant-frequency measurements and using data from all available transducers (or a single scan), presentations showing casing cross sections are used to highlight casing damage such as thinning, corrosion metal loss, and collapse (Fig. 3C.19).
Interpretation of acoustic-log data begins with the slowness values obtained from processing the recorded waveforms. Slowness, or porosity derived from slowness, may be corrected for additional factors (discussed below) before use in applications. Today, log analysis and interpretation is routinely performed by computers during data acquisition (in real time) or in post-acquisition processing in offices and computing centers. Modern interpretation software is built on algorithms derived from the empirical relationships originally used to create a variety of graphical solutions contained in crossplots and nomograms. Because porosity tools vary in their sensitivity to lithology, porosity, and fluid properties, the combination of different measurements allows more-accurate determination of porosity, petrophysical, and geological characteristics. The logging service companies issue chartbooks with tool-specific charts to facilitate rapid manual determination of porosity, lithology, shaly-sand analysis, saturation, mechanical properties, and cement-bond evaluation. These charts, which are accompanied by instructions for their use, form the basis for quick-look interpretations and for quality control of computer interpretations.
The earliest applications of acoustic-logging measurements made use of compressional slowness, used alone or in combination with resistivity logs and other porosity logs. This group of applications, which includes velocity, porosity, gas identification, lithology, fluid saturation, and rock strength, still comprises the basic set of applications for acoustic-log data. In the 1960s, laboratory studies indicated that shear-wave data could also prove useful in formation evaluation. The development and introduction of reliable shear-wave acquisition and full-waveform logging in the late 1970s and 1980s resulted in improved accuracy in these conventional applications. Further advances throughout the 1980s and 1990s, in both basic and applied research as well as in tool technology, led to new and advanced applications, e.g., permeability estimation, anisotropy determination, and direct hydrocarbon indicators, that use shear- and Stoneley-wave amplitude and attenuation (Table 3C.11).
In terms of velocity, v:
where ϕ = fractional porosity of the rock, v = velocity of the formation (ft/sec), vf = velocity of interstital fluids (ft/sec) and, vma = velocity of the rock matrix (ft/sec).
In terms of transit time (Δt):
where Δt = acoustic transit time (μsec/ft), Δtf = acoustic transit time of interstitial fluids (μsec/ft), and Δtma = acoustic transit time of the rock matrix (μsec/ft). (See Table 3C.3 for typical values of Δtma and Δtf.)
In terms of travel time:
The velocity of most borehole and reservoir fluids (except gas) does not vary greatly and a fluid velocity (Δtf) of 189 μsec/ft (5,300 ft/sec) is generally assumed for fresh drilling fluids; a slightly lower value, 185 μsec/ft, is used for salt muds. Fluid type becomes more of a concern when oil-based mud (OBM) is used if the formation of interest is not invaded or if invasion is very shallow. The lithology must be known or estimated in order to select the appropriate matrix velocity.
The Wyllie equation represents consolidated and compacted formations. In poorly consolidated or unconsolidated rocks, a correction factor is necessary (Eq. 3C.8). Also, the presence of shale or clay within the sand matrix will increase Δt by an amount proportional to the bulk-volume fraction of the clay. An empirical equation is used for calculating porosity in sandstones in which adjacent shale values (Δtsh) exceed 100 μsec/ft (Eq. 3C.9):
where the compaction correction factor Cp is
where Δtsh = specific acoustic transit time in adjacent shales (μsec/ft), and 100 = acoustic transit time in compacted shales (μsec/ft). The shale compaction coefficient (C) generally ranges from 1.0 to 1.3, depending on the regional geology.
The highest velocities observed in sandstones approach 20,000 ft/sec (50 μsec/ft), but most sandstones have a lower matrix velocity. Velocities in adjacent shales are used to adjust the matrix velocity for sands with velocities lower than 18,000 ft/sec. Table 3C.12 provides guidelines for selecting the appropriate value of Δtma. If the lithology of carbonate rocks can be reasonably estimated and if the porosity distribution is fairly uniform, the Wyllie time-average formula can provide reliable determination of porosity for this group.
where α = (Δtma/2Δtf) − 1.
Graphical solutions for both algorithms for a sandstone matrix are compared in Fig. 3C.20. One caveat regarding the use of empirically derived porosity transforms: they do not account for all the factors influencing acoustic velocity. Consequently, these relationships may not be valid for all reservoirs.
These velocity/porosity methods are for clean (shale free), water-filled formations. The calculated apparent porosity must still be corrected for the volume of pore-filling material (shale). If the formation contains shale or dispersed-clay particles, or is hydrocarbon bearing and invaded to only a very shallow depth, corrections to the basic log data are necessary before reasonable porosity values can be calculated.
Because shale transit times range from 62 to 167 μsec/ft, failure to correct for the presence of shale may result in overly optimistic porosity calculations. The acoustic measurement is also influenced by the way the shale is distributed within the sandstone reservoirs. The fraction of shale, or shale volume, can be estimated using a combination of log measurements that are influenced by shale, such as neutron porosity, density, gamma ray, or spontaneous potential (SP). Chartbook nomograms developed for porosity determination include graphical solutions for both undercompaction and shale volume.
In some producing regions, producibility indexes based on the volume of shale in producing sandstone reservoirs have been developed. The fraction of total porosity occupied by dispersed clay (q factor) is empirically related to effective and total porosity and production characteristics. Local experience is used to create permeability cutoffs using the q factor (Fig. 3C.21).
An additional empirical velocity-porosity predictive model has recently been proposed and is still in the experimental phase.
Carbonate and complex lithology reservoirs are generally comprised of varying proportions of limestone, dolomite, chert, quartzite, and occasionally, evaporites. The primary influences on porosity in these rocks are lithology and pore type. Generally, any shale present is in dispersed form and in small amounts that do not significantly impact porosity calculations. Acoustic porosity is a measure of the primary or intergranular (matrix) porosity. In dual-porosity reservoirs, the secondary porosity (e.g., isolated pores, vugs, and fractures) may significantly influence the rock-pore distribution, but may be overlooked by acoustic-log measurements. This topic is the subject of ongoing research. In contrast, nuclear-porosity devices, such as density and neutrons, measure total porosity. The difference between the nuclear-porosity and acoustic-porosity measurements is an approximation of the secondary porosity.
Lithology Identification. Acoustic velocity is primarily a function of the rock matrix and can be used to identify different lithologies and for stratigraphic correlations. A variety of crossplot techniques, using acoustic measurements alone, or in combination with other porosity logs (neutron and density), have been devised to assist in lithologic identification (Fig. 3C.22). In particular, the M-N and mineral-identification-plot (MID) techniques use all three porosity logs in different combinations. Before lithology determination, the individual log measurements must be corrected for influences of gas effect, secondary porosity, bad hole conditions, and shaliness. In general, the MID plot is more sensitivity to lithology, gas, and secondary porosity and provides superior results to M-N plots. Crossplots, using a variety of log measurements or combinations of measurements can be used to resolve specific lithologic problems related to local or regional geology.
Fig. 3C.23 - Δtp vs. Δts crossplot demonstrates sensitivity to variations in lithology (courtesy of SPE).
Fig. 3C.24 – Shear slowness as a lithology indicator. The cased-hole apparent compensated-neutron porosity (Track 2) appears to increase with depth, while the dipole shear-velocity data decrease with depth. This indicates a lithology change rather than a porosity change, which is confirmed in the lithology track (courtesy of Baker Atlas).
Fig. 3C.25 - Vp/Vs vs. Δtp crossplot. In addition to the low-saturation sands identified in (a), the high-resolution result (b) also resolves turbidite and dry sand intervals (courtesy of SPWLA).
Compressional velocities are affected (slowed) by the compressive fluids in the pore space, while shear velocity is affected only by the rock matrix. Consequently, the presence of gas is especially noticeable in compressional-wave slowness. The combination of compressional and shear slowness, either as a ratio or as a log overlay, provides a quick-look gas indicator (Fig. 3C.27).
Fig. 3C.28 – Crossplot of shear slowness vs. Vp/Vs serves as a quick-look hydrocarbon indicator. Compressional-slowness travel times increase where the rock pore space is occupied by more compressible fluid (modified from original, courtesy of Baker Atlas).
Oil. Recent work suggests that the Vp/Vs ratio may also serve as an indicator of bypassed oil in cased wells. Research on the acoustic properties of heavy oils indicates that under the proper conditions of temperature and viscosity, these oils may behave as solids and generate shear waves that may be detectable at logging-tool frequencies.
The higher operating frequency of acoustic-logging tools and the smaller TR distances allows for higher-quality velocity data and finer vertical resolution than surface reflection techniques. Acoustic-velocity logs were originally developed for calibrating surface seismic velocities and reflectors. Acoustic-log interval travel time or transit time, Δt, can be summed, i.e., integrated, over the entire logged interval to provide the equivalent of seismic one-way time which is compared to borehole seismic surveys and reflection seismic two-way time.
Acoustic-log data are commonly calibrated using checkshot (velocity) or vertical seismic profile (VSP) surveys prior to use in geophysical applications. Data from these surveys, which use downhole receivers and surface acoustic sources, are used to adjust the log data for drift and borehole conditions and result in improved time-depth correlation. Acoustic-log data are combined with density-log data, to generate an impedance log that in turn is used to produce a synthetic seismogram. Synthetic seismograms are artificial seismic records that tie seismic time to log depth and are also used to match well-log quantities to seismic attributes for distinguishing primary seismic events (geologic structure and stratigraphy). It is possible, however, that a synthetic seismogram may not provide a very good match to the seismic field data. Disagreements commonly result from the differences scale and acquisition physics used in seismic and well-log measurement; for example, operational frequency (wavelength), borehole condition, and angle of measurement (particularly in the presence of anisotropy). Acoustic-log data provide a fundamental and essential element of modern seismic reservoir characterization. The chapter on Fundamentals of Geophysics in this volume of the Handbook contains more information on the determination and use of these types of analyses.
Drilling and Reservoir Engineering ApplicationsPore Pressure and Overpressure Detection and Evaluation. Abnormal pressure is defined as any departure from normal hydrostatic pressure at a given depth. Abnormal subsurface pressures, either overpressure (geopressure) or underpressure, are encountered in hydrocarbon basins throughout the world in all lithologies, all geologic ages, and at all depths. Early and reliable detection of geopressure is vital to avoid or mitigate potential drilling and safety hazards, e.g., shallow water flow, blowouts, and shale instability. During drilling, advanced warning of approaching geopressuring enables the mud weight to be adjusted to avoid well and reservoir damage and to determine casing points. This is a particular concern in deepwater wells in which the pressure difference; i.e., the operating window, between the hydrostatic gradient and the fracture gradient can be very narrow.
Geopressuring in hydrocarbon reservoirs may result from a variety of geologic and tectonic processes. Borehole-acoustic detection methods using compressional and shear slowness can identify abnormally pressured zones before they are drilled and can quantify pressure gradients. These methods, discussed below, are used in conventional borehole logging (wireline and LWD), new seismic-while-drilling techniques, and more recently, surface seismic data.
Undercompaction is the primary mechanism for creating overpressure, particularly in deltaic basins in which high rates of deposition commonly prevent the escape of pore water trapped in shales. Undercompacted shales have higher acoustic transit times (i.e., higher apparent porosity) than normally pressured shales at the same depth. With the onset of overpressuring, a semi-logarithmic plot of acoustic slowness with depth will diverge from a normal (hydrostatic) straight-line trend of decreasing slowness (increasing velocity) with depth (Fig. 3C.29).
Fig. 3C.29 – Semi-log plot illustrating acoustic detection of geopressure by use of shale slowness (courtesy of SPE).
Fig. 3C.30 – Relationship between fluid-pressure gradient (FPG) and the acoustic slowness difference for U.S. Gulf Coast (courtesy of SPE).
Recent investigations into the effects of pressure on shale porosity suggest that the relationship is more complex than previously thought. While additional study is necessary, the results to-date suggest that it may be necessary to reconsider or revise the well-log methods currently used in pore pressure and exhumation analysis (see section on Geological Applications below).
Fracture Identification. Locating fractures, recognizing fracture morphology, and identifying fluid-flow properties in the fracture system are important criteria in characterizing reservoirs that produce predominantly from fracture systems. However, fracture identification and evaluation using conventional resistivity and compressional-wave acoustic logs is difficult, in part because fracture recognition is very dependent on the dip angle of fractures with respect to the borehole.
Fractures are physical discontinuities that generate acoustic reflection, refraction, and mode conversion—all of which contribute to a loss of transmitted acoustic energy. In particular, compressional- and shear-wave amplitude and attenuation and Stoneley-wave attenuation are significantly affected by the presence of fractures. Compressional waves are primarily affected by oblique fractures—those with dip angles between 15° and 85°—while shear waves are primarily affected by horizontal or near-horizontal fractures. On conventional-velocity logs, fracture-induced attenuation may be evidenced as cycle skipping, variations in the Vp/Vs ratio and on VDL presentations, or as chevron (crisscross) patterns caused by mode-conversion interference. Borehole-televiewer-type imaging devices provide a higher degree of success in identifying fractures and determining whether or not they are open (producible) or closed. The development of reliable full-waveform shear- and borehole-imaging devices enabled enhanced fracture identification and evaluation. Aguilera summarizes the use of conventional acoustic-log methods for fracture identification.
Recently developed anisotropy-analysis methods use crossed-dipole shear, Stoneley-wave, and acoustic-imaging data—individually or in combination—to provide reliable identification and evaluation of in-situ and induced fractures (see the Crossed-Dipole Anisotropy Analysis section of this chapter).
Geological ApplicationsEstimates of Erosion and Uplift (Exhumation). The amount of erosion that has occurred in a region that has been uplifted can be estimated from the degree of shale compaction measured by acoustic travel time. This technique assumes that shale compaction is irreversible and that the shale retains the degree of compaction it gained at its maximum burial depth. Uplift and erosion will result in lower porosities than expected for the current burial depth (i.e., a shale will appear to be overcompacted).
Determination of Organic Richness and Source-Rock Potential. Acoustic slowness, used alone or in conjunction with formation resistivity, can provide qualitative indications and quantitative determination of source-rock potential (when calibrated to laboratory data). The identification of potential petroleum-source rocks and characterizing the thermal maturity of these rocks is important for assessing petroleum potential (risking) and for basin modeling.
Studies of coals and organic-rich shales have demonstrated that acoustic velocity is reduced by the presence of organic material, that changes in velocity are proportional to the volume of organic material present, and that increases in thermal maturity (largely a function of burial depth and temperature) are accompanied by increases in acoustic velocity (decreases in transit time, Δt). Total organic carbon and the level of organic thermal maturation, expressed in terms of vitrinite reflectance, are two key parameters used for determining the potential of a formation to source hydrocarbons, and each can be mathematically related to Δt. Because acoustic velocity is influenced by a number of factors in addition to organic carbon content, a combination of log measurements can provide improved results when other factors do not mask responses. In particular, acoustic-resistivity crossplot techniques (Fig. 3C.31) and log overlays (Fig. 3C.32) have proved successful.
Fig. 3C.31 – Identification of source-rock potential using a sonic transit-time-resistivity crossplot – points plotted above the line are source rocks, while those below are not. R75°F is well-log resistivity corrected to 75° (courtesy of AAPG).
Fig. 3C.32 – Source-rock identification and assessment using an acoustic/resistivity overlay technique. Separation of the acoustic and resistivity curves (Track 2), labeled Δ log R, indicates organic-rich intervals as shown by the core analyses in Tracks 3 and 4. The amount of separation is directly related to the amount of total organic carbon and is a function of thermal maturity (courtesy of AAPG).
Geomechanical ApplicationsRock Mechanical Properties. The determination of a reservoir’s mechanical properties is critical to reducing drilling risk and maximizing well and reservoir productivity. Estimates of rock mechanical properties are central to drilling programs, well placement, and well-completion design. Mechanical properties include the elastic properties (Young’s modulus, shear modulus, bulk modulus, and Poisson’s ratio) and the inelastic properties (fracture gradient and formation strength).
Elasticity is the property of matter that causes it to resist deformation in volume or shape. Hooke’s law describes the behavior of elastic materials and states that for small deformations, the resulting strain is proportional to the applied stress. Stress is the force applied per unit area, and strain is the fractional distortion that results because of the acting force. The modulus of elasticity is the ratio of stress to strain. Depending on the mode of the acting geological force and type of geological media the force is acting upon, three types of deformation can result as well as three elastic moduli that correspond to each type of deformation. Young’s modulus, E, is the ratio of uniaxial compressive (tensile) stress to the resultant strain. Bulk modulus, K, is the change in volume under hydrostatic pressure (i.e., the ratio of stress to strain). (K is the reciprocal of compressibility.) Shear modulus, μ, is the ratio of shearing (torsional) stress to shearing strain. An additional parameter, Poisson’s ratio, σ, is a measure of the geometric change of shape under uniaxial stress. These four elastic parameters are interrelated such that any one can be expressed in terms of two others and can also be expressed in terms of acoustic-wave velocity and density (Table 3C.13).
The data needed to compute mechanical rock properties are compressional and shear velocities (slowness) and density. Shear and compressional velocities are a function of the bulk modulus, shear modulus, and density of the formation being measured. The Vp/Vs ratio, combined with formation density, ρ, is used to calculate Poisson’s ratio, Young’s modulus, the bulk modulus, and the shear modulus. Whenever possible, log-derived, dynamic rock properties should be calibrated to core-derived static (laboratory) properties, because the static measurements more accurately represent the in-situ reservoir mechanical properties. Rock mechanical properties can be determined using conventional empirical charts or computer programs. The elastic moduli and Poisson’s ratio are used in a variety of applications. These applications include predictions of formation strength, well stimulation (fracture pressure and fracture height), borehole and perforation stability, sand production and drawdown limits in unconsolidated formations, coal evaluation, and determining the roof-rock-strength index for underground mining operations. Rock-mechanics applications of modern multipole tools are discussed in the Anisotropy Analysis section following.
Near-Well ImagingAcoustic data acquired using modern array tools can provide high-resolution (0.5 m), microscale "seismic" 2D and 3D images of structural features in the near-borehole region (10 to 15 m). Conventional seismic-processing techniques, including filtering and migration, are used to extract compressional and shear reflections from the acoustic data. The reflections are then used to image geological features near the borehole. This technique allows the imaging of bed boundaries, thin beds (stringers), fractures, and faults in openhole and cased wells (Fig. 3C.33).
Fig. 3C.33 – Near-well acoustic image of a fractured reservoir using direct (compressional) (right panel) and converted (shear) (left panel) acoustic modes. The fractures indicated in the converted mode data are more apparent than in the direct data and also correlate better with fractures seen on the televiewer image (not shown) (courtesy of SPE).
Advanced Data Analysis And Applications
Processing acoustic data downhole as well as at the surface is necessary to transform the raw acoustic signals recorded by modern logging instruments into data suitable for interpretation and analysis. Data processing takes place during acquisition, in the logging tool itself and in the surface acquisition unit, and also in post-acquisition processing at computing centers. There are a variety of sources of noise in the downhole environment that contaminate the recorded acoustic signal: tool ("road") noise, measurement error, reflection and scattering from rough borehole or bed boundaries, mode conversion, and interferences that occur in the downhole environment. The goal of acoustic-data processing is to minimize the data noise while maximizing the petrophysical information. Data preprocessing reduces the influences of these sources, thus allowing extraction of the true formation signal.
Following the rapid theoretical advances in acoustic-wave propagation made during the 1980s and 1990s, significant advances in data processing provided improved quality in slowness measurements and enabled a number of new applications using Stoneley and dipole-shear wave in open and cased holes. The combined interpretation of Stoneley and dipole-shear acoustic measurements with NMR and borehole imaging enhances formation evaluation.
Slowness AnalysisOne of the primary goals in borehole acoustic logging is to obtain formation slowness from array measurements. Accurate slowness analysis is vital to most petrophysical, geophysical, and seismic applications. A variety of techniques are used for computing slowness from array tools (Table 3C.14).
Fig. 3C.34 – Diagram showing the gather of array acoustic data into either common-source subarrays or common-receiver subarrays. These subarrays cover the same depth interval. Data grouping is in the vertical direction for the source gather, and in the horizontal direction for the receiver gather (courtesy of SPWLA).
The objective of the semblance technique is to find the slowness that maximizes the coherence (time domain) among the wave power spectra over the receiver array. The Nth-root stacking technique is very similar to the semblance technique, except that the waveform amplitudes are modified in a different way to produce a so-called "pseudo" power spectrum instead of the true power spectrum.
Although the algorithms used in semblance and Nth-root stacking differ, the processing procedure is similar. A time window of fixed length is set up on each signature in the array. The windows are offset by a specified time interval on the successive signatures and a correlation is computed over these windows. The windows on all but the first signature are then stepped a certain distance out in time and another correlation value is computed. The process continues for the range of anticipated slowness in the well. This process is repeated as the window is progressively stepped on the first signature.
The group of correlation values obtained is known as a correlogram, in which the waveform coherence or correlation values either from semblance or Nth root method are displayed as a function of slowness and arrival time creating a 3-D surface. The correlation values in a range of arrival times are further combined to form a "combined correlogram," which is a projection of the 3-D surface on the slowness axis. The correlogram peak at each measured-depth level is used to obtain the Δt slowness for the wave type. These individual slowness values are then displayed as a continuous log curve. Filtering for the compressional wave or for the dispersive shear wave is normally required before using either method. As a quality check, the slowness results of semblance or Nth-root stack processing can also be plotted in the combined correlogram together with a "computed travel time" which is plotted against the waveform. The computed travel time represents the sum of the transit-time measurements from the transmitter to the receiver along with the mud travel time. Dipole shear-slowness processing must also take into account dispersion of the flexural wave. Depending on the frequency of the flexural mode, a correction may be required to obtain shear-wave slowness.
Resolution Enhancement. The aperture (vertical resolution) of array tools is typically 3.5 ft. This means that a bed must be at least that thickness to measure true acoustic slowness, although such tools can detect (resolve) beds down to 2-ft thickness. To meet the need for the higher resolution necessary for thin-bed evaluation, waveform-matching (multishot) processing techniques use the redundant information contained in overlapping receiver subarrays to improve vertical resolution (Fig. 3C.34). Recently introduced processing techniques reduce the aperture to 0.5 ft and thus achieve a true vertical resolution of 1.0 ft (Fig. 3C.35). However, reducing the aperture makes the measurements more susceptible to noise.
Fig. 3C.35 – Compressional-wave-slowness curves obtained for various configurations (apertures) of possible overlapping subarrays based on an array tool having four transmitters and eight receivers. Note the increasing resolution of the curves with decreasing subarray aperture. Track 9 is a consistency check obtained by averaging the curves to 3.5-ft aperture and overlaying the results (courtesy of SPWLA).
Anisotropy AnalysisFormation anisotropy—the directional variation of physical properties—can be the result of depositional processes (intrinsic) or tectonic processes (stress-induced). Formation anisotropy is evidenced through variations in permeability, rock strength, fractures, and borehole failure. In acoustic/seismic terms, intrinsic anisotropy is structural in nature and is commonly seen as transverse isotropy (TI or vertical transverse isotropy, VTI), in which properties differ in the vertical or horizontal planes, such as in shales or thinly bedded intervals. Stress-induced anisotropy is known as azimuthal anisotropy (or horizontal transverse isotropy, HTI), in which acoustic parameters in a vertical borehole vary with azimuthal orientation, such as the case of fractures parallel to the borehole.
Analyses of in-situ anisotropy (primarily stress induced) are made using direct or derived shear-wave velocity and provide the magnitude and azimuth of anisotropy (i.e., direction of the maximum and minimum horizontal stresses as well as an indication of their difference). Anisotropy analysis has been widely used in solid-earth seismology, geothermal studies, and more recently, in exploration geophysics (see summaries in Crampin and Chastin and Helbig and Thomsen). In the petroleum industry, these results are used in well design and well placement for optimum reservoir drainage, to detect and characterize faults and fractures in openhole and cased hole, to predict borehole instability and sand production, and for optimizing the design and evaluation of well completions (perforations and hydraulic fracturing) (see following section on Crossed-Dipole Anisotropy Analysis).
Compressional-, shear- and Stoneley-wave properties are each affected, to some degree, by the presence and type of formation anisotropy. While the shear-wave response to azimuthal anisotropy (see the following section on Crossed-Dipole Anisotropy Analysis) and Stoneley waves to VTI anisotropy is well known, the effect on compressional-wave energy is less well characterized. In some situations, information from additional measurements, such as borehole images, dip logs, or both, may be necessary for relating the measured anisotropy to geological features.
Borehole acoustic azimuthal and VTI anisotropy analysis is an advance made possible by the recent introduction of new inversion methods. These methods use the crossed-dipole shear to derive azimuthal anisotropy and the Stoneley wave to derive TI anisotropy in slow formations, or a combination of these modes in deviated wells. A reasonable shear velocity can be derived using inversion techniques with low-frequency Stoneley-wave dispersion which is sensitive to the horizontal shear (in contrast to the dipole’s sensitivity to the vertical shear).
Several new applications have been made possible by anisotropy analysis: identification of formation alteration using dipole-shear dispersion, stress estimation, and distinguishing between intrinsic and stress-induced anisotropy using dispersion crossover. Anisotropy analysis can also be conducted using shear-wave parameters derived from Stoneley-wave dispersion.
Crossed-Dipole Anisotropy Analysis. In anisotropic media, shear waves (both monopole refracted and dipole flexural) split into orthogonally polarized components having different velocities. This is known as shear-wave splitting, shear-wave birefringence, or shear-wave velocity anisotropy. The difference in fast and slow shear-wave slowness provides a measure of the magnitude of anisotropy. Shear-wave splitting is useful for evaluating fractures, faults, bedding planes that intersect the borehole at an angle, and unbalanced tectonic stresses perpendicular to the borehole.
Because the flexural waves induced by dipole-acoustic sources have a directional component, the use of mutually perpendicular (crossed) pairs of dipole transmitters and receivers can detect and measure dipole shear-wave splitting. In isotropic formations, both receiver components (X and Y) will measure the same shear-wave arrival time. However, in an azimuthally anisotropic formation, the shear wave measured at the X-component of the receiver pair will be different from the one measured at the Y-component; one is called the fast shear component, the other, the slow shear component (Fig. 3C.36). To measure the fast and slow shear-wave slowness, four components are measured: two inline components, X-X and Y-Y, and two crossline components, XY and YX (Fig. 3C.37).
Fig. 3C.36 – Diagram illustrating the principle and configuration of crossed-dipole logging: two orthogonal dipole transmitter and receiver array systems. The tool acquires four array data sets, two in-line, XX and YY; two crossed-line, XY and YX, and uses the physics of shear-wave splitting in azimuthally anisotropic rocks to determine the fast-shear polarization azimuth (courtesy of Baker Atlas).
Initially, mathematical rotation methods originally developed for use in processing surface seismic data were used to determine the fast-wave direction together with the slowness values for the fast and slow shear wave. More recently, inversion methods are being applied to the four data sets to simultaneously determine the azimuth and magnitude of the anisotropy (Fig. 3C.38).
The stresses operating on a rock formation are described by a triaxial coordinate system that consists of two principal horizontal stresses, σx and σy, and a vertical stress component, σz, which is the overburden. In a borehole, these downhole stresses are expressed as radial components at the borehole wall: the vertical component, σz, the radial component, σr, and the tangential component, σθ, and the (azimuthal) shear component, σrθ.
Unbalanced-formation-stress components produce distortion around the borehole (stress-induced anisotropy) that results in shear-wave splitting. The azimuth of the fast-shear wave parallels the direction of maximum horizontal stress and the azimuth of the slow-shear wave parallels the direction of the minimum horizontal stress. Hydraulic-fracture azimuth is parallel to the direction of maximum horizontal stress.
Hydraulic well stimulation consists of perforating an interval, then packing and pressuring these perforations to create fractures behind casing to allow increased production. Crossed-dipole anisotropy logging can estimate the vertical extent of the formation-stress fracture along the borehole and its azimuth in the formation (Fig. 3C.40).
Attenuation AnalysisAcoustic-wave attenuation correlates with a variety of petrophysical parameters, including formation lithology and pore-fluid permeability, the degree of fluid saturation, and fractures type. A recent study suggests that the combination of compressional- and shear-wave attenuation logs may provide a potential formation-evaluation tool (Figs. 3C.41 and 3C.42).
Fig. 3C.41 – Comparison between shear-wave attenuation logs and core permeability in a permeable oil zone. The two attenuation curves agree in high-permeability zones (courtesy of SEG).
Fig. 3C.42 – Correlation between compressional- and shear-wave attenuation and fracture location. The interval of 2920 to 2930 m in which P- and S-wave attenuation curves overlay corresponds to the interval of high acoustic reflectivity (Track 4) and fractures identified by the down-going Stonely-wave reflectivity (Track 5) (courtesy of SEG).
Stonely-Wave AnalysisPermeability. Stoneley-wave velocity and attenuation are sensitive to formation and fracture permeability, particularly at low frequencies. Stoneley-wave velocity decreases, and its attenuation increases, as permeability increases. Initial efforts (begun in the 1970s) to derive permeability information from Stoneley data were unsuccessful because neither the necessary low-frequency tools nor the appropriate processing methods had been developed. The parallel development of modern multipole array tools and sophisticated semblance- and inversion-processing methods enable computation of continuous profiles of formation permeability from monopole Stoneley-wave data. Typically, these methods first model the nonpermeability effects using the elastic-wave theory and then relate differences between the modeled and the measured data to formation permeability. One approach to Stoneley-wave processing is comprised of three parts: slowness analysis, reflectance mapping, and permeability estimation (Fig. 3C.43).
Wave-separation processing minimizes the effects of nonpermeability-related influences (e.g., road noise and borehole scattering) and yields reflectance logs for the direct and reflected Stoneley-wave data. The center-frequency log for the reflected wave data characterizes the Stoneley-wave attenuation and can be used to indicate fractures, vugs, and bed boundaries. The center-frequency log for the direct (transmitted) data is used to estimate formation permeability. Knowledge of the formation-fluid properties (viscosity and compressibility from core or NMR) enables quantitative estimates. Without this information log-derived permeability estimates are only qualitative.
These models require sophisticated computer processing. A simplified, field-oriented technique based on Stoneley amplitude has so far provided good results in ideal conditions and when calibrated to core or NMR data.
Stoneley-derived formation permeability compares well with estimates obtained by nonacoustic methods, including core analysis and NMR logs (Fig. 3C.44). Additional improvement in permeability estimation is possible when anisotropy information is incorporated in the process. Formation evaluation is further enhanced when these diverse measurements are integrated through joint interpretation. As indicated by the plots in Fig. 3C.44, calibration of log responses to core data would improve log-predicted values in noncored intervals.
Gas Detection. Stoneley-wave properties, used in combination with other data, facilitate identification of gas-bearing intervals. A gas-saturated formation interval has drastically different fluid mobility and compressibility compared to those of the surrounding formations. In the presence of gas, Stoneley-derived permeability is overestimated because of increased pore-fluid mobility (decreased viscosity) and compressibility, and NMR permeability may be underestimated because of a decreased hydrogen index. If first calibrated in a nongas-bearing interval, separation of the two permeability curves indicates the presence of gas (Fig. 3C.45). The plot also shows the neutron-density crossplot is effective in detecting the gas zone and may be more sensitive in some intervals.
Fig. 3C.46 – Gas detection in a high-porosity, laminated formation using thin-bed slowness and Stoneley reflection. The increase in P-wave slowness (Track 2) across the interval of thin gas ssands, X740 to X800 m, is accompanied by low density values (Track 2) and numerous Stoneley-wave reflections (Track 3) (courtesy of Baker Atlas).
The Stoneley wave responds to fracture permeability while NMR does not, thus, the Stoneley-wave permeability is greater in the presence of fractures than NMR permeability. The combined use of Stoneley- and NMR-derived permeability provides a fracture indicator in both sandstones (Fig. 3C.47) and carbonates (Fig. 3C.48).
Fig. 3C.47 – Stonely-wave and NMR permeability profiles in the presence of borehole fractures. Stoneley-wave results are dominated by the fracture-system contribution, while the NMR results are dominated by the matrix-based permeability. The fractured sandstone interval, indicated by separation of the two permeability curves (Tracks 2 and 6), is confirmed by Stoneley-wave reflection data (Track 3) and the acoustic image log (Track 5) (courtesy of Baker Atlas).
Fig. 3C.48 – Correlation between azimuthal anisotropy and Stoneley reflections. A fracture zone in a U.S. mid-continent well causes both shear-wave splitting and Stoneley-wave reflections – in the anisotropy map (left), brighter colors represent higher anisotropy. The azimuth of the fast shear arrival is east/west (center). The fractures intersecting the borehole cause significant up- and down-going Stoneley-wave reflections seen in Stoneley-waveform amplitudes (right) (courtesy of Baker Atlas).
This chapter has summarized the characteristics and capabilities of acoustic-logging tools. Although determination of geophysical (seismic) and petrophysical (porosity) properties have traditionally been the most widely used applications for acoustic-log data, new and advanced applications in geomechanics (anisotropy, fracture, and stress evaluation) and petrophysics (permeability) have been made possible by the multiarray tools and inversion processing methods now available. Interpretation methods that combine acoustic data with core and other log data, such as resistivity, nuclear, NMR, and borehole images, facilitate determination of fluid saturation and typing when logging conditions provide reliable data. For all applications of acoustic logging careful planning is critical to ensure successful acquisition of the desired data.
|C||=||shale compaction coefficient|
|Cp||=||compaction correction factor|
|R1||=||receiver one in a two-receiver tool configuration|
|R2||=||receiver two in a two-receiver tool configuration|
|Sw||=||water saturation, fraction|
|t1||=||travel time between T2 and R1|
|t2||=||travel time between T2 and R2|
|t3||=||travel time between T1 and R1|
|t4||=||travel time between T1 and R2|
|TR||=||transmitter-to-receiver spacing (ft)|
|T1||=||transmitter one in a two-transmitter tool (BHC) configuration|
|T2||=||transmitter two in a two-transmitter tool (BHC) configuration|
|v||=||velocity of the formation|
|x||=||distance between R1 and R2|
|X12||=||distance between R1 and R2|
|α||=||porosity correction factor|
|Δt||=||transit time, sec|
|Δtc||=||compressional travel time|
|ϕ||=||fractional porosity of rock|
|ρ||=||formation density, lbm/ft3|
|σr||=||borehole radial stress|
|σx||=||principal horizontal stress x-direction|
|σy||=||principal horizontal stress y-direction|
|σz||=||principal vertical stress z-direction|
|σθ||=||borehole tangential stress|
- Hornby, B.E. 1995. Use of Full-Wave Sonic Data to Image Near-Borehole Structural Features. Petroleum Geoscience 1 (2): 109–114.
- Coates, R., Kane, M., Chang, C. et al. 2000. Single-well Sonic Imaging: High-Definition Reservoir Cross-sections from Horizontal Wells. Presented at the SPE/CIM International Conference on Horizontal Well Technology, Calgary, Alberta, Canada, 6-8 November 2000. SPE-65457-MS. http://dx.doi.org/10.2118/65457-MS.
- Cheng, C.H., Paillet, F.L., and Pennington, W.D. 1992. Acoustic-Waveform Logging--Advances In Theory and Application. The Log Analyst 33 (3): 239. SPWLA-1992-v33n3a2.
- Paillet, F.L. and Cheng, C.H. 1991. Acoustic Waves in Boreholes, 1–264. Boca Raton, Florida: CRC Press.
- Mavko, G., Mukerji, T., and Dvorkin, J. 1998. The Rock Physics Handbook—Tools for Seismic Analysis in Porous Media, 1-329. Cambridge, England: Cambridge University Press.
- Hearst, J.R., Nelson, P.H., and Paillet, F.L. 2000. Acoustic Logging, Well Logging for Physical Properties, second edition, Chap. 8, 257–303. New York City: Wiley and Sons Inc.
- Wu, P.T., Darling, H.L., and Scheibner, D. 1995. Low-Frequency P-Wave Logging for Improved Compressional Velocity in Slow Formation Gas Zones, paper BG1.3, Expanded Abstracts, 1995 Annual Meeting Technical Program, SEG, 9–12.
- Tang, X.-M. and Cheng, A. 2004. Quantitative Borehole Acoustic Methods. In Handbook of Geophysical Exploration, Seismic Exploration, Vol. 24. Oxford, England: Pergamon Press (Elsevier).
- Cheng, C.H. and Toksoz, M.N. 1983. Determination of Shear Wave Velocities in ‘Slow’ Formation, paper V, Trans., 1983 Annual Logging Symposium, Soc. of Professional Well Log Analysts (SPWLA) 1–18.
- Valero, H.P., Peng, L., Yamamoto, M. et al. 2004. Processing of Monopole Compressional In Slow Formation. Presented at the SEG Annual Meeting, Denver, Colorado, 10-15 October.
- Carmichael, R.S. ed. 1982. Handbook of Physical Properties of Rocks, Vol. 2, 1-228. Boca Raton, Florida: CRC Press Inc.
- Leonardon, E.G. 1961. Logging, Sampling, and Testing. In History of Petroleum Engineering, Ch. 8, 493-578. Dallas, Texas: American Petroleum Inst., Div. of Production.
- Snyder, D.D. and Fleming, D.B. 1985. Well-Logging—A Twenty-Five-Year Perspective. Geophysics 50 (12): 2,504–2,529. http://dx.doi.org/10.1190/1.1441881.
- Jorden, J.R. and Campbell, F.L. 1986. Well Logging II—Electric and Acoustic Logging, Vol. 10, 95-151. Richardsson, Texas: SPE Monograph Series.
- Timur, A. 1987. Acoustic Logging. Petroleum Engineering Handbook, H.B. Bradley ed., Ch. 7. Richardson, Texas: SPE.
- Liu, O.Y. 1987. The Sources of Errors in Slowness Measurements and an Evaluation of Full Waveform Compensation Techniques. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 27-30 September 1987. SPE-16772-MS. http://dx.doi.org/10.2118/16772-MS.
- Aron, J. and Murray, J. 1978. Formation Compressional and Shear Interval Transit Time Logging by Means of Long Spacings and Digital Techniques. Presented at the SPE Annual Fall Technical Conference and Exhibition, Houston, Texas, 1-3 October 1978. SPE-7446-MS. http://dx.doi.org/10.2118/7446-MS.
- Koerperich, E.A. 1980. Shear Wave Velocities Determined From Long-and Short-Spaced Borehole Acoustic Devices. SPE J 20 (5): 317-326. SPE-8237-PA. http://dx.doi.org/10.2118/8237-PA.
- Morris, C.F., Little, T.M., and W., L.I. 1984. A New Sonic Array Tool for Full Waveform Logging. Presented at the SPE Annual Technical Conference and Exhibition, Houston, Texas, 16-19 September 1984. SPE-13285-MS. http://dx.doi.org/10.2118/13285-MS.
- Williams, D.M. et al. 1984. The Long Spaced Acoustic Logging Tool, paper T. Trans., 1984 Annual Logging Symposium, SPWLA, 1–15.
- Chen, S.T. and Eriksen, E.A. 1991. Compressional- and Shear-Wave Logging in Open and Cased Holes Using a Multipole Tool. Geophysics 56 (4): 550–557. http://dx.doi.org/10.1190/1.1443071.
- Georgi, D.T. et al. 1991. Application of Shear and Compressional Transit-Time Data to Cased-Hole Carbonate Reservoir Evaluation. The Log Analyst 32 (3): 129.
- Tang, X., and Patterson, D. 2005. Analyzing and Processing Acoustic Logging Data for Poorly Bonded Cased Boreholes, paper OO. Trans., 2005 Annual Logging Symposium, SPWLA, 1–11.
- Stevens, J.L. and Day, S.M. 1986. Shear Velocity Logging in Slow Formations Using the Stoneley Wave. Geophysics 51 (1): 137–147. http://dx.doi.org/10.1190/1.1442027.
- Liu, O.Y. 1984. Stoneley Wave-Derived Delta-T Shear Log, paper ZZ. Trans., 1984 Annual Logging Symposium, SPWLA, 1–14.
- White, J.E. 1967. The HULA Log—A Proposed Acoustic Tool, paper I. Trans., 1967 Annual Logging Symposium, SPWLA, 1–30.
- Zemanek, J. et al. 1984. Continuous Acoustic Shear Wave Logging, paper U. Trans., 1984 Annual Logging Symposium, SPWLA, 1–14
- Kessler, C. et al. 2004. New Development in Monopole Acoustic Logging, paper Z. Trans., 2004 Annual Logging Symposium, SPWLA, 1–12.
- Harrison, A.R., Randall, C.J., Aron, J.B. et al. 1990. Acquisition and Analysis of Sonic Waveforms From a Borehole Monopole and Dipole Source for the Determination of Compressional and Shear Speeds and Their Relation to Rock Mechanical Properties and Surface Seismic Data. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 23-26 September 1990. SPE-20557-MS. http://dx.doi.org/10.2118/20557-MS.
- Tello, L.N., Blankinship, T.J., Roberts, E.K. et al. 1999. A Dipole Array Sonic Tool for Vertical and Deviated Wells. Presented at the SPE Annual Technical Conference and Exhibition, Houston, Texas, 3-6 October 1999. SPE-56790-MS. http://dx.doi.org/10.2118/56790-MS.
- Oden, C.P., Stowell, J.R., and LoCoco, J.J. 2000. Variable Frequency Monopole-Dipole Sonic Logging for Shear Velocity—Applications and Test Results. Proc., 2000 Intl. Symposium on Borehole Geophysics for Minerals, Geotechnical, and Groundwater Applications, Minerals, and Geotechnical Logging Soc., SPWLA, 71–76.
- Kessler, C. and Varsamis, G.L. 2001. A New Generation Crossed Dipole Logging Tool: Design and Case Histories. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 30 September-3 October 2001. SPE-71740-MS. http://dx.doi.org/10.2118/71740-MS.
- Pistre, V. et al. 2005. A Modular Wireline Sonic Tool for Measurements of 3D (Azimuthal, Radial, and Axial) Formation Acoustic Properties, paper P. Trans., 2005 Annual Logging Symposium, SPWLA, 1–13.
- Plona, T. et al. 2002. Mechanical Damage Detection and Anisotropy Evaluation Using Dipole Sonic, paper F. Trans., 2002 Annual Logging Symposium, SPWLA, 1–14.
- Sinha, B.K., Vissapragada, B., Kisra, S. et al. 2005. Optimal Well Completions Using Radial Profiling of Formation Shear Slownesses. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 9-12 October 2005. SPE-95837-MS. http://dx.doi.org/10.2118/95837-MS.
- Aron, J. et al. 1994. Sonic Compressional Measurements While Drilling, paper SS. Trans., 1994 Annual Logging Symposium, SPWLA, 1–17.
- Minear, J. et al. 1995. Compressional Slowness Measurements While Drilling, paper VV. Trans., 1995 Annual Logging Symposium, SPWLA, 1–12.
- Hsu, K. et al. 1997. Interpretation and Analysis of Sonic While Drilling Data in Overpressured Formations, paper FF. Trans., 1997 Annual Logging Symposium, SPWLA, 1–14.
- Hashem, M. et al. 1999. Seismic Tie Using Sonic-While-Drilling Measurements, paper I. Trans., 1999 Annual Logging Symposium, SPWLA, 1–13.
- Market, J., Schmitt, D., and Deady, R. 2004. LWD Sonic Logging in Cased Holes, paper U. Trans., 2004 Annual Logging Symposium, SPWLA, 1–14.
- Varsamis, G. et al. 1999. A New MWD Full Wave Dual Mode Sonic Tool Design and Case Histories, paper F. Trans., 1999 Annual Logging Symposium, SPWLA, 1–12.
- Varsamis, G.L. et al. 2000. LWD Shear Velocity Logging in Slow Formations—Design Decisions and Case Histories, paper O. Trans., 2000 Annual Logging Symposium, SPWLA, 1–13.
- Joyce, B. et al. 2001. Introduction to a New Omni-Directional Acoustic System for Improved Real Time LWD Sonic Logging—Tool Design and Field Test Results, paper SS. Trans., 2001 Annual Logging Symposium, SPWLA, 1–14.
- III, J.V.L., Dubinsky, V., Patterson, D. et al. 2001. Field Test Results Demonstrating Improved Real-Time Data Quality in an Advanced LWD Acoustic System. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 30 September-3 October 2001. SPE-71732-MS. http://dx.doi.org/10.2118/71732-MS.
- Tang, X. et al. 2003. Shear-Velocity Measurement in the Logging-While-Drilling Environment—Modeling and Field Evaluations. Petrophysics 44 (2): 79–90.
- Market, J. et al. 2002. Processing and Quality Control of LWD Dipole Sonic Measurements, paper PP. Trans., 2002 Annual Logging Symposium, SPWLA, 1–14.
- Dubinsky, V., Tang, X.M., Bolshakov, A. et al. 2003. Engineering Aspects of LWD Quadrupole Measurements and Field Test Results. Presented at the SPE Annual Technical Conference and Exhibition, Denver, Colorado, 5-8 October 2003. SPE-84248-MS. http://dx.doi.org/10.2118/84248-MS.
- Goldberg, D. et al. 2003. Analysis of LWD Sonic Data in Low-Velocity Formations, paper BH 1.5. Expanded Abstracts, 2003 Annual Meeting Technical Program, SEG, 301–304.
- Boonen, P. and Yogeswaren, E. A Dual-Frequency LWD Sonic Tool Expands Existing Unipolar Transmitter Technology to Supply Shear Wave Data in Soft Formations, paper X. Trans., 2004 Annual Logging Symposium, SPWLA, 1–11.
- Endo, T., Yoneshima, S., and Valero, H.-P. 2004. Analysis of LWD Sonic Data in Very Slow Formation—Compressional Processing and Indirect Vp/Vs Estimation, paper G. Proc., 2004 Formation Evaluation Symposium of Japan, SPWLA Japan Chapter, 1–8.
- Tang, X., Zheng, Y., and Dubinsky, V. 2005. Logging While Drilling Acoustic Measurement in Unconsolidated Slow Formations, paper R. Trans., 2005 Annual Logging Symposium, 1–13.
- Zhu, Z. et al. 2005. Experimental Studies of Multipole Acoustic Logging with Scaled Borehole, paper BG 2.6. Expanded Abstracts, 2005 Annual Meeting Technical Program, SEG, 376–379.
- Tang, X. et al. 2003. Logging-While-Drilling Shear and Compressional Measurements in Varying Environments, paper II. Trans., 2003 Annual Logging Symposium, SPWLA, 1–13.
- Wang, T., and Tang, X. 2003. Investigation of LWD Quadrupole Shear Measurement in Real Environments, paper KK. Trans., 2003 Annual Logging Symposium, SPWLA, 1–12.
- Haugland, S.M. 2004. Frequency Dispersion Effects on LWD Shear Sonic Measurements in Acoustically Slow Environments. Presented at the SPE Annual Technical Conference and Exhibition, Houston, Texas, 26-29 September 2004. SPE-90505-MS. http://dx.doi.org/10.2118/90505-MS.
- Huang, X., Zheng, Y., and Toksoz, N.M. 2004. Effects of Tool Eccentricity on Acoustic Logging While Drilling (LWD) Measurements, paper BG 1.4. Expanded Abstracts, 2004 Annual Meeting Technical Program, SEG, 290–293.
- Underhill, W., Esmersoy, C., Hawthorn, A. et al. 2001. Demonstrations of Real-Time Borehole Seismic From an LWD Tool. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 30 September-3 October. SPE-71365-MS. http://dx.doi.org/10.2118/71365-MS.
- Haldorsen, J.B.U., Esmersoy, C., Hawthorn, A. et al. 2003. Optimizing the Well Construction Process: Full-Waveform Data From While-Drilling Seismic Measurements in the South Caspian Sea. Presented at the SPE/IADC Drilling Conference, Amsterdam, Netherlands, 19-21 February 2003. SPE-79844-MS. http://dx.doi.org/10.2118/79844-MS.
- Esmersoy, C., Hawthorn, A., Durrand, C., and Armstrong, P. 2005. Seismic MWD: Drilling in time, on time, it's about time. The Leading Edge 24(1): 56–62. http://dx.doi.org/10.1190/1.1859702.
- Althoff, G., Cornish, B., Varsamis, G. et al. 2004. New Concepts for Seismic Surveys While Drilling. Presented at the SPE Annual Technical Conference and Exhibition, Houston, Texas, 26-29 September 2004. SPE-90751-MS. http://dx.doi.org/10.2118/90751-MS.
- Walker, T. 1968. Origin of the ‘W’ Pattern on Cased Hole Micro-Seismogram Logs. The Log Analyst 9 (2): 30–32.
- Zemanek, J., Caldwell, R.L., Glenn Jr., E.E. et al. 1969. The Borehole TeleviewerA New Logging Concept for Fracture Location and Other Types of Borehole Inspection. J Pet Technol 21 (6): 762-774. SPE-2402-PA. http://dx.doi.org/10.2118/2402-PS. Cite error: Invalid
<ref>tag; name "r61" defined multiple times with different content
- Prensky, S. 1999. Advances in Borehole Imaging Technology and Applications. In Borehole Imaging—Applications and Case Histories, M.A. Lovell, G. Williamson, and P.K. Harvey eds., 1-44. London: Geological Soc., Special Publication No. 159.
- Bigelow, E.L. 1990. Cement Evaluation, 1-142. Houston: Atlas Wireline Services.
- Hill, A.D. 1990. Cement-Quality Logging. Production Logging—Theoretical and Interpretative Elements, Vol. 14, 123-140. Richardson, Texas: Monograph Series, SPE.
- Jutten, J. and Morriss, S.L. 1990. Cement Job Evaluation. Well Cementing, E.B. Nelson ed., Ch. 16, 16-1-16-44. Amsterdam: Elsevier, Developments in Petroleum Science No. 28.
- Rouillac, D. 1994. Cement Evaluation Logging Handbook, 1-176. Paris: Editions Technip.
- Smolen, J.J. 1996. Cased Hole and Production Log Evaluation, Ch. 10 and 11, 161-213. Tulsa, Oklahoma: PennWell Publishing Co.
- Baker, L.J. 1984. The Effect of the Invaded Zone on Full Wavetrain Acoustic Logging. Geophysics 49 (6): 796–809. http://dx.doi.org/10.1190/1.1441708.
- Alberty, M. 1996. The Influence of the Borehole Environment upon Compressional Sonic Logs. The Log Analyst 37 (4): 30–44.
- Chi, S., Torres-Verdín, C., Wu, J. et al. 2006. Assessment of mud-filtrate invasion effects on borehole acoustic logs and radial profiling of formation elastic parameters. SPE Res Eval & Eng 9 (5): 553–564. SPE-90159-PA. http://dx.doi.org/10.2118/90159-PA.
- Pilkington, P.E. 1992. Cement Evaluation - Past, Present, and Future. J Pet Technol 44 (2): 132-140. SPE-20314-PA. http://dx.doi.org/10.2118/20314-PA.
- Hayman, A.J., Hutin, R., and Wright, P.V. 1991. High-Resolution Cementation and Corrosion Imaging by Ultrasound, paper KK. Trans., 1991 Annual Logging Symposium, SPWLA, 1–25.
- Graham, W.L., Silva, C.I., Leimkuhler, J.M. et al. 1997. Cement Evaluation and Casing Inspection With Advanced Ultrasonic Scanning Methods. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 5-8 October 1997. SPE-38651-MS. http://dx.doi.org/10.2118/38651-MS.
- Froelich, B., Dumont, A., Pittman, D. et al. 1982. Cement Evaluation Tool: A New Approach to Cement Evaluation. J Pet Technol 34 (8): 1835-1841. SPE-10207-PA. http://dx.doi.org/10.2118/10207-PA.
- Sheives, T.C., Tello, L.N., Maki Jr., V.E. et al. 1986. A Comparison of New Ultrasonic Cement and Casing Evaluation Logs With Standard Cement Bond Logs. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 5-8 October 1986. SPE-15436-MS. http://dx.doi.org/10.2118/15436-MS.
- Bigelow, E.L., Domangue, E.J., and Lester, R.A. 1990. A New and Innovative Technology for Cement Evaluation. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 23-26 September 1990. SPE-20585-MS. http://dx.doi.org/10.2118/20585-MS.
- Schmidt, M.G. 1989. The Micro CBL—A Second Generation Radial Cement Evaluation Instrument, paper Z. Trans., 1989 Annual Logging Symposium, SPWLA, 1–16.
- Tello, L.N. et al. 1995. New Efficient Method for Radial Cement Bond Evaluation. Canadian Well Logging Soc. J. 20: 60–69.
- Butsch, R.J. 1995. Overcoming Interpretation Problems of Gas-Contaminated Cement Using Ultrasonic Cement Logs. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 22-25 October 1995. SPE-30509-MS. http://dx.doi.org/10.2118/30509-MS.
- Frisch, G.J., Graham, W.L., and Griffith, J. 1999. Assessment of Foamed - Cement Slurries Using Conventional Cement Evaluation Logs and Improved Interpretation Methods. Presented at the SPE Rocky Mountain Regional Meeting, Gillette, Wyoming, 15-18 May 1999. SPE-55649-MS. http://dx.doi.org/10.2118/55649-MS.
- Frisch, G., Graham, L., and Griffith, J. 2000. A Novel and Economical Processing Technique Using Conventional Bond Logs and Ultrasonic Tools for Enhanced Cement Evaluation, paper EE. Trans., 2000 Annual Logging Symposium, SPWLA, 1–14.
- Frisch, G.J., Fox, P.E., Hunt, D.A. et al. 2005. Advances in Cement Evaluation Tools and Processing Methods Allow Improved Interpretation of Complex Cements. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 9-12 October 2005. SPE-97186-MS. http://dx.doi.org/10.2118/97186-MS.
- Kuijk, R.v., Zeroug, S., Froelich, B. et al. 2005. A Novel Ultrasonic Cased-Hole Imager for Enhanced Cement Evaluation. Presented at the International Petroleum Technology Conference, Doha, Qatar, 21-23 November 2005. IPTC-10546-MS. http://dx.doi.org/10.2523/10546-MS.
- Dumont, A., Patin, J.-B., and Le Floch, G. 1984. A Single Tool for Corrosion and Cement Evaluation. Presented at the SPE Annual Technical Conference and Exhibition, Houston, Texas, 16-19 September 1984. SPE-13140-MS. http://dx.doi.org/10.2118/13140-MS.
- Roberts, D.L. and Richards, J.W. 1987. Casing Corrosion Evaluation Using Ultrasonic Techniques: A Unique Approach for West Texas Wells. Proc., 1987 Annual Short Course, Southwestern Petroleum Short Course Association, Lubbock, Texas, 20–33.
- Frisch, G. and Mandal, B. 2001. Advanced Ultrasonic Scanning Tool and Evaluation Methods Improve and Standardize Casing Inspection. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 30 September-3 October 2001. SPE-71399-MS. http://dx.doi.org/10.2118/71399-MS.
- Fertl, W.H. 1981. Openhole Crossplot Concepts A Powerful Technique in Well Log Analysis. J Pet Technol 33 (3): 535-549. SPE-8115-PA. http://dx.doi.org/10.2118/8115-PA.
- Bigelow, E.L. 1995. Introduction to Wireline Log Analysis, 1-312. Houston: Baker Atlas.
- Boyer, S. and Mari, J.-L. 1997. Seismic Surveying and Well Logging Oil and Gas Exploration Techniques, 1-192. Paris: Editions Technip. Cite error: Invalid
<ref>tag; name "r90" defined multiple times with different content
- Pickett, G.R. 1963. Acoustic Character Logs and Their Applications in Formation Evaluation. J Pet Technol 15 (6): 659-667. SPE-452-PA. http://dx.doi.org/10.2118/452-PA.
- Medlin, W.L. and Alhllall, K.A. 1992. Shear-Wave Porosity Logging in Sands. SPE Form Eval 7 (1): 106-112. SPE-20558-PA. http://dx.doi.org/10.2118/20558-PA.
- Tixier, M.P., Alger, R.P., and Doh, C.A. 1959. Sonic Logging. J Pet Technol Trans., AIME, 216: 106.
- Wyllie, M.R.J., Gregory, A.R., and Gardner, L.W. 1956. Elastic Wave Velocities in Heterogeneous and Porous Media. Geophysics 21 (1): 41–70. http://dx.doi.org/10.1190/1.1438217.
- Raymer, L.L., Hunt, E.R., and Gardner, J.S. 1980. An Improved Sonic Transit Time-to-Porosity Transform, paper P. Trans., 1980 Annual Logging Symposium, SPWLA, 1–12.
- Dvorkin, J. and Nur, A. 1998. Time-Average Equation Revisited. Geophysics 63 (2): 460–464. http://dx.doi.org/10.1190/1.1444347.
- Krief, M. et al. 1990. A New Petrophysical Interpretation Using the Velocities of P and S Waves (Full Waveform Sonic). The Log Analyst 31 (6): 355–369.
- Spears, R., and Nicosia, W.: "Porosity Estimation from Shear Wave Interval Transit Time in the Norphlet Aeolian Jurassic Sandstone of Southern and Offshore Alabama," paper N, Trans., 2003 Annual Logging Symposium, SPWLA, 1–12.
- Valero, H.-P., Skelton, O., and Almeida, C.M. 2003. Processing of Monopole Sonic Waveforms Through Cased Hole, paper BH 1.1. Expanded Abstracts, 2003 Annual Meeting Technical Program, SEG, 285–288.
- Tello, L.N. et al. 2004. Through-Tubing Hostile-Environment Acoustic Logging Tool, paper Y. Trans., 2004 Annual Logging Symposium, SPWLA, 1–11.
- Smith, J., Fisburn, T., and Bigelow, E. 1996. Acoustic Porosity Corrections for Gas and Light Hydrocarbon-Bearing Sandstones, paper AAA. Trans., 1996 Annual Logging Symposium, SPWLA, 1–9.
- Knackstedt, M.A., Arns, C.H., and Pinczewski, W.V. 2003. Velocity-Porosity Relationships, Part 1: Accurate Velocity Model for Clean Consolidated Sandstones. Geophysics 68 (6): 1822–1834. http://dx.doi.org/10.1190/1.1635035.
- Knackstedt, M.A., Arns, C.H., and Pincsewski, W.V. 2005. Velocity-Porosity Relationship: Predictive Velocity Model for Cemented Sands Composed of Multiplemineral Phases. Geophysical Prospecting 53 (3): 349–372.
- Kazatchenko, E., Markov, M., and Mousatov, A. 2003. Determination of Primary and Secondary Porosity in Carbonate Formations Using Acoustic Data. Presented at the SPE Annual Technical Conference and Exhibition, Denver, Colorado, 5-8 October 2003. SPE-84209-MS. http://dx.doi.org/10.2118/84209-MS.
- Burke, J.A., Campbell Jr., R.L., and Schmidt, A.W. 1969. The Litho-Porosity Cross Plot a Method of Determining Rock Characteristics for Computation of Log Data. Presented at the SPE Illinois Basin Regional Meeting, Evansville, Indiana, 30-31 October. SPE-2771-MS. http://dx.doi.org/10.2118/2771-MS.
- Clavier, C. and Rust, D.H. 1976. MID Plot: A New Lithology Technique. The Log Analyst 17 (6): 16–24.
- Robert, H.V. and Campbell, R.L. Jr. 1976. Applications of CORIBAND to the Evaluation of Sandstones Containing Mica. The Log Analyst 17 (1): 33–40.
- Nations, J.F. 1974. Lithology and Porosity from Acoustic Shear and Compressional Wave Transit Time Relationships. The Log Analyst 15 (6): 3–8.
- astwood, R.L. and Castagna, J.P. 1987. Interpretation of V p /V s Ratios from Sonic Logs. Shear-Wave Exploration, S.H. Danbom and S.N. Domenico eds., 139-153. Tulsa, Oklahoma: SEG, Geophysical Development Series No. 1.
- Johnston, J.E. and Christensen, N.I. 1993. Compressional to Shear Velocity Ratios in Sedimentary Rocks. Intl. J. of Rock Mechanics and Mining Sciences and Geomechanical Abstracts 30 (7): 751–754.
- Assefa, S., McCann, C., and Sothcott, J. 2003. Velocities of Compressional and Shear Waves in Limestone. Geophysical Prospecting 51 (1): 1–13.
- Hsu, K. and Chang, S.K. 1987. Multiple-Shot Processing of Array Sonic Waveforms for High Resolution Sonic Logs. Geophysics 52 (10): 1,376–1,390. http://dx.doi.org/10.1190/1.1442250.
- Tsuneyama, F. et al. 2003. Vp/Vs Ratio as a Rock Frame Indicator for a Carbonate Reservoir. First Break 21 (7): 53–58.
- Zhang, T., Tang, X.M., and Patterson, D. 2000. Evaluation of Laminated Thin Beds in Formations Using High-Resolution Acoustic Slowness Logs, paper XX. Trans., 2000 Annual Logging Symposium, SPWLA, 1–14.
- Hingle, A.T. Jr. 1959. The Use of Logs in Exploration Problems. Technical Program, 1959 Annual Meeting, SEG, 38.
- Berry, J.E. 1959. Acoustic Velocity in Porous Media. Trans., AIME 216: 262.
- Pickett, G.R. 1966. A Review of Current Techniques for Determination of Water Saturation From Logs. J Pet Technol 18 (11): 1425-1433. SPE-1446-PA. http://dx.doi.org/10.2118/1446-PA.
- Pickett, G.R. 1973. Pattern Recognition as a Means of Formation Evaluation. The Log Analyst 14 (4): 3–11.
- Kamel, M.H., and Mabrouk, W.M. 2002. An Equation for Estimating Water Saturation in Clean Formations Utilizing Resistivity and Sonic Logs: Theory and Application. J. of Petroleum Science and Engineering 36 (3–4): 59–168.
- Souder, W.W. 2002. Using Sonic Logs to Predict Fluid Type. Petrophysics 43 (5): 412–419.
- Chardac, O., Brie, A., and Chouker, A.C. 2003. Correlations of Shear vs. Compressional in Shaly Sands and Application to Quicklook Hydrocarbon Detection. Presented at the SPE Annual Technical Conference and Exhibition, Denver, Colorado, 5-8 October 2003. SPE-84205-MS. http://dx.doi.org/10.2118/84205-MS.
- Hamada, G.M. 2004. Vp/Vs Crossplots Identify Reservoir Fluids. Oil & Gas J. 102 (24): 41–44.
- Khazanehdari, J., and McCann, C. 2005. Acoustic and Petrophysical Relationships in Low-Shale Sandstone. Geophysical Prospecting 53 (4): 447–462.
- Williams, D.M. 1990. The Acoustic Log Hydrocarbon Indicator, paper W. Trans., 1990 Annual Logging Symposium, SPWLA, 1–22.
- Brie, A., Pampuri, F., Marsala, A.F. et al. 1995. Shear Sonic Interpretation in Gas-Bearing Sands. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 22-25 October 1995. SPE-30595-MS. http://dx.doi.org/10.2118/30595-MS.
- Rogen, B. et al. 2005. Ultrasonic Velocities of North Sea Chalk Samples: Influence of Porosity. Geophysical Prospecting 53 (4): 481–496.
- Watson, I.A., Lines, L.R., and Brittle, K.F. 2002. Heavy-Oil Reservoir Characterization Using Elastic Wave Properties. The Leading Edge 21 (8): 736–739. http://dx.doi.org/10.1190/1.1503191.
- Gutierrez, M., Dvorkin, J., and Nur, A. 2000. In-Situ Hydrocarbon Identification and Reservoir Monitoring Using Sonic Logs, La Cira-Infantas Oil Field (Colombia), paper RPB 2.6. Expanded Abstracts, 2000 Annual Meeting Technical Program, SEG, 1727–1730.
- Gutierrez, M.A., Dvorkin, J., and Nur, A. 2001. Theoretical Rock Physics for Bypassed Oil Detection Behind the Casing—La Cira-Infantas Oil Field. The Leading Edge 20 (2): 192–198. http://dx.doi.org/10.1190/1.1438910.
- Batzle, M. et al. 2004. Heavy Oils: Seismic Properties, paper RP P2.2. Expanded Abstracts, 2004 Annual Meeting Technical Program, SEG, 1762–1765.
- Thomas, D.H. 1978. Seismic Applications of Sonic Logs. The Log Analyst 19 (1): 23–32.
- Liner, C.L. Interpreting Seismic Data. 2000. Exploring for Oil and Gas Traps, Treatise of Petroleum Geology, Handbook of Petroleum Geology E.A. Beaumont, and N.H. Foster eds., Chap. 12, 12-13–12–17. Tulsa, Oklahoma: AAPG.
- Bork, J. and Wood, L.C. 2001. Seismic Interpretation of Sonic Logs, paper INT 1.4. Expanded Abstracts, 2001 Annual Meeting Technical Program, SEG, 510–513.
- Box, R., and Lowrey, P. 2003. Reconciling Sonic Logs with Check-Shot Surveys; Stretching Synthetic Seismograms. The Leading Edge 22 (6): 510–517. http://dx.doi.org/10.1190/1.1587672.
- White, R.E. 2003. Tying Well-Log Synthetic Seismograms to Seismic Data: The Key Factors, paper W-3.7. Expanded Abstracts, 2003 Annual Meeting Technical Program, SEG, 2449–2452.
- Walls, J., Dvorkin, J., and Carr, M. 2004. Well Logs and Rock Physics in Seismic Reservoir Characterization. Presented at the Offshore Technology Conference, Houston, Texas, 3-6 May. OTC-16921-MS. http://dx.doi.org/10.4043/16921-MS.
- Bruce, B. and Bowers, G. 2002. Pore Pressure Terminology. The Leading Edge 21 (2): 170–173. http://dx.doi.org/10.1190/1.1452607.
- Fertl, W.H., Chapman, R.E., and Holz, R.F. eds. 1994. Studies in Abnormal Pressure, 1-454. Amsterdam: Elsevier, Developments in Petroleum Science No. 38.
- Bowers, G.L. 2002. Detecting High Pressure. The Leading Edge 21 (2): 174–177. http://dx.doi.org/10.1190/1.1452608.
- Chilingar, G.V., Serebryakov, V.A. and Robertson, J.O. Jr. eds. 2002. Origin and Prediction of Abnormal Formation Pressures, 1-390. Amsterdam: Elsevier, Developments in Petroleum Science No. 50.
- Evans, B.J. 1999. Developments in Pore Pressure Prediction Using Seismic-While-Drilling Technology. Australian Petroleum Production and Exploration Assn. J. 39 (Part 1): 461–474.
- Walls, J., Dvorkin, J., Mavko, G. et al. 2000. Use of Compressional and Shear Wave Velocity for Overpressure Detection. Presented at the Offshore Technology Conference, Houston, Texas, 1-4 May. OTC-11912-MS. http://dx.doi.org/10.4043/11912-MS.
- Badri, M.A., Sayers, C., Hussein, R.A. et al. 2001. Pore Pressure Prediction Data Using Seismic Velocities and Log Data in the Offshore Nile Delta, Egypt. Presented at the SPE Middle East Oil Show, Bahrain, 17-20 March 2001. SPE-68195-MS. http://dx.doi.org/10.2118/68195-MS.
- Carcione, J.M. and Tinivella, U. 2001. The Seismic Response to Overpressure—A Modelling Study Based on Laboratory, Well and Seismic Data. Geophysical Prospecting 49 (5): 523–539.
- Huffman, A.R. 2001. The Future of Pore-Pressure Prediction Using Geophysical Methods. Presented at the Offshore Technology Conference, Houston, Texas, 30 April-3 May. OTC-13041-MS. http://dx.doi.org/10.4043/13041-MS.
- Sayers, C.M. and Woodward, M.J. 2002. Seismic Pore-Pressure Prediction Using Reflection Tomography and 4-C Seismic Data. The Leading Edge 21 (2): 188–192. http://dx.doi.org/10.1190/1.1452611.
- Citta, F. et al. 2004. Deepwater Hazard Avoidance in a Large Top-Hole Section Using LWD Acoustic Data. The Leading Edge 23 (6): 566–573. http://dx.doi.org/10.1190/1.1766236.
- Althoff, G., Cornish, B., Varsamis, G. et al. 2004. New Concepts for Seismic Surveys While Drilling. Presented at the SPE Annual Technical Conference and Exhibition, Houston, Texas, 26-29 September 2004. SPE-90751-MS. http://dx.doi.org/10.2118/90751-MS.
- Mancilla-Castillo, J., Mendez-Hernández, E., and Santana-Fernández, J. 2003. Geopressure Evaluation from seismic data and its application for exploratory wells in Mexico. Presented at the Offshore Technology Conference, Houston, Texas, 5-8 May. OTC-15250-MS. http://dx.doi.org/10.4043/15250-MS.
- Storvoll, V., Bjorlykke, K., and Mondol, N.M. 2005. Velocity-Depth Trends in Mesozoic and Cenozoic Sediments from the Norwegian Shelf. AAPG Bulletin 89 (3): 359–381.
- Hottmann, C.E. and Johnson, R.K. 1965. Estimation of Formation Pressures from Log-Derived Shale Properties. J Pet Technol 17 (6): 717-722. SPE-1110-PA. http://dx.doi.org/10.2118/1110-PA.
- Ham, H.H. 1966. A Method of Estimating Formation Pressures from Gulf Coast Well Logs. Trans., Gulf Coast Assn. of Geological Soc. 16: 185–197.
- Japsen, P. 1999. Overpressured Cenozoic Shale Mapped from Velocity Anomalies Relative to a Baseline for Marine Shale, North Sea. Petroleum Geoscience 5: 321–336.
- Draou, A. and Osisanya, S.O. 2000. New Methods for Estimating of Formation Pressures and Fracture Gradients from Well Logs. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 1-4 October 2000. SPE-63263-MS. http://dx.doi.org/10.2118/63263-MS.
- Bowers, G.L. 1995. Pore Pressure Estimation From Velocity Data: Accounting for Overpressure Mechanisms Besides Undercompaction. SPE Drill & Compl 10 (2): 89–95. SPE-27488-PA. http://dx.doi.org/10.2118/27488-PA.
- Swarbrick, R.E. 2001. Pore-Pressure Prediction: Pitfalls in Using Porosity. Presented at the Offshore Technology Conference, Houston, Texas, 30 April-3 May. OTC-13045-MS. http://dx.doi.org/10.4043/13045-MS.
- Eaton, B.A. 1975. The Equation for Geopressure Prediction from Well Logs. Presented at the Fall Meeting of the Society of Petroleum Engineers of AIME, Dallas, 28 September–1 October. SPE 5544. http://dx.doi.org/10.2118/5544-MS.
- Eaton, B.A. and Eaton, T.L. 1997. Fracture Gradient Prediction for the New Generation. World Oil 218 (10): 93–94 and 96–100.
- Lahann, R.W., Conoco, McCarty, D.K. et al. 2001. Influence of Clay Diagenesis on Shale Velocities and Fluid-Pressure. Presented at the Offshore Technology Conference, Houston, Texas, 30 April-3 May. OTC-13046-MS. http://dx.doi.org/10.4043/13046-MS.
- Domnesteanu, P., McCann, C., and Sothcott, J. 2002. Velocity Anisotropy and Attenuation of Shale in Under- and Overpressured Conditions. Geophysical Prospecting 50 (5): 487–503.
- Dewhurst, D.N., Siggins, A.F., and Raven, M.D. 2002. Influence of Pore Pressure, Composition and Microstructure on the Acoustic Properties of Shales. Presented at the SPE/ISRM Rock Mechanics Conference, Irving, Texas, 20-23 October 2002. SPE-78197-MS. http://dx.doi.org/10.2118/78197-MS.
- Morris, R.L., Grine, D.R., and Arkfeld, T.E. 1964. Using Compressional and Shear Acoustic Amplitudes for The Location of Fractures. J Pet Technol 16 (6): 623-632. SPE-723-PA. http://dx.doi.org/10.2118/723-PA.
- Minne, J.-C. and Gartner, J. 1979. Fracture Detection in the Middle East. Presented at the Middle East Technical Conference and Exhibition, Bahrain, 25-28 February 1979. SPE-7773-MS. http://dx.doi.org/10.2118/7773-MS.
- Cheung, P.S.V. 1984. Fracture Detection Using the Sonic Tool, paper 42. Trans., 1984 Intl. Formation Evaluation Symposium, SPWLA, Paris Chapter (SAID) 1–8.
- Paillet, F.L. 1991. Qualitative and Quantitative Interpretations of Fracture Permeability by Means of Acoustic Full-Waveform Logs. The Log Analyst 32 (3): 256–270.
- Aguilera, R. 1995. Formation Evaluation by Well Log Analysis. In Naturally Fractured Reservoirs, second edition, Ch. 3, 181-315. Tulsa, Oklahoma: PennWell Publishing Co.
- Magara, K. 1976. Thickness of Removed Sedimentary Rocks, Paleopore Pressure, and Paleotemperature, Southwestern Part of Western Canada Basin. AAPG Bulletin 60 (4): 554–565.
- Magara, K. 1986. Thickness of Erosion. In Geological Models of Petroleum Entrapment, Ch. 7, 129-151. London: Elsevier Applied Science Publishers.
- Bulat, J., and Stoker, J.S. 1987. Uplift Determination from Interval Velocity Studies, UK Southern North Sea. In Petroleum Geology of North West Europe, Brooks, J. and Glennie, K. eds., Vol. 1, 293-305. London: Graham and Trottman.
- Hillis, R.R. 1993. Quantifying Erosion in Sedimentary Basins from Sonic Velocities in Shales and Sandstones. Exploration Geophysics 24: 561–566.
- Hillis, R.R., Thomson, K., and Underhill, J.R. 1994. Quantification of Tertiary Erosion in the Inner Moray Firth Using Velocity Data from the Chalk and the Kimmeridge Clay. Marine and Petroleum Geology 11 (3): 283–293.
- Hansen, S. 1996. A Compaction Trend for Cretaceous and Tertiary Shales on the Norwegian Shelf Based on Sonic-Transit Times. Petroleum Geoscience 2 (1): 59–68.
- Heasler, H.P. and Kharitonova, N.A. 1996. Analysis of Sonic Well Logs Applied to Erosion Estimates in the Bighorn Basin, Wyoming. AAPG Bulletin 80 (5): 630–646.
- Evans, D.J. 1997. Estimates of the Eroded Overburden and the Permian-Quaternary Subsidence History of the Area West of Orkney. Scottish J. of Geology 33 (Part 2): 169–181.
- Japsen, P. 1998. Regional Velocity-Depth Anomalies, North Sea Chalk—A Record of Overpressure and Neogene Uplift and Erosion. AAPG Bulletin 82 (11): 2,031–2,074.
- Densley, M.R., Hillis, R.R., and Redfearn, J.E.P. 2000. Quantification of Uplift in the Carnarvon Basin Based on Interval Velocities. Australian J. of Earth Sciences 47 (1): 111–122.
- Fuh, S.-C. 2000. Magnitude of Cenozoic Erosion from Mean Sonic Transit Time, Offshore, Taiwan. Marine and Petroleum Geology 17 (9): 1,011–1,028.
- Al-Chalabi, M. 2001. The Use of Instantaneous Velocity in Uplift Investigation. Geophysical Prospecting 49: 645–655.
- Ware, P.D. and Turner, J.P. 2002. Sonic Velocity Analysis of the Tertiary Denudation of the Irish Sea Basin. In Exhumation of The North Atlantic Margin; Timing, Mechanisms and Implications for Petroleum Exploration, Dore, A.G., Cartwright, J.A., Stoker, M.S., Turner, J.P., and White, N.J. eds., 355-370. London: Special Publication 196, Geological Society.
- Mavromatidis, A., and Hillis, R. 2005. Quantification of Exhumation in the Eromanga Basin and Its Implications for Hydrocarbon Exploration. Petroleum Geoscience 11 (1): 79–92.
- Stocks, A.E. and Lawrence, S.R. 1990. Identification of Source Rocks from Wireline Logs. In Geological Applications of Wireline Logs, A. Hurst, M.A. Lovell, and A.C. Morton eds., 241–252. Geological Soc. of London Special Publication No. 48.
- Herron, S.L. 1991. In Situ Evaluation of Potential Source Rocks by Wireline Logs. In Source and Migration Processes and Evaluation Techniques, Treatise of Petroleum Geology--Handbook of Petroleum Geology, R.K. Merrill ed. Ch. 13, 127-134. Tulsa, Oklahoma: AAPG.
- Lang, W.H. 1994. The Determination of Thermal Maturity in Potential Source Rocks Using Interval Transit Time/Interval Velocity. The Log Analyst 35 (6): 47–59.
- Meyer, B.L., and Nederlof, M.H. 1984. Identification of Source Rocks on Wireline Logs by Density/Resistivity and Sonic Transit Time/Resistivity Crossplots. AAPG Bulletin 68 (2): 121–129.
- Passey, Q.R. et al. 1990. A Practical Model for Organic Richness from Porosity and Resistivity Logs. AAPG Bulletin 74 (12): 1,777–1,794.
- Huang, Z. et al. 1996. Cyclicity in the Egret Member (Kimmeridgian) Oil Source Rock, Jeanne d’Arc Basin, Offshore Eastern Canada. Marine and Petroleum Geology 13 (1): 91–105.
- Fjaer, E. et al. 1992. Petroleum Related Rock Mechanics, 1-338. Amsterdam: Developments in Petroleum Science No. 33, Elsevier, Amsterdam.
- Montmayeur, H. and Graves, R.M. 1985. Prediction of Static Elastic/Mechanical Properties of Consolidated and Unconsolidated Sands From Acoustic Measurements: Basic Measurements. Presented at the SPE Annual Technical Conference and Exhibition, Las Vegas, Nevada, 22-26 September 1985. SPE-14159-MS. http://dx.doi.org/10.2118/14159-MS.
- Montmayeur, H. and Graves, R.M. 1986. Prediction of Static Elastic/Mechanical Properties of Consolidated and Unconsolidated Sands From Acoustic Measurements: Correlations. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 5-8 October 1986. SPE-15644-MS. http://dx.doi.org/10.2118/15644-MS.
- Holt, R.M., Ingsoy, P., and Mikkelson, M. 1989. Rock Mechanical Analysis of North Sea Reservoir Formations. SPE Form Eval 4 (1): 33-37. SPE-16796-PA. http://dx.doi.org/10.2118/16796-PA.
- Gatens III, J.M., Harrison III, C.W., Lancaster, D.E. et al. 1990. In-Situ Stress Tests and Acoustic Logs Determine Mechanical Propertries and Stress Profiles in the Devonian Shales. SPE Form Eval 5 (3): 248-254. SPE-18523-PA. http://dx.doi.org/10.2118/18523-PA.
- Yale, D.P. 1994. Static and Dynamic Rock Mechanical Properties in the Hugoton and Panoma Fields, Kansas. Presented at the SPE Mid-Continent Gas Symposium, Amarillo, Texas, 22-24 May 1994. SPE-27939-MS. http://dx.doi.org/10.2118/27939-MS.
- Kowalski, J.J. 1978. Formation Strength parameters from Well Logs, paper N. Trans., 1978 Annual Logging Symposium, SPWLA, 1–19.
- Sethi, D.K. 1981. Well Log Applications in Rock Mechanics. Presented at the SPE/DOE Low Permeability Gas Reservoirs Symposium, Denver, Colorado, 27-29 May 1981. SPE-9833-MS. http://dx.doi.org/10.2118/9833-MS.
- Tixier, M.P., Loveless, G.W., and Anderson, R.A. 1975. Estimation of Formation Strength From the Mechanical-Properties Log(incudes associated paper 6400 ). J Pet Technol 27 (3): 283-293. SPE-4532-PA. http://dx.doi.org/10.2118/4532-PA.
- Stein, N. 1976. Mechanical Properties of Friable Sands From Conventional Log Data (includes associated papers 6426 and 6427 ). J Pet Technol 28 (7): 757-763. SPE-5500-PA. http://dx.doi.org/10.2118/5500-PA.
- Onyia, E. 1988. Relationships Between Formation Strength, Drilling Strength, and Electric Log Properties. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 2–5 October. SPE-18166-MS. http://dx.doi.org/10.2118/18166-MS.
- Stein, N. 1992. Sonic Log Data Help Determine Formation Strength. Oil & Gas J. (28 December): 96.
- Raaen, A.M., Hovem, K.A., Joranson, H. et al. 1996. FORMEL: A Step Forward in Strength Logging. Presented at the SPE Annual Technical Conference and Exhibition, Denver, Colorado, 6-9 October 1996. SPE-36533-MS. http://dx.doi.org/10.2118/36533-MS.
- Anderson, T. and Walker, T. 1972. Log Derived Rock Properties for Use in Well Stimulation Design. Presented at the Fall Meeting of the Society of Petroleum Engineers of AIME, San Antonio, Texas, 8-11 October 1972. SPE-4095-MS. http://dx.doi.org/10.2118/4095-MS.
- Anderson, R.A., Ingram, D.S., and Zanier, A.M. 1973. Determining Fracture Pressure Gradients From Well Logs. J Pet Technol 25 (11): 1259-1268. SPE-4135-PA. http://dx.doi.org/10.2118/4135-PA.
- Newberry, B.M., Nelson, R.F., and Ahmed, U. 1985. Prediction of Vertical Fracture Migration Using Compression and Shear Wave Slowness. Presented at the SPE/DOE Low Permeability Gas Reservoirs Symposium, Denver, Colorado, 19–22 May. SPE-13895-MS. http://dx.doi.org/10.2118/13895-MS.
- Stein, N. 1988. How to Calculate Fracture Pressures from Well Logs. Petroleum Engineer Intl. 60 (8): 36–38.
- Bruce, S. 1990. A Mechanical Stability Log. Presented at the SPE/IADC Drilling Conference, Houston, Texas, 27 February-2 March 1990. SPE-19942-MS. http://dx.doi.org/10.2118/19942-MS.
- Stein, N. and Hilchie, D.W. 1972. Estimating the Maximum Production Rate Possible from Friable Sandstones Without Using Sand Control. J Pet Technol 24 (9): 1157-1160. SPE-3499-PA. http://dx.doi.org/10.2118/3499-PA.
- Edwards, D.P., Sharma, Y., and Charron, A. 1983. Zones of Sand Production Identified by Log-Derived Mechanical Properties—A Case Study, paper S. Trans., 1983 European Formation Evaluation Symposium, SPWLA, London Chapter, 1–23.
- Ong, S., May, A., George, I. et al. 2000. An Accurate Characterization of Sand Strength in Weak and Unconsolidated Formations Aids Offshore Production Test Designs - A Bohai Bay Case Study. Presented at the International Oil and Gas Conference and Exhibition in China, Beijing, China, 7-10 November 2000. SPE-64738-MS. http://dx.doi.org/10.2118/64738-MS.
- Fertl, W.H. and DeVries, M.R. 1997. Coal Evaluation Using Geophysical Well Logs, paper F. Trans., 1997 Formation Evaluation Symposium, Canadian Well Logging Society, 1–17.
- Bond, L.O., Alger, R.P., and Schmidt, A.W. 1971. Well Log Applications in Coal Mining and Rock Mechanics. Trans., AIME 250: 355–362.
- Kowalski, J. and Fertl, W.H. 1977. Application of Geophysical Well Logging to Coal Mining Operations. Energy Sources 3 (2): 133–147.
- Hornby, B. 1989. Imaging of near‐borehole structure using full‐waveform sonic data. Geophysics 54 (6): 747–757. http://dx.doi.org/10.1190/1.1442702.
- Esmeroy, C. et al. 1997. Sonic Imaging—A Tool for High-Resolution Reservoir Description, paper BH 2.7. Expanded Abstracts, 1997 Annual Meeting Technical Program, SEG, 1, 278–281.
- Chabot, L. et al. 2001. Single-Well Imaging Using the Full Waveform of an Acoustic Sonic, paper BH 3.7. Expanded Abstracts, 2001 Annual Meeting Technical Program, SEG, 420–423.
- Yamamoto, H., Watanabe, S., Koelman, J.M.V. et al. 2000. Borehole Acoustic Reflection Survey Experiments in Horizontal Wells for Accurate Well Positioning. Presented at the SPE/CIM International Conference on Horizontal Well Technology, Calgary, Alberta, Canada, 6-8 November 2000. SPE-65538-MS. http://dx.doi.org/10.2118/65538-MS.
- Tang, X.M. 2004. Imaging Near-Borehole Structure Using Directional Acoustic-Wave Measurement. Geophysics 69 (6): 1378–1386. http://dx.doi.org/10.1190/1.1836812.
- Zheng, Y., and Tang, X. 2005. Imaging Near-Borehole Structure Using Acoustic Logging Data with Prestack F-K Migration, paper BG 2. Expanded Abstracts, 2005 Annual Meeting Technical Program, SEG, 360–363.
- Mari, J.-L., Glangeaud, F., and Coppens, F. 1999. Signal Processing for Geologists and Geophysicists, 1-458. Paris: Editions Technip.
- Neidell, N.S. and Taner, M.T. 1971. Semblance and Other Coherency Measures for Multichannel Data. Geophysics 36 (3): 482–497. http://dx.doi.org/10.1190/1.1440186.
- Kimball, C.B. and Marzetta, T.M. 1984. Semblance Processing of Borehole Acoustic Array Data. Geophysics 49 (3): 274–281. http://dx.doi.org/10.1190/1.1441659.
- Kimball, C.V. 1998. Shear Slowness Measurements by Dispersive Processing of the Borehole Flexural Mode. Geophysics 63 (2): 337–344. http://dx.doi.org/10.1190/1.1444333.
- McFadden, P.L., Drummond, B.J., and Kravis, S. 1986. The Nth-Root Stack—Theory, Applications, and Examples. Geophysics 51 (10): 1,879–1,892. http://dx.doi.org/10.1190/1.1442045.
- Kravis, S. 1980. The Nth Root Slant Stack—A New Method of Coherency Enhancement. First Break 8 (9): 339–344.
- Brie, A. et al. 1997. Shear Slowness Determination from Dipole Measurements, paper F. Trans., 1997 Annual Logging Symposium, SPWLA, 1–14.
- Wang, Z. 2002. Seismic Anisotropy in Sedimentary Rocks, Part 1—A Single-Plug Laboratory Method; Part 2—Laboratory Data. Geophysics 67 (5): 1415–1440. http://dx.doi.org/10.1190/1.1512787.
- Crampin, S. and Chastin, S. 2003. A Review of Shear Wave Splitting in the Crack-Critical Crust. Geophysical J. Intl. 155: 221–240.
- Helbig, K. and Thomsen, L. 2005. 75-plus Years of Anisotropy in Exploration and Reservoir Seismics: A Historical Review of Concepts and Methods. Geophysics 70 (6): 9ND–23ND. http://dx.doi.org/10.1190/1.2122407.
- Franco, J.L.A., de la Torre, H.G., Ortiz, M.A.M. et al. 2005. Using Shear-Wave Anisotropy To Optimize Reservoir Drainage And Improve Production In Low-Permeability Formations in the North of Mexico. Presented at the SPE Annual Conference and Technical Exhibition, Dallas, Texas, 9–12 October. SPE-96808-MS. http://dx.doi.org/10.2118/96808-MS.
- Klimentos, T. 2003. Shear-Wave Anisotropy Applications for Perforation Strategy and Production Optimization in Oil Bearing Porous Sands, paper LL. Trans., 2003 Annual Logging Symposium, SPWLA, 1–12.
- Klimentos, T., Farghaly, A., and Qleibo, M. 2003. Finding Faults with Shear-Wave Anisotropy, paper F. Trans., 2003 Annual Logging Symposium, SPWLA, 1–12.
- Furre, A.-K. and Brevik, I. 1998. Characterization of Angle Dependency in Sonic Logs, paper BH 2.1. Expanded Abstracts, 1998 Annual Meeting Technical Program, SEG, 292–295.
- Hornby, B., Howie, J., and Ince, D. 1999. Anisotropy Correction for Deviated Well Sonic Logs—Application to Seismic Well Tie, paper BH/RP 4.7. Expanded Abstracts, 1999 Annual Meeting Technical Program, SEG, 112–115.
- Wang, Z. 2001. Seismic Anisotropy in Sedimentary Rocks, paper RP 2.4. Expanded Abstracts, 2001 Annual Meeting Technical Program, SEG, 1740–1743.
- Sun, Y. et al. 2003. Effects of Stress-Induced Anisotropy on Monopole and Dipole Logging, paper RBG P1.2. Expanded Abstracts, 2003 Annual Meeting Technical Program, SEG, 1282–1285.
- Koster, K. et al. 1994. Applied Production Geophysics Using Shear-Wave Anisotropy—Production Applications for the Dipole Shear Imager and the Multicomponent VSP, paper DP1.1. Expanded Abstracts, 1994 Annual Meeting Technical Program, SEG, 233–235.
- Esmersoy, C. et al. 1995. Fracture and Stress Evaluation Using Dipole-Shear Anisotropy Logs, paper J. Trans., 1995 Annual Logging Symposium, SPWLA, 1–12.
- Tang, X. 1996. Processing Dipole Waveform Logs for Formation Alteration Identification, paper BG 3.4. Expanded Abstracts, 1996 Annual Meeting Technical Program, SEG, 1, 162–165.
- Huang, X. et al. 1999. Formation Stress Estimation Using Standard Acoustic Logging, paper BH/RP 2.6. Expanded Abstracts, 1999 Annual Meeting Technical Program, SEG, 53–56.
- Tang, X., Cheng, N., and Cheng, A. 1999. Identifying and Estimating Formation Stress from Borehole Monopole and Cross-Dipole Acoustic Measurements, paper QQ. Trans., 1999 Annual Logging Symposium, SPWLA, 1–14.
- Badri, M., Sousa, S., and Klimentos, T. 2000. Shear Anisotropy Applications in Production Optimization, Western Desert, Egypt, paper RPB 1.5. Expanded Abstracts, 2000 Annual Meeting Technical Program, SEG, 1695–1698.
- Tang, X., Patterson, D., and Markovic, M. 2000. Fracture Measurements in Open/Cased Holes Using Cross-Dipole Logging—Theory and Field Results, paper RPB 1.6. Expanded Abstracts, 2000 Annual Meeting Technical Program, SEG, 1699–1702.
- Tang, X.M. 2003. Determining Shear-Wave Transverse Isotropy from Borehole Stoneley Waves. Geophysics 68 (1): 118–126. http://dx.doi.org/10.1190/1.1543199.
- Cheng, N. and Cheng, C.H. 1996. Estimations of Formation Velocity, Permeability, and Shear-Wave Velocity Using Acoustic Logs. Geophysics 61 (2): 437–443. http://dx.doi.org/10.1190/1.1443971.
- Murray, D., Plona, T., and Valero, H.P. 2004. Case Study of Borehole Sonic Dispersion Curve Analysis, paper BB. Trans., 2004 Annual Logging Symposium, SPWLA, 1–14.
- Sinha, B.K., Bratton, T.R., Cryer, J.V. et al. 2005. Near-Wellbore Alteration and Formation Stress Parameters Using Borehole Sonic Data. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 9-12 October 2005. SPE-95841-MS. http://dx.doi.org/10.2118/95841-MS.
- Plona, T.J., Winkler, K.W., Sinha, B.K. et al. 1998. Measurement of Stress Direction and Mechanical Damage Around Stressed Boreholes Using Dipole and Microsonic Techniques. Presented at the SPE/ISRM Rock Mechanics in Petroleum Engineering, Trondheim, Norway, 8-10 July 1998. SPE-47234-MS. http://dx.doi.org/10.2118/47234-MS.
- Plona, T. et al. 1999. Stress-Induced Dipole Anisotropy—Theory, Experiment and Field Data, paper RR. Trans., presented at the 1999 Annual Logging Symposium, SPWLA, 1–14.
- Plona, T.J. et al. 2000. Using Acoustic Anisotropy, paper H. Trans., 2000 Annual Logging Symposium, SPWLA, 1–12.
- Sinha, B.K., Kane, M.R., and Frignet, B. 2000. Dipole Dispersion Crossover and Sonic Logs in Limestone Reservoir. Geophysics 65 (2): 390–407. http://dx.doi.org/10.1190/1.1444734.
- Xu, S., Wu, X., Huang, X. et al. 2005. Evaluation of Anisotropic Rock Properties in Sedimentary Rocks From Well Logs. Presented at the Offshore Technology Conference, Houston, Texas, 2-5 May. OTC-17251-MS. http://dx.doi.org/10.4043/17251-MS.
- Crampin, S. 1985. Evaluation of Anisotropy by Shear-Wave Splitting. Geophysics 50 (1): 142–152. http://dx.doi.org/10.1190/1.1441824.
- Rai, C.S. and Hanson, K.E. 1988. Shear-Wave Velocity Anisotropy in Sedimentary Rocks—A Laboratory Study. Geophysics 53 (6): 800–806. http://dx.doi.org/10.1190/1.1442515.
- Yale, D.P. and Sprunt, E.S. 1989. Prediction of Fracture Direction Using Shear Acoustic Anisotropy. The Log Analyst 30 (2): 65–70.
- Winkler, K.W. 1996. Azimuthal Velocity Variations Caused by Borehole Stress Concentrations. J. of Geophysical Research 101 (B4): 8,615–8,621.
- Esmersoy, C. et al. 1994. Dipole Shear Anisotropy Logging, paper SL 3.7. Expanded Abstracts, 1994 Annual Meeting Technical Program, SEG, 1,139–1,142.
- Hatchell, P.J. et al. 1995. Quantitative Comparison Between a Dipole Log and VSP in Anisotropic Rocks from Cymric Oil Field, California, paper BG1.4. Expanded Abstracts, 1995 Annual Meeting Technical Program, SEG, 13–16.
- Beckham, W.E. 1996. Seismic Anisotropy and Natural Fractures from VSP and Borehole Sonic Tools—A Field Study. Geophysics 61 (2): 456–466. http://dx.doi.org/10.1190/1.1443973.
- Walls, J.D., Dvorkin, J., and Mavko, G. 1996. In-Situ Stress Magnitude and Azimuth from Well Log Measurements, Final Report, Report GRI-95/0356, 1-147. Chicago, Illinois: Gas Research Inst.
- Grandi, S.K. et al. 2003. In Situ Stress Modeling at a Borehole—A Case Study, paper BH 1.4. Expanded Abstracts, 2003 Annual Meeting Technical Program, SEG, 297–300.
- De, G.S. and Schmitt, D.P. 2005. Issues With Shear-wave Azimuthal Anisotropy in Highly Deviated Wells. Presented at the Offshore Technology Conference, Houston, Texas, 2-5 May. OTC-17647-MS. http://dx.doi.org/10.4043/17647-MS.
- Tang, X.M. and Patterson, D. 2005. Characterizing Seismic Anisotropy Using Cross-Dipole Measurement in Deviated Wells, paper BG 2.5. Expanded Abstracts, 2005 Annual Meeting Technical Program, SEG, 372–375.
- Patterson, D.J. and Tang, X. 2005. Pit Falls In Dipole Logging - Anisotropy: Cause Of Discrepancy In Borehole Acoustic Measurements. Presented at the Offshore Technology Conference, Houston, Texas, 2-5 May. OTC-17644-MS. http://dx.doi.org/10.4043/17644-MS.
- Hornby, B.E., Howie, J.M., and Ince, D.W. 2003. Anisotropy Correction for Deviated-Well Sonic Logs—Application to Seismic Well Tie. Geophysics 68 (2): 464–471. http://dx.doi.org/10.1190/1.1567212
- Alford, R.M. 1986. Shear Data in the Presence of Azimuthal Anisotropy—Dilley, Texas, paper S9.6. Expanded Abstracts, 1986 Annual Technical Meeting Program, SEG.
- Tang, X.M. and Chunduru, R.K. 1999. Simultaneous Inversion of Formation Shear-Wave Anisotropy parameters from Cross-Dipole Acoustic Array Waveform Data. Geophysics 64 (5): 1,502–1,511. http://dx.doi.org/10.1190/1.1444654.
- Gelinsky, S. et al. 1998. Anisotropic Permeability in Fractured Reservoirs from Integrated Acoustic Measurements, paper RC 3.5. Expanded Abstracts, 1998 Annual Meeting Technical Program, SEG, 956–959.
- Smith, M.B. 1979. Effect of Fracture Azimuth on Production With Application to the Wattenburg Gas Field. Presented at the SPE Annual Technical Conference and Exhibition, Las Vegas, Nevada, 23-26 September 1979. SPE-8298-MS. http://dx.doi.org/10.2118/8298-MS.
- Hubbert, M.K. and Willis, D.G. 1957. Mechanics of Hydraulic Fracturing. Trans., AIME 210: 153–166.
- Cipolla, C.L., Liu, D., and Kyte, D.G. 1994. Practical Application of In-Situ Stress Profiles. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 25-28 September 1994. SPE-28607-MS. http://dx.doi.org/10.2118/28607-MS.
- Cipolla, C.L. 1996. Hydraulic Fracture Technology in the Ozona Canyon and Penn Sands. Presented at the Permian Basin Oil and Gas Recovery Conference, Midland, Texas, 27-29 March 1996. SPE-35196-MS. http://dx.doi.org/10.2118/35196-MS.
- De, G.S., Winterstein, D.F., Johnson, S.J. et al. 1998. Predicting Natural or Induced Fracture Azimuths From Shear-Wave Anisotropy. SPE Res Eval & Eng 1 (4): 311-318. SPE-50993-PA. http://dx.doi.org/10.2118/50993-PA.
- Garg, A., Desai, A.M., and Towler, B.F. 1997. Horizontal Stresses in Anisotropic Formation and Prediction of Hydraulic Fracture Direction. Presented at the SPE Western Regional Meeting, Long Beach, California, 25-27 June 1997. SPE-38342-MS. http://dx.doi.org/10.2118/38342-MS.
- Aquila, F.J., Barajas, J.S., Mesa, H. et al. 2003. Using Cross Dipole Sonic Anistropy Data to Improve Reservoir Understanding in the Southern/Marine Areas of Mexico. Presented at the SPE Annual Technical Conference and Exhibition, Denver, Colorado, 5-8 October 2003. SPE-84204-MS. http://dx.doi.org/10.2118/84204-MS.
- Barajas, J.S., Patino, A.H., Garcia, E.R. et al. 2004. Case History - Cased Hole Dipole Sonic Applications in Mexico. Presented at the SPE Annual Technical Conference and Exhibition, Houston, Texas, 26-29 September 2004. SPE-90703-MS. http://dx.doi.org/10.2118/90703-MS.
- Klimentos, T. and McCann, C. 1990. Relationships Between Compressional Wave Attenuation, Porosity, Clay Content, and Permeability in Sandstones. Geophysics 55(8): 998–1,014. http://dx.doi.org/10.1190/1.1442928.
- Klimentos, T. 1995. Attenuation of P- and S-Waves as a Method of Distinguishing Gas from Oil and Water. Geophysics 60 (2): 447–458. http://dx.doi.org/10.1190/1.1443782.
- Akbar, N., Dvorkin, J., and Nur, A. 1993. Relating P-Wave Attenuation to Permeability. Geophysics 58 (1): 20–29. http://dx.doi.org/10.1190/1.1443348.
- Sun, X. et al. 2000. P- and S-Wave Attenuation Logs from Monopole Sonic Data. Geophysics 65 (3): 755–765. http://dx.doi.org/10.1190/1.1444774.
- Burns, D.R. 1991. Predicting Relative and Absolute Variations of In-Situ Permeability from Full-Waveform Acoustic Logs. The Log Analyst 32 (3): 246–255.
- Tang, X.M., Cheng, C.H., and Toksoz, M.N. 1991. Dynamic Permeability and Borehole Stoneley Waves—A Simplified Biot-Rosenbaum Model. J. of the Acoustical Soc. of America 90 (3): 1,632–1646.
- Tang, X., Cheng, C.H., and Toksoz, M.N. 1991. Stoneley-Wave Propagation in a Fluid Filled Borehole with a Vertical Fracture. Geophysics 56(4): 447–460. http://dx.doi.org/10.1190/1.1443062.
- Saxena, V. 1994. Permeability Quantification from Borehole Stoneley Waves, paper T. Trans., 1994 Annual Logging Symposium, SPWLA, 1–22.
- Kimball, C.V. and Endo, T. 1998. Quantitative Stoneley Mobility Inversion, paper BH 1.1. Expanded Abstracts, 1998 Annual Meeting Technical Program, SEG, 252–255.
- Tang, X.M. et al. 1997. Permeability from Borehole Acoustic Logs—An Overview with Recent Advances, paper BH 2.6. Expanded Abstracts, 1997 Annual Meeting Technical Program, SEG, 274–277.
- Sinha, A. et al. 1998. A New Method for Deriving Permeability from Borehole Stoneley Waves and Its Application in the North Mongas Field of Eastern Venezuela, paper P. Trans., 1998 Annual Logging Symposium, SPWLA, 1–12.
- Canady, W., Spooner, P., and Vasquez, R.B. 2005. Permeability Estimation From Stoneley Amplitude, Corrected for Borehole Geometry and Rugosity. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 9-12 October 2005. SPE-96598-MS. http://dx.doi.org/10.2118/96598-MS.
- Geerits, T. et al. 1999. Comparison Between Stoneley, NMR, and Core-Derived Permeabilities, paper TT. Trans., 1999 Annual Logging Symposium, SPWLA, 1–14.
- Qobi, L., Kuijper, A.D., Tang, X.M. et al. 2000. Permeability Determination from Stoneley Waves in the Ara Group. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 1-4 October 2000. SPE-63140-MS. http://dx.doi.org/10.2118/63140-MS.
- Tang, X.M., and Patterson, D. 2004. Estimating Formation Permeability and Anisotropy from Borehole Stoneley Waves, paper W. Trans ., 2004 Annual Logging Symposium, SPWLA, 1–14.
- Tang, X.M., Altunbay, M., and Shorey, D. 1998. Joint Interpretation of Formation Permeability from Wireline Acoustic NMR, and Image Log Data, paper KK. Trans., 1998 Annual Logging Symposium, SPWLA, 1–13.
- Aladeddin, E., Tchambaz, M., and Al-Adani, N. 2004. Combining NMR and Stoneley Analysis for a Better Estimation of Permeability in Carbonate Reservoirs, paper E. Proceedings, 2004 Formation Evaluation Symposium of Japan, Society of Petrophysicists and Well Log Analysts, Japan Chapter, 7 p.
- Tang, X.M. and Patterson, D. 2001. Detecting Thin Gas Beds in Formations Using Stoneley Wave Reflection and High-Resolution Slowness Measurements, paper OO. Trans., 2001 Annual Logging Symposium, SPWLA, 1–14.
- Hornby, B.E. et al. 1989. Fracture Evaluation Using Reflected Stoneley-Wave Arrivals. Geophysics 54 (10): 1,274–1,288. http://dx.doi.org/10.1190/1.1442587.
- Medlin, W.L. and Schmitt, D.P. 1994. Fracture Diagnostics With Tube-Wave Reflection Logs. J Pet Technol 46 (3): 239-248. SPE-22872-PA. http://dx.doi.org/10.2118/22872-PA.
SI Metric Conversion Factors
|dB/ft||×||3.048*||E − 01||=||dB/m|
|ft||×||3.048*||E − 01||=||m|
|ft/sec||×||3.048*||E − 01||=||m/sec|
|in.||×||2.54*||E + 00||=||cm|
|psi/ft||×||2.262 059||E + 01||=||kPa/m|
|μsec/ft||×||3.280 840||E + 00||=||μsec/m|
Conversion factor is exact.