Gassmann Fluid Substitution Calculator: Rock Physics Vp and Vs Prediction

Predict how Vp and Vs change when a different pore fluid replaces the original one using Gassmann’s equation, model a multi-layer fluid substitution profile, or solve backward for the fluid bulk modulus a target Vp requires. Everything updates live.

How Gassmann’s Equation Separates the Frame From the Fluid

Gassmann’s equation states that the saturated rock bulk modulus depends on the dry frame modulus, the mineral modulus, the fluid modulus, and porosity, while the shear modulus stays exactly the same regardless of which fluid fills the pores. This separation is the key to fluid substitution, since it means the frame properties can be backed out from a measurement at one fluid saturation and then combined with a different fluid modulus to predict the rock’s behavior at a new saturation, all without ever measuring the second fluid case directly.

The full workflow runs in three steps: first invert the measured Vp, Vs, and density at the initial fluid saturation to get the dry frame modulus and shear modulus, second apply Gassmann’s equation forward with the new fluid’s bulk modulus to get the new saturated bulk modulus, and third combine that with the new fluid’s density to predict the new Vp and Vs.

Why Vp Moves More Than Vs

Because shear modulus is fluid independent, Vs changes only through the bulk density term in the denominator, typically a small effect. Vp responds to both the bulk modulus change, often large between a liquid and a gas, and the same density change, so the two effects can partly offset or reinforce each other depending on the rock’s stiffness. In stiffer, more consolidated rocks, the frame modulus dominates and a brine to gas substitution produces only a modest Vp drop, while in softer, more unconsolidated sands the fluid modulus contribution is large relative to the frame and the Vp drop from brine to gas can be dramatic, which is exactly the contrast that produces a bright seismic amplitude anomaly over a gas sand.

Building a Multi-Layer Substitution Profile

Real intervals rarely have uniform porosity or frame stiffness from top to bottom, so a multi-layer profile backs out the dry frame modulus independently at each depth from that depth’s own measured velocities, rather than applying a single average frame modulus across the whole interval. This matters because the dry frame modulus derived from a stiffer, lower porosity layer will respond differently to the same fluid substitution than a softer, higher porosity layer just above or below it.

Solving Backward From a Target Velocity

Since shear modulus and bulk density fix Vs2 once porosity and mineral properties are set, a target Vp2 pins down the required saturated bulk modulus directly, and inverting Gassmann’s equation from there solves for the fluid bulk modulus that specific target requires, no iteration needed. This is useful for checking whether an observed amplitude anomaly is consistent with a particular fluid type, or for testing how sensitive a velocity prediction is to fluid modulus uncertainty.

Common Mistakes

Applying a single dry frame modulus across a layer with varying porosity instead of backing it out independently at each depth, which can meaningfully bias the predicted velocities away from what the actual frame stiffness would produce.

Forgetting that Gassmann’s equation assumes a fully saturated, homogeneous, isotropic rock with fluid in pressure communication throughout the pore space, assumptions that break down at very high frequencies or in rocks with strong patchy saturation.

Mixing up bulk modulus and density effects when reasoning about a fluid substitution result, since the two act through completely different parts of the velocity equation and a small net Vp change can hide a large modulus change mostly offset by a large density change.

Using dry frame properties derived from an unrelated rock or a generic empirical relation instead of backing them out from the interval’s own measured velocities whenever those measurements are available.

Who Uses This

Rock physicists and geophysicists predicting how seismic velocities respond to a hypothetical fluid change ahead of amplitude versus offset or direct hydrocarbon indicator analysis.

Reservoir engineers checking whether a fluid substitution result is consistent with observed time lapse seismic changes during production.

Students and early career geoscientists learning how Gassmann’s equation separates frame and fluid effects on seismic velocity.

Analysts building multi-layer rock physics models to test how substitution results vary with porosity and frame stiffness across an interval.

Try It Above

Use the quick convert tab for a single fluid substitution calculation from measured Vp, Vs, and density, move to the layered model tab to see how the substitution plays out across several depths, or use the solve target Vp tab when a target velocity is set and you need the fluid modulus it requires.

Related Calculators and References

Continue building your workflow with the Pressure Gradient Calculator. For background on the formulas used above, see Wikipedia: Gassmann’s equation.