geophysics · avo
Gassmann Fluid Substitution
Ksat = Kdry + (1 − Kdry/Kmin)² / (φ/Kfl + (1−φ)/Kmin − Kdry/Kmin²)
click formula to derive ↑
Inputs
GPa
GPa
GPa
GPa
fraction
g/cc
Description
Performs Gassmann fluid substitution to predict how the saturated bulk modulus, P- and S-wave velocities, acoustic impedance, and Vp/Vs change when the pore fluid is replaced (e.g., brine→gas, oil→brine). The foundation for 4D seismic interpretation and AVO forward modeling.
Variables
| Symbol | Unit | Description |
|---|---|---|
| Ksat | GPa | Saturated Bulk Modulus |
| Vp | ft/s | P-wave Velocity |
| Vs | ft/s | S-wave Velocity |
| AI | (g/cc)·(m/s) | Acoustic Impedance |
| Vp/Vs | — | Vp/Vs Ratio |
| Kdry | GPa | Dry-frame (drained) bulk modulus of the rock skeleton with no pore fluid — typically derived by inverting Gassmann from brine-saturated log measurements. |
| Kmin | GPa | Mineral (grain) bulk modulus of the rock matrix — common values: quartz≈37 GPa, calcite≈71 GPa, dolomite≈95 GPa. |
| Kfl | GPa | Bulk modulus of the target pore fluid to substitute in — common values: brine≈2.5 GPa, oil≈0.7-1.2 GPa, gas≈0.02-0.1 GPa. |
| G | GPa | Shear modulus of the saturated rock. Gassmann assumes shear modulus is unaffected by pore fluid, so this value is held constant through fluid substitution. |
| φ | fraction | Effective porosity of the rock, from a porosity log or core measurement. |
| ρ | g/cc | Bulk density of the rock saturated with the target fluid, from the density log or computed from mineral/fluid mixing. |
Assumptions
- The rock frame is fluid-independent (Gassmann's low-frequency assumption) — valid for well-log and seismic frequencies in most conventional reservoirs but can break down in low-permeability or heavy-oil rocks
- Pore fluid is uniformly distributed (patchy saturation effects are not modeled)
- Dry-frame moduli (Kdry, G) are accurately known or have been correctly back-calculated from the in-situ saturated-rock measurements
Limitations
- Gassmann assumes full equilibration of pore pressure at the frequency of interest — high-frequency laboratory or some sonic measurements may show dispersion not captured here
- Highly viscous fluids (heavy oil, tar) and very low permeability rocks can violate the fluid-independent-frame assumption, biasing results
- Mineral modulus Kmin is sensitive to mineralogy assumptions — mixed lithologies require an effective Kmin from a mineral mixing model
Use Cases
- → 4D seismic (production monitoring): Predict the expected impedance change as reservoir fluid is depleted or replaced (e.g., waterflood front, gas cap expansion) for time-lapse seismic interpretation.
- → AVO forward modeling: Generate synthetic Vp, Vs, and density logs for different fluid scenarios to forward-model expected AVO response.
- → Fluid substitution for DHI analysis: Test whether an observed amplitude anomaly is consistent with a hydrocarbon-fluid Gassmann substitution from a brine baseline.
- → Rock physics template construction: Build Vp/Vs vs. AI templates spanning a range of fluid and porosity scenarios for use in seismic interpretation workflows.
Related Calculations
Region Notes
Permian Basin
Wolfcamp mineralogy is quartz/clay/carbonate mixed — use an effective Kmin around 50-65 GPa (Voigt-Reuss-Hill average) rather than pure quartz (37 GPa) for realistic Ksat predictions.
Gulf of Mexico
Unconsolidated Miocene/Pliocene sands have low Kdry (often <5 GPa), making the Gassmann fluid term ΔKf a large fraction of Ksat — brine-to-gas substitution produces some of the strongest 4D signals in GOM.
North Sea
Brent Group sands are commonly used as Gassmann case studies for waterflood monitoring (e.g., Schiehallion, Magnus); typical Kdry≈8-15 GPa with quartz-dominated Kmin≈37-40 GPa.
References
Primary source
Primary source
Need geoscience support? BauerCalc is just one of the tools we build. If you're looking for independent expertise in well planning, geosteering, reservoir characterization, or opportunity evaluation, visit BauerSubsurface.com.