I have a seismic dataset that contains a depth-migrated seismic (amplitude) volume in depth and then the exact same volume converted to time using the migration velocities. Unfortunately, I do not have the velocity volume used to do the conversion.
How can I extract the velocity volume using only these two (amplitude) volumes?
Since the signals are exactly the same, just stretched and squeezed, I believe the solution will require something like dynamic time warping but I don't have enough experience with this to make it work.
Since each trace in the volume is converted independently, this reduces to a 1D problem. An optimal solution for a single trace can be used to solve the entire 3D volume (unless there is a lot of error and then we can have lateral instability, but I can live with that).
I will eventually need to code it in MATLAB for application on sgy-formatted data.