Biot Theory

Biot theory is the mathematical framework developed by Maurice Biot in 1956 to describe how elastic (acoustic) waves propagate through porous, fluid-saturated rock. The key insight of Biot's work was that porous rock saturated with a compressible fluid (water, oil, or gas) supports not one but two distinct compressional (P-wave) modes: a fast P-wave in which the solid frame and pore fluid move together, and a slow P-wave in which the solid and fluid move in opposite directions (the Biot slow wave). Biot theory is the theoretical foundation for acoustic well log interpretation, seismic rock physics, and the Gassmann fluid substitution equations used to model how seismic velocities change when reservoir fluids change. It is one of the most practically important theoretical contributions in applied geophysics.

Key Takeaways

The Biot slow wave (also called the Biot wave or the Type II P-wave) travels much more slowly than the fast P-wave and is heavily attenuated. At seismic frequencies (10 to 100 Hz), the slow wave is essentially diffusive rather than propagating: it dies out within millimetres of its source. At ultrasonic frequencies (megahertz range) used in laboratory core measurements, the slow wave can be observed directly. At sonic logging frequencies (1 to 20 kHz), the slow wave contributes to attenuation of the fast P-wave through a mechanism called Biot global flow.
At low frequencies (seismic frequencies), Biot theory reduces to the Gassmann equations, which relate the saturated rock bulk modulus to the dry frame modulus, fluid modulus, and mineral modulus. The Gassmann equations are the workhorse of seismic rock physics: they allow geophysicists to calculate how the P-wave velocity of a rock changes when one fluid is replaced by another (for example, brine replaced by gas), which is the basis for 4D seismic interpretation of reservoir fluid changes.
Biot theory predicts that attenuation (loss of wave energy during propagation) depends on the permeability of the rock and the viscosity of the pore fluid. High-permeability rock with low-viscosity fluid (gas) shows significant attenuation due to fluid flow driven by the acoustic wave (squirt flow in the pore space). Low-permeability rock with high-viscosity fluid shows less Biot-type attenuation. This permeability-attenuation relationship is the basis for attempts to measure formation permeability from acoustic logs.
The dynamic elastic moduli measured by acoustic logs (compressional and shear wave velocities) differ from the static elastic moduli measured in core tests (deformation under slow applied stress). Biot theory provides part of the explanation for this discrepancy: at the fast loading rates of acoustic waves, the pore fluid has no time to redistribute, giving a stiffer (higher modulus) response than the slow loading of a static core test, where fluid can flow and equalize pressure.
Biot's original 1956 papers assumed isotropic rock with a single connected pore fluid. Real reservoir rocks are often anisotropic (different properties in different directions) due to layering, aligned fractures, or stress anisotropy. Anisotropic extensions of Biot theory are used in seismic processing of fractured reservoirs, where P-wave and S-wave velocities vary with the direction of wave propagation relative to the fracture orientation (seismic anisotropy).

What Is Biot Theory and Why Does It Matter?

Squeeze a wet sponge. Water comes out through the pores as the solid frame of the sponge compresses. Release the sponge and water flows back in. Now squeeze the sponge very quickly, faster than the water can move. The sponge feels stiffer because the incompressible water cannot escape fast enough and helps support the load. The stiffness of the sponge depends on how fast you squeeze it and how easily water can flow out through the pores.

Rock containing oil, gas, or water in its pores behaves the same way when an acoustic wave passes through it. The wave alternately compresses and expands the rock, and the pore fluid tries to flow in response. Whether the fluid has time to flow or not depends on the frequency of the wave and the permeability of the rock. Biot theory is the mathematical description of this interaction between the elastic rock frame and the compressible pore fluid across all frequencies.

In practical oilfield terms, Biot theory is important because it tells you how changing the pore fluid changes the seismic wave velocity. This is the central question in 4D seismic monitoring of producing reservoirs: if oil is replaced by water during waterflooding, how does the seismic P-wave velocity change? Biot theory (through the Gassmann equations) provides the answer, and that answer is what makes 4D seismic interpretable as a tool for tracking fluid fronts in a reservoir.

Fast Facts

Maurice Biot (1905-1985) was a Belgian-American applied mathematician and physicist who published his landmark papers on acoustic wave propagation in porous media in the Journal of the Acoustical Society of America in 1956. At the time, the papers were not widely read outside of mechanical engineering and applied mathematics. It was not until the 1970s and 1980s, when the oil and gas industry began using seismic data quantitatively for reservoir characterization rather than just for structural mapping, that Biot's framework became central to geophysics. Today Biot's 1956 papers are among the most cited works in applied geophysics and rock physics. His name is attached to the Biot slow wave, Biot global flow, the Biot-Gassmann equations, and Biot-Willis coefficients.

The Gassmann Equations: Biot Theory in Practice

At low frequencies (below the Biot critical frequency, which for typical reservoir rocks and fluids is between 1 and 100 kHz), Biot theory simplifies to the Gassmann equations. These equations give the saturated rock bulk modulus K_sat in terms of:

