Skip to content

Instantly share code, notes, and snippets.

@bradlipovsky
Created March 15, 2022 21:55
Show Gist options
  • Save bradlipovsky/5b43bfacd3f76fdb000ddf48bc9b0524 to your computer and use it in GitHub Desktop.
Save bradlipovsky/5b43bfacd3f76fdb000ddf48bc9b0524 to your computer and use it in GitHub Desktop.
dV Back of the envelope
dx_apparent = 14 # Known source-receiver distance
vtrue = 6e3 # Known p-wave speed
t0 = 380/48e3 # Measured travel time
dt = (5/48e3) # Measured travel time change
dx_true = t0 * vtrue # Actual propagation path >> source-receiver distance
print('True propagation path length: %f'%dx_true)
v1 = dx_true/t0
v2 = dx_true/(t0-dt)
dv = (v2 - v1)/v1
print('Velocity change: %f'%dv)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment