Seismic Wave Velocity Calculator
Calculate P-wave and S-wave velocities from elastic moduli and density. Convert between bulk/shear modulus and Young modulus inputs.
Seismic Wave Velocity Calculator
Calculator
Adjust values & calculateEnter your values below. Every result is computed in your browser โ no data is sent to any server.
Formula: Vp = sqrt((K + 4G/3) / rho) | Vs = sqrt(G / rho)
Worked example โ Vp = 5828 m/s, Vs = 3365 m/s, Vp/Vs = 1.732
Formula
Vp = sqrt((K + 4G/3) / rho) | Vs = sqrt(G / rho)
P-wave velocity depends on bulk modulus K, shear modulus G, and density rho. S-wave velocity depends only on shear modulus and density. These fundamental relationships connect elastic properties of rocks to seismic wave propagation speeds.
Worked Examples
Example 1: Granite P and S-wave Velocity
Problem:Calculate Vp and Vs for granite with K = 50 GPa, G = 30 GPa, density = 2650 kg/m^3.
Solution:Vp = sqrt((50e9 + 4*30e9/3) / 2650) Vp = sqrt((50e9 + 40e9) / 2650) Vp = sqrt(3.396e7) = 5828 m/s Vs = sqrt(30e9 / 2650) = sqrt(1.132e7) = 3365 m/s
Result:Vp = 5828 m/s, Vs = 3365 m/s, Vp/Vs = 1.732
Example 2: Limestone from Young Modulus
Problem:Find velocities for limestone: E = 60 GPa, Poisson ratio = 0.28, density = 2700 kg/m^3.
Solution:K = 60e9 / (3*(1-0.56)) = 45.45 GPa G = 60e9 / (2*(1.28)) = 23.44 GPa Vp = sqrt((45.45e9 + 4*23.44e9/3) / 2700) = 5395 m/s Vs = sqrt(23.44e9 / 2700) = 2946 m/s
Result:Vp = 5395 m/s, Vs = 2946 m/s
Frequently Asked Questions
What determines seismic wave velocity?
Seismic wave velocity depends on the elastic properties and density of the medium. P-wave velocity Vp depends on the bulk modulus K (resistance to compression), shear modulus G (resistance to shearing), and density: Vp = sqrt((K + 4G/3)/rho). S-wave velocity depends only on shear modulus and density: Vs = sqrt(G/rho). Higher elastic moduli increase velocity, while higher density decreases it. In practice, the modulus effect dominates, so denser rocks generally have higher velocities because they also tend to be stiffer.
What is the difference between P-waves and S-waves?
P-waves (primary or compressional waves) involve particle motion parallel to the direction of wave propagation, alternately compressing and expanding the material. S-waves (secondary or shear waves) involve particle motion perpendicular to the propagation direction. P-waves travel faster than S-waves because they involve both bulk and shear deformation, while S-waves involve only shear. S-waves cannot propagate through fluids (liquids and gases) because fluids have zero shear modulus, which is why the outer core of the Earth blocks S-waves.
What are typical seismic velocities for common rocks?
Common P-wave velocities include: water 1500 m/s, unconsolidated sand 200-800 m/s, sandstone 2000-4500 m/s, limestone 3500-6500 m/s, granite 5000-6500 m/s, basalt 5500-6500 m/s, and the Earth mantle 7800-13700 m/s. S-wave velocities are typically 50-70% of P-wave velocities for solid rocks. Velocities generally increase with depth due to increased confining pressure, which closes cracks and increases stiffness. Fluid saturation increases Vp but has little effect on Vs.
How is Poisson ratio related to wave velocities?
Poisson ratio (nu) is directly calculable from the Vp/Vs ratio: nu = ((Vp/Vs)^2 - 2) / (2*(Vp/Vs)^2 - 2). For most rocks, Poisson ratio ranges from 0.1 to 0.45. A ratio of 0.25 corresponds to Vp/Vs = 1.732. Fluid-saturated rocks tend to have higher Poisson ratios (0.3-0.4) than dry rocks. Gas-bearing sands often show anomalously low Vp/Vs ratios and Poisson ratios, making this measurement valuable for hydrocarbon detection in seismic exploration.
How are Young modulus and Poisson ratio converted into bulk and shear modulus?
The isotropic elastic constants are interchangeable in pairs. Bulk modulus follows from K = E / (3 x (1 - 2 x nu)) and shear modulus from G = E / (2 x (1 + nu)), after which the standard velocity expressions apply. Note that K blows up as nu approaches 0.5, the incompressible limit, which is why Seismic Wave Velocity Calculator rejects a Poisson ratio of 0.5 or above. Use this route when a laboratory report gives static or dynamic Young modulus rather than the bulk and shear pair.
How does pore fluid saturation change Vp and Vs differently?
Filling pore space with brine raises the saturated bulk modulus because the fluid resists compression, so Vp increases. Shear modulus is unchanged by a fluid that cannot support shear, so the numerator of Vs stays fixed while the bulk density rises, and Vs actually drops slightly. Gassmann's 1951 relations formalise this fluid substitution. The divergent response is exactly why Vp/Vs and Poisson ratio are diagnostic: a gas-filled sand has a low fluid modulus and shows a distinctly depressed Vp/Vs.
How do I get P-wave velocity from a sonic log?
A sonic tool records interval transit time, or slowness, rather than velocity, so you invert it. With slowness DT in microseconds per foot, velocity in feet per second is 1,000,000 / DT, meaning a reading of 100 microseconds per foot corresponds to 10,000 ft/s, about 3048 m/s. With slowness in microseconds per metre the same reciprocal gives metres per second directly. Pairing that Vp with a bulk density log lets you compute acoustic impedance and, where a dipole shear log exists, the Vp/Vs ratio.
What is the difference between interval velocity, average velocity and RMS velocity?
Interval velocity is the true propagation speed within a single layer, the quantity these elastic formulas give. Average velocity is total depth divided by total one-way travel time. RMS velocity is a root-mean-square quantity estimated from the curvature of reflection moveout on a common midpoint gather and used for normal moveout correction and stacking. Interval velocities are recovered from stacking velocities through the Dix equation, which differences the squared RMS velocities weighted by two-way time between reflectors.
Why does seismic velocity increase with depth of burial?
Rising confining pressure closes microcracks and compliant grain contacts, and compaction reduces porosity, so both the bulk and shear moduli stiffen faster than density increases. The effect is strongest in the first few hundred metres, where velocity can rise steeply, then flattens once most crack porosity is closed. Overpressured intervals reverse the trend: high pore pressure holds grains apart, lowering the effective stress and producing a velocity drop that is used to predict drilling hazards ahead of the bit.
What is seismic anisotropy and why does a single velocity value not capture it?
Layered shales and aligned fractures make velocity depend on propagation direction, so a rock has no single Vp. The common case is transverse isotropy with a vertical symmetry axis, where horizontal propagation is faster than vertical. Thomsen's 1986 parameters epsilon, delta and gamma quantify the departure from isotropy, and shales commonly show 10 to 30 percent anisotropy. Ignoring it causes mis-tied well depths, under-migrated images and moveout that a single-velocity hyperbola cannot flatten at long offsets.
References
Background & Theory
History
Reviewed for accuracy by Daniel Agrici, Founder & Lead Developer ยท Editorial policy
Related Calculators
๐งฎPlate Motion Velocity Calculator
Calculate plate motion velocity with inputs, formulas, and instant results.
๐งฎSeismic Attenuation Calculator
Calculate seismic attenuation with inputs, formulas, and instant results.
๐งฎSeismic Energy Release Calculator
Calculate seismic energy release with inputs, formulas, and instant results.
๐งฎSeismic Intensity Calculator
Calculate seismic intensity with inputs, formulas, and instant results.
๐งฎSeismic Moment Tensor Components Calculator
Calculate seismic moment tensor components with inputs, formulas, and instant results.
๐งฎSeismic Reflection Coefficient Calculator
Calculate seismic reflection coefficient with inputs, formulas, and instant results.
๐งฎGroundwater Flow Velocity Calculator
Calculate groundwater flow velocity with inputs, formulas, and instant results.
๐งฎGlacial Flow Velocity (glenโs Law) Calculator
Calculate glacial flow velocity (glenโs law) with inputs, formulas, and instant results.