[Seisunix] How to creat a wedge spike model?

131 views
Skip to first unread message

WKurt

unread,
Feb 21, 2012, 8:28:42 AM2/21/12
to seis...@mailman.mines.edu
Hi, everyone,

I am a new user to SU, does anyone know how to creat a wedge spike model in SU? I need to use it as a model reflectivity profile.

Thanks a lot!


Wietze Eckhardt

unread,
Feb 21, 2012, 10:23:50 AM2/21/12
to WKurt, seis...@mailman.mines.edu
Use suhelp and sufind and search for something with "syn" for synthetic. That should give you a shortlist of modules. Then use the self doc to find out how to parameterize.

Greetings,
Wietze

For example:
>sufind syn

or

> suhelp | grep susyn

you need to enter a few returns and then you get:

> suhelp | grep susyn

suaddevent              suimp2d                 susyncz 
suea2df                 suimp3d                 susynlv 
sufctanismod            suimpedance             susynlvcw 
sufdmod1                sukdsyn2d               susynlvfti 
sufdmod2                sunhmospike             susynvxz 
sufdmod2_pml            sunull                  susynvxzcs 

>


And then read the self docs by typing the program name without arguments and choose the program that best fits your problem. For instance:

> susyncz
                                                                       
 SUSYNCZ - SYNthetic seismograms for piecewise constant V(Z) function  
           True amplitude (primaries only) modeling for 2.5D           
                                                                       
  susyncz > outfile [parameters]                                       
                                                                       
 Required parameters:                                                  
 none                                                                  
                                                                       
 Optional Parameters:                                                  
 ninf=4        number of interfaces (not including upper surface)      
 dip=5*i       dips of interfaces in degrees (i=1,2,3,4)               
 zint=100*i    z-intercepts of interfaces at x=0 (i=1,2,3,4)           
 v=1500+ 500*i velocities below surface & interfaces (i=0,1,2,3,4)     
 rho=1,1,1,1,1 densities below surface & interfaces (i=0,1,2,3,4)      
 nline=1        number of (identical) lines                            
 ntr=32        number of traces                                
 dx=10         trace interval                                          
 tdelay=0      delay in recording time after source initiation         
 dt=0.004      time interval                                           
 nt=128        number of time samples                                  
                                                                       
 Notes:                                                        
 The original purpose of this code was to create some nontrivial       
 data for Brian Sumner's CZ suite.                                     
                                                                       
 The program produces zero-offset data over dipping reflectors.
                                                                       
 In the original fortran code, some arrays had the index               
 interval 1:ninf, as a natural way to index over the subsurface
 reflectors.  This indexing was preserved in this C translation.       
 Consequently, some arrays in the code do not use the 0 "slot".
                                                                       
 Example:                                                              
        susyncz | sufilter | sugain tpow=1 | display_program           
                                                                       
 Trace header fields set: tracl, ns, dt, delrt, ntr, sx, gx            
>








_______________________________________________
seisunix mailing list
seis...@mailman.mines.edu
https://mailman.mines.edu/mailman/listinfo/seisunix
Unsubscribe: seisunix-u...@mailman.mines.edu


Robert Hardy

unread,
Feb 21, 2012, 10:23:12 AM2/21/12
to WKurt, seis...@mailman.mines.edu
Greetings,

There are several ways to create wedge models in su. Two methods I have used in the past are..

1. Use any of the modelling programs e.g. susynlv to create the wedge in a constant velocity background and then migrate it with sustolt to get the reflectivity model
2. Use suplane to create the wedge reflectivity and then convolve it with chosen wavelet e.g. ricker generated by suwaveform

Cheers

Rob

________________________________________
From: seisunix...@mailman.mines.edu [seisunix...@mailman.mines.edu] On Behalf Of WKurt [kurt_a...@hotmail.com]
Sent: 21 February 2012 13:28
To: seis...@mailman.mines.edu
Subject: [Seisunix] How to creat a wedge spike model?

Hi, everyone,

I am a new user to SU, does anyone know how to creat a wedge spike model in SU? I need to use it as a model reflectivity profile.

Thanks a lot!


Jinlong Yang

unread,
Feb 21, 2012, 3:23:44 PM2/21/12
to seis...@mailman.mines.edu
Hi everyone,

Does anybody know how to convert a binary file to video or movie file?
I can make a movie using xmovie or suxmovie to xwindow, but I don't
how to output the movie.

Thanks a lot.

Jinlong Yang