K_dry (the dry frame bulk modulus, measured on dry core), K_mineral (the bulk modulus of the rock mineral grains, a material property), K_fluid (the bulk modulus of the pore fluid, which depends strongly on fluid type), and the porosity (φ). The shear modulus G_sat equals G_dry because pore fluids do not resist shear deformation (fluids have no shear strength).

With K_sat from Gassmann and the density ρ_sat (computed from grain density, fluid density, and porosity), the saturated P-wave velocity is Vp = sqrt((K_sat + 4/3 G) / ρ). This calculation, from dry frame properties to saturated P-wave velocity, is called forward Gassmann substitution. The reverse (from measured saturated velocity to dry frame properties) is inverse substitution, which is used to strip out the fluid effect from a log measurement and characterize the rock frame independently of the current fluid.

In the North Sea Brent reservoir group, geophysicists use Gassmann substitution to model how the P-wave impedance of oil-saturated sands would change if the oil were replaced by water or gas. These forward models are compared to observed amplitude changes on 4D seismic surveys to map where waterflood has swept and where bypassed oil remains.

Biot Theory and Acoustic Well Log Interpretation

Acoustic (sonic) logs measure the travel time of P-waves and S-waves through the formation adjacent to the borehole. The measured velocities depend on the rock frame properties, the porosity, and the fluid. Separating these three contributions is the core problem of log-based rock physics.

Biot theory predicts that at the frequencies of sonic logs (1 to 20 kHz), the measurement is in the "Biot undrained" regime for typical reservoir permeabilities: the pore fluid does not have time to fully equilibrate during the acoustic cycle. This means the measured sonic velocity is not the same as would be measured at seismic frequencies (where full equalization occurs). The Biot frequency correction adjusts for this difference when upscaling log velocities to seismic frequencies or downscaling seismic velocities to log frequencies.

In tight gas sands (Montney, Duvernay, Falher) where permeability is in the microdarcy range, Biot frequency corrections are small because the pore fluid barely moves at any practical frequency. In high-permeability gravel pack completions or highly permeable carbonates, the Biot correction is larger and can affect the accuracy of porosity estimation from sonic logs.

Biot theory is also called Biot's poroelastic theory or Biot's theory of acoustic propagation in porous media. Related terms include Gassmann equations (the low-frequency limit of Biot theory relating saturated rock elastic moduli to dry frame moduli, fluid moduli, and porosity; the standard rock physics tool for fluid substitution calculations), Biot slow wave (the second compressional wave mode predicted by Biot theory, in which the solid frame and pore fluid move in opposite directions; heavily attenuated at seismic frequencies but measurable in laboratory ultrasonic experiments), fluid substitution (the calculation of how seismic velocities and acoustic impedance change when one pore fluid is replaced by another; performed using the Gassmann equations derived from Biot theory), rock physics (the study of the relationships between physical properties of rocks (velocity, density, electrical resistivity) and their geological and petrophysical properties (porosity, mineralogy, fluid); Biot theory is the fundamental framework of rock physics), and seismic anisotropy (the variation of seismic wave velocity with direction, caused by rock layering, aligned fractures, or stress; anisotropic extensions of Biot theory describe wave propagation in these more complex media).

How a Biot Gassmann Calculation Identified Bypassed Oil in a Norwegian North Sea Field

A reservoir engineering team at a producing Brent sandstone field on the Norwegian Continental Shelf was analyzing a 4D seismic survey acquired 8 years after the baseline. The field had been on waterflood for 6 years. The 4D difference map showed a large amplitude increase in the eastern sector of the field, which the initial interpretation attributed to waterflood sweep: water replacing oil in the reservoir sands should produce an amplitude increase in the near-stack seismic data (because water has higher acoustic impedance than oil in this reservoir).

A rock physicist on the team ran Gassmann fluid substitution calculations on the log data from the nearest well to the anomaly. The calculation showed that replacing oil with brine in the Brent sand at this location would increase the P-wave impedance by 12 percent, which matched the observed amplitude increase. So far, the interpretation seemed confirmed.

However, the rock physicist also ran the calculation for gas-replacing-oil. Gas substitution gave an amplitude decrease, not an increase, because gas reduces density more than it reduces P-wave modulus, net-decreasing impedance. The anomaly was clearly not gas. Running the calculation for different brine salinities, the team found that the high amplitude was also consistent with a 20 percent increase in water saturation rather than full displacement, which would mean the anomaly was the oil-water transition zone rather than fully swept rock. The interpretation shifted: the high-amplitude anomaly was the migrating waterflood front, with oil still present above it and potentially trapped by a subtle fault not visible on the 3D seismic.

A targeted infill well was drilled 1.8 kilometres updip from the anomaly to a structural high that the new interpretation suggested was unswept. The well encountered 22 metres of net pay at 71 percent oil saturation in the Brent sand. It produced 2,400 cubic metres of oil per day on natural flow. The Biot-Gassmann framework that turned an "already swept" interpretation into a "look updip" recommendation was the enabling calculation. The infill well added USD 180 million in net present value to the field development.