JohnWStockwellJr / SeisUnix

The CWP/SU: Seismic Un*x Package - a free open seismic processing, research, and educational software package. Please seek distribution gzipped tar files at https://wiki.Seismic-Unix.org not the Github version.
https://wiki.Seismic-Unix.org
Other
262 stars 124 forks source link

suztot issues? #64

Closed mithcell closed 2 years ago

mithcell commented 4 years ago

Anyone had luck using this module successfully? I get a result with the correct sample rate but the stretching hasn't worked. Perhaps the interpolation has a bug....

JohnWStockwellJr commented 4 years ago

The program works just fine. It is for stretching data in depth to time. It is the inverse of “suttoz”

John Stockwell Visit: https://wiki.Seismic-Unix.org The new home of Seismic Un*x

On Feb 8, 2020, at 1:48 AM, mithcell notifications@github.com wrote:

 Anyone had luck using this module successfully? I get a result with the correct sample rate but the stretching hasn't worked. Perhaps the interpolation has a bug....

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or unsubscribe.

mithcell commented 4 years ago

Thank you for the reply.

I find that I can convert a v(z) model to v(t) using: velconv intype=vintz outtype=vintt nz=300 dz=10 dt=0.004 nt=1001 <marmousi.bin | suaddhead ns=1001 | suximage

Now that I know what this model will look like in time, I try to use suztot but can't acheive the desired result: suztot vfile=marmousi.bin <marmousi.su nt=1001 dt=0.004 verbose=1 | suaddhead ns=1001 | suximage

suztot: Input: suztot: number of depth samples = 300 suztot: depth sampling interval = 10 suztot: first depth sample = 0

suztot: Output: suztot: number of time samples = 1001 suztot: time sampling interval = 0.004 suztot: first time sample = 0

Any suggestions?

Thank you

JohnWStockwellJr commented 4 years ago

When it says v(t) or v(z) it means that the velocity profile is one dimensional. While it is possible, with care, to use velconv to convert a v(z,x) file to a v(t,x), and vice versa, the programs suttoz and suztot do not support v(t,x) or v(z,x) files. These support only v(t) or v(z) velocity profiles.

For your Marmousi profiles, the best you could do would be to use suaddhead to put su headers on the Marmousi models and use sustack to average the models to create average v(z) and v(t) files. These could then be used with suttoz and suztot.

John Stockwell Visit: https://wiki.Seismic-Unix.org The new home of Seismic Un*x

On Feb 8, 2020, at 9:22 PM, mithcell notifications@github.com wrote:

 Thank you for the reply.

I find that I can convert a v(z) model to v(t) using: velconv intype=vintz outtype=vintt nz=300 dz=10 dt=0.004 nt=1001 <marmousi.bin | suaddhead ns=1001 | suximage

Now that I know what this model will look like in time, I try to use suztot but can't acheive the desired result: suztot vfile=marmousi.bin <marmousi.su nt=1001 dt=0.004 verbose=1 | suaddhead ns=1001 | suximage

suztot: Input: suztot: number of depth samples = 300 suztot: depth sampling interval = 10 suztot: first depth sample = 0

suztot: Output: suztot: number of time samples = 1001 suztot: time sampling interval = 0.004 suztot: first time sample = 0

Any suggestions?

Thank you

— You are receiving this because you commented. Reply to this email directly, view it on GitHub, or unsubscribe.