unread,
Feb 21, 2012, 3:58:27 PM2/21/12
to Mike Stone, seis...@mailman.mines.edu
Hi Mike,

Thanks for your reply. There is still some issues. The following is my test:

#create the data junk.su
suplane | suaddnoise sn=20 >> junk.su
suplane | suaddnoise sn=15 >> junk.su
suplane | suaddnoise sn=10 >> junk.su
suplane | suaddnoise sn=5 >> junk.su
suplane | suaddnoise sn=3 >> junk.su
suplane | suaddnoise sn=2 >> junk.su
suplane | suaddnoise sn=1 >> junk.su

#plot junk.su using suxmovie
suxmovie < junk.su n2=32 title="frame=%g" loop=1 &

#convert junk.su to ps file
supsmovie < junk.su n2=32 > junk.ps

#convert junk.ps to a movie file
Here I can't find the command imagemagick and I tried the command
convert, but it is not work.

convert junk.ps junk.avi

Could you help me to find where I am wrong and how to convert the ps
file to avi?
Thanks again.

Jinlong


2012/2/21 Mike Stone <mst...@act3solutions.com>:
> Jinlong,
>
> I prefer to use psmovie or supsmovie. The colorbar settings are a little
> different so be careful.
>
> After you output it to a ps file, use imagemagick. "Convert movie.ps
> movie.avi"
>
>
>
> Alternatively you can record the output with the walk-through here:
> http://www.seismicunix.com/w/Wiki_Guide_Plotting:Making_movies
>
>
>
> Mike Stone
> 832-768-0884
> MSt...@ACT3Solutions.com
>
>
>
> 2012/2/21 Jinlong Yang <uhj...@gmail.com>

John Stockwell

unread,
Feb 21, 2012, 4:59:08 PM2/21/12
to Jinlong Yang, seis...@mailman.mines.edu
You can read Jeff Barker's instructions for making Gif movies at:

http://bingweb.binghamton.edu/~jbarker/animations.html

likely you will have to install "ImageMagic" which does come as
an option on most Linux systems. You have not mentioned what your
operating system is.

If you want to convert the gif movie to an .avi format, then you
may install ffmpeg to do this.

-John

> Hi everyone,
>
> Does anybody know how to convert a binary file to video or movie file?
> I can make a movie using xmovie or suxmovie to xwindow, but I don't
> how to output the movie.
>
> Thanks a lot.
>
>
>
>

> ? 2012?2?21? ??7:28?WKurt <kurt_a...@hotmail.com> ???


>> Hi, everyone,
>>
>> I am a new user to SU, does anyone know how to creat a wedge spike model in
>> SU? I need to use it as a model reflectivity profile.
>>
>> Thanks a lot!
>>
>>
>>
>> _______________________________________________
>> seisunix mailing list
>> seis...@mailman.mines.edu
>> https://mailman.mines.edu/mailman/listinfo/seisunix
>> Unsubscribe: seisunix-u...@mailman.mines.edu
>>
>
> _______________________________________________
> seisunix mailing list
> seis...@mailman.mines.edu
> https://mailman.mines.edu/mailman/listinfo/seisunix
> Unsubscribe: seisunix-u...@mailman.mines.edu
>

John Stockwell | jo...@dix.Mines.EDU
Center for Wave Phenomena (The Home of Seismic Un*x)
Colorado School of Mines
Golden, CO 80401 | http://www.cwp.mines.edu/cwpcodes
voice: (303) 273-3049

Our book:
Norman Bleistein, Jack K. Cohen, John W. Stockwell Jr., [2001],
Mathematics of multidimensional seismic imaging, migration, and inversion,
(Interdisciplinary Applied Mathematics, V. 13.), Springer-Verlag, New York.

Glenn Reynolds

unread,
Feb 22, 2012, 6:34:48 AM2/22/12
to Jinlong Yang, seis...@mailman.mines.edu
ImageMagick is a separate package from SU, and is worth having. It's open source. 

Glenn Reynolds
(mobile)

ryana...@gmail.com

unread,
Feb 28, 2016, 6:36:51 AM2/28/16
to osdeve.mirror.geophysics.seisunix, kurt_a...@hotmail.com, seis...@mailman.mines.edu, HAR...@tcd.ie
Greetings,

What if I want to make.. let say homogeneous subsurface with single velocity. But there is on point with different (higher) velocity.
how to generate stacked seismic section for this kind of problem? I may have an idea with suplane and using len1=ntr but I think we don't use any velocity and suplane

Thank in advance
Reply all
Reply to author
Forward
0 new messages