From Moduli to Velocities

Part 1, Part 1: Elasticity and the Moduli

Learning objectives

  • Write the two velocity equations and identify which moduli each one feels
  • Use the unit convention in which moduli in GPa and density in g/cc give velocity directly in km/s
  • Explain why a pore-fluid change moves Vp and Vs differently, the germ of fluid discrimination
  • Trace how density in the denominator lets a stiffening fluid and a heavier rock fight over Vp

The Two Equations

With the two moduli in hand, the velocities are immediate. A P-wave both compresses and shears the rock as it passes, so it feels both moduli: VP=K+43μρV_P = \sqrt{\dfrac{K + \tfrac{4}{3}\mu}{\rho}}. An S-wave only changes shape, so it feels the shear modulus alone: VS=μρV_S = \sqrt{\dfrac{\mu}{\rho}}. Both carry the density ρ\rho in the denominator, because inertia resists the motion regardless of which stiffness is driving it. These two equations are the hinge of the whole course: composition and texture set KK, μ\mu, and ρ\rho on the right, and the seismic velocities fall out on the left.

A Convenient Set of Units

The equations look cleanest in a unit system that geophysics adopted for exactly this reason. Put the moduli in gigapascals and the density in grams per cubic centimeter, and the square root comes out directly in kilometers per second, with no conversion factor to remember, because GPa/(g/cc)\sqrt{\text{GPa}/(\text{g/cc})} is exactly km/s. Take mineral quartz as the standard check: K=36.6K = 36.6 GPa, μ=45\mu = 45 GPa, ρ=2.65\rho = 2.65 g/cc. Then VP=(36.6+43⋅45)/2.65=6.04V_P = \sqrt{(36.6 + \tfrac{4}{3}\cdot 45)/2.65} = 6.04 km/s and VS=45/2.65=4.12V_S = \sqrt{45/2.65} = 4.12 km/s. Those two numbers reappear in every part of this course; if a frame or fluid model ever predicts a clean quartz sand faster than mineral quartz itself, the model is wrong.

Moduli to velocitiesK = 36.6μ = 45ρ = 2.65GPa, GPa, g/ccVₚ6.04 km/sVₛ4.12 km/sVₚ carries K and the shear modulus; Vₛ carries the shear modulus alone.GPa over g/cc gives km/s directly: quartz reads Vp 6.04, Vs 4.12.

Why the Fluid Shows Up Twice

Now let a fluid enter the pores and watch the two velocities part ways. The fluid raises the bulk modulus KK, because it resists the squeeze, and it raises the density ρ\rho, because it has mass, but it leaves μ\mu exactly where it was. So VS=μ/ρV_S = \sqrt{\mu/\rho} can only fall: its numerator is fixed while its denominator grows, and a saturated rock always has a slightly lower shear velocity than the same rock dry. VPV_P is the subtle one, because KK and ρ\rho pull it in opposite directions: a stiffer KK speeds the P-wave up, a heavier ρ\rho slows it down, and which wins depends on the fluid. That is why the choice of fluid, not merely its presence, changes the seismic answer.

The two headline cases follow directly. Swapping brine for gas softens the saturated KK far more than it lightens ρ\rho, so the drop in KK wins and VPV_P falls sharply, the reason gas often lights up, while VSV_S actually edges up a little as the density drops. Swapping brine for oil nudges KK and ρ\rho by similar small amounts whose effects on VPV_P nearly cancel, so VPV_P barely moves and oil is close to invisible on velocity alone. That the same fluid reaches VPV_P through two competing channels and VSV_S through only one is why no single velocity tells the whole story. The cleanest way to read the pair is their ratio, which cancels the density and exposes the balance of the two moduli. That ratio, and the Poisson's ratio it is equivalent to, is the next section.

References

  • Mavko, G., Mukerji, T., & Dvorkin, J. (2009). The Rock Physics Handbook (2nd ed.). Cambridge University Press.
  • Simm, R., & Bacon, M. (2014). Seismic Amplitude: An Interpreter's Handbook. Cambridge University Press.

This page is prerendered for SEO and accessibility. The interactive widgets above hydrate on JavaScript load.