Rarefaction
Rarefaction in seismic wave physics is the phase of a compressional (P-wave) or acoustic wave cycle in which the propagating disturbance causes the particles of the elastic medium to move apart from their equilibrium positions, reducing the local density and pressure below the ambient value and creating a momentary zone of lower pressure (negative pressure perturbation) that travels through the medium at the P-wave velocity of that medium; rarefaction is the opposite of compression (also called condensation), which is the phase in which particles are pushed together, increasing local density and pressure above ambient, and together these two complementary phases constitute the sinusoidal pressure variation that defines the compressional wave -- the compression phase carries positive acoustic pressure while the rarefaction phase carries negative acoustic pressure, with the amplitude of each phase equal in magnitude for a symmetric wave but potentially asymmetric for nonlinear wave propagation at high amplitude in fluids near their tensile strength limit; in seismic reflection interpretation, the polarity of the wavelet at a reflecting interface is determined by the acoustic impedance contrast (the product of velocity and density) across the boundary, with a positive reflection coefficient (impedance increases downward, such as at a shale-to-carbonate interface) generating a compression-first wavelet (positive peak first in SEG normal polarity convention) and a negative reflection coefficient (impedance decreases downward, as at a gas-sand interface where the low-impedance gas sand underlies a higher-impedance shale) generating a rarefaction-first wavelet (negative trough first in SEG normal polarity, indicating a soft reflector), making the identification of rarefaction-first arrivals a key diagnostic for fluid-related anomalies in seismic interpretation.
Key Takeaways
- The distinction between compression and rarefaction phases in a seismic wave is fundamental to seismic polarity convention and amplitude versus offset (AVO) analysis: the Society of Exploration Geophysicists (SEG) normal polarity convention defines a positive reflection coefficient (hard reflector, such as a dense carbonate below a less dense shale) as producing a positive number on the seismic trace, displayed as a peak (conventionally shaded black or red), while a negative reflection coefficient (soft reflector, such as a gas sand below a shale) produces a negative number displayed as a trough (conventionally shaded white or blue); the physical wave that arrives at the receiver after reflecting from a gas-sand top is a rarefaction-first wavelet because the downgoing compression phase reflects off the lower-impedance gas sand and returns to the surface as a rarefaction (the reflected wave's polarity is inverted relative to the downgoing wave by the negative reflection coefficient), so interpreters look for trough-first arrivals at the top of a gas sand as a direct hydrocarbon indicator (DHI), which requires knowing the acquisition and processing polarity conventions with certainty before making this interpretation.
- Rarefaction waves play a critical role in reservoir pressure dynamics during production: when a well is put on production and bottomhole flowing pressure drops below the initial reservoir pressure, a pressure rarefaction (depletion wave) propagates outward through the formation at a velocity determined by the pore compressibility and fluid compressibility (the acoustic velocity of the reservoir fluid-rock system, typically 1,200 to 2,500 m/s in liquid-saturated reservoirs and 400 to 900 m/s in gas-saturated reservoirs), with the depletion wave eventually reaching the reservoir boundary or adjacent wells; in tight reservoirs with very low permeability (less than 0.1 millidarcy), the diffusivity coefficient (k/phi*mu*ct) is so low that the pressure rarefaction propagates outward only tens to hundreds of meters in the first year of production, explaining why tight gas and tight oil wells exhibit rapid early decline as the near-wellbore zone depletes before the pressure wave has propagated to a significant drainage radius.
- Acoustic rarefaction in liquids (cavitation) occurs when the negative pressure excursion of a rarefaction wave exceeds the tensile strength of the liquid (approximately 0.1 to 1.0 MPa for degassed water and much lower for saturated water or crude oil with dissolved gas), causing the liquid to rupture and form vapor-filled cavities that subsequently collapse violently during the following compression phase; cavitation in pump impellers (caused by high-velocity fluid flow creating local rarefaction zones near the impeller blades) causes physical damage to metal surfaces through the high-pressure water hammer generated by bubble collapse (estimated at 300 to 1,000 MPa locally), reducing pump efficiency and causing erosive pitting that limits pump life; in oilfield operations, cavitation is a concern in high-rate centrifugal pumps used for produced water injection and in multistage hydraulic fracturing pumps where the suction-side pressure can fall below the vapor pressure of the fluid being pumped, particularly when pumping gas-saturated fluid or when suction pressure drops due to restricted intake or high fluid viscosity.
- Seismic rarefaction identification in amplitude interpretation requires careful attention to processing history and wavelet phase: a zero-phase processed seismic dataset (the industry standard for interpretation) contains both a compression phase and a rarefaction phase at every reflecting interface, with the peak of the zero-phase wavelet centered on the reflector; a minimum-phase dataset (the output of deconvolution without phase rotation) has the energy front-loaded so the first arriving energy at a hard reflector is a compression and at a soft reflector is a rarefaction, but the wavelet shape differs from zero-phase; bright-spot DHI identification relies on recognizing that a high-amplitude trough (rarefaction) at the top of a gas sand is physically consistent with a negative reflection coefficient, while a high-amplitude peak (compression) at the same location would represent a polarity reversal inconsistent with gas and requiring a different explanation such as carbonate cementation or tight sandstone; AVO analysis extends this polarity discrimination by examining how the reflection amplitude changes with offset (incidence angle), with gas sands typically showing amplitude that increases with offset for both the compression phase at the base of the sand (positive reflection) and the rarefaction phase at the top (negative reflection), while brine sands show more complex offset-dependent behavior.
- Rarefaction waves in gas-dominated reservoirs and pipelines are the physical mechanism behind the blowdown pressure pulse that propagates through a gas system when a valve is rapidly opened or a pipeline segment is vented: when a high-pressure gas column is suddenly depressurized at one end (by opening a blowdown valve), a rarefaction wave travels back through the gas at the gas acoustic velocity (approximately 350 to 500 m/s for natural gas at typical pipeline conditions), reducing the pressure sequentially along the pipeline from the open end inward; the rarefaction wave velocity in gas is approximately equal to the gas sound speed modified by the Mach number of the flowing gas, and in retrograde gas condensate systems, the rarefaction can locally reduce pressure below the dew point, causing condensate dropout along the pipeline wall that was not present at initial static conditions, a phenomenon relevant to pipeline blowdown design and hydrocarbon recovery from gas condensate reservoirs during pressure depletion below the dew point.
Fast Facts
The physical concept of rarefaction as the low-pressure phase of a compressional wave was established by Robert Boyle and later formalized by Isaac Newton in his Principia Mathematica (1687), which contained the first analytical derivation of the speed of sound in air as proportional to the square root of pressure divided by density -- though Newton's calculation underestimated the actual speed because he assumed isothermal rather than adiabatic compression and rarefaction. Pierre-Simon Laplace corrected Newton's formula in 1816 by incorporating the adiabatic gamma factor, yielding the familiar formula c = sqrt(gamma*P/rho) that is still used for the speed of sound in ideal gases. In seismic exploration, the polarity conventions for displaying compression and rarefaction as peaks and troughs on seismic sections were not standardized across the industry until the SEG established the normal polarity convention in 1975, prior to which different acquisition companies and processing vendors used inconsistent polarity assignments that made cross-company seismic comparison hazardous for DHI interpretation.
What Is Rarefaction?
Rarefaction is the phase of a compressional (P-wave) seismic or acoustic wave during which the medium's particles separate and local pressure drops below ambient, creating the negative-pressure half-cycle of the wave. It is the physical complement of the compression phase. In seismic reflection interpretation, rarefaction-first arrivals (troughs in SEG normal polarity) indicate negative reflection coefficients at interfaces where impedance decreases downward, such as at the top of a gas sand below shale. Identifying rarefaction-first versus compression-first arrivals is essential for DHI analysis and AVO interpretation, requiring precise knowledge of acquisition and processing polarity conventions.
Synonyms and Related Terminology
Rarefaction is also called dilatation, dilation, or the negative pressure phase of a compressional wave. In seismic contexts it may be called a trough or soft reflection. Related terms include compressional wave (P-wave, the seismic body wave in which particle motion is parallel to the propagation direction, alternating between compression (particle convergence, positive pressure) and rarefaction (particle divergence, negative pressure) as the wave propagates through solid rock, fluid, or gas; P-waves travel faster than S-waves in all media and are the primary wave type used in conventional seismic reflection surveys for subsurface imaging), reflection coefficient (the ratio of reflected to incident seismic wave amplitude at an impedance boundary, equal to (Z2 - Z1)/(Z2 + Z1) where Z1 and Z2 are the acoustic impedances of the upper and lower media; a positive reflection coefficient produces a compression-first (peak) reflected wavelet and a negative reflection coefficient produces a rarefaction-first (trough) wavelet, the latter being the expected DHI response at the top of a gas sand), polarity (the convention defining whether a positive reflection coefficient is displayed as a peak or a trough on a seismic section; in SEG normal polarity convention, a positive impedance increase (hard reflector) generates a peak and a negative impedance decrease (soft reflector, typical of a gas sand) generates a trough representing a rarefaction; incorrect polarity identification causes DHI interpretations to be reversed, misidentifying gas as brine and brine as gas), amplitude versus offset (AVO, the analysis of how seismic reflection amplitude varies with source-receiver offset (incidence angle), used to discriminate hydrocarbon-related reflectors from lithological reflectors; gas sands typically display characteristic AVO behavior (Class II or III) in which the rarefaction amplitude at the top of the sand increases with offset, providing a hydrocarbon indicator independent of the absolute amplitude value), and direct hydrocarbon indicator (DHI, a seismic attribute that is physically predictable from rock physics for a hydrocarbon-bearing reservoir, including bright spots (high amplitude), flat spots (horizontal reflectors at fluid contacts), and polarity reversals (where a rarefaction-first top-sand reflection in a hydrocarbon case inverts to compression-first in the equivalent brine-saturated case); DHI identification requires careful polarity analysis to correctly assign rarefaction versus compression phase to each candidate anomaly).
Why Polarity and Rarefaction Identification Are Non-Negotiable in DHI Interpretation
A seismic interpreter who inverts the polarity convention misidentifies every gas sand as brine and every brine sand as a DHI anomaly. This is not a theoretical risk -- polarity inversions in processing workflows, phase rotations applied for aesthetic wavelet shape, and inconsistent vendor conventions have caused real exploration wells to be drilled on what turned out to be polarity-reversed noise anomalies while genuine gas sands were dismissed as "hard" reflectors. The rarefaction phase of a wave is not just an abstract physics concept; it is the physical signature that a negative impedance contrast produces in the earth, and recognizing it correctly requires the interpreter to track polarity from the field acquisition parameters through every processing step to the final migrated stack. Getting that chain right is the difference between a discovery well and a dry hole drilled on a misread wavelet.