[Seisunix] How to create a 3d seismic data

1,023 views
Skip to first unread message

Glauco Lira

unread,
Dec 13, 2011, 6:31:48 PM12/13/11
to seis...@mailman.mines.edu
Dear

How to create a 3D seismic data, and view this data cube? If the "Su" does not make, which free software to obtain a 3d data?

Best regards

Karl Schleicher

unread,
Dec 14, 2011, 9:49:56 AM12/14/11
to glp...@yahoo.com.br, seismic unix
I have the teapot dome 3D dataset from RMOTC.  You can contact them and get FTP access.  Let me know if you are interested and cannot find them in google.

I think the migrated volume has the same number of xline traces for each inline, so after running segyread you can view it in suxmovie.  The teapot data has iline and xline in header locations su use for their own extensions.  I moved these bytes into another location when running segyread.

I am not sure what you are trying to do.  Processing 3D seismic data with open source software is a broad request.  I am interested, but just getting started.  I think some software updates might be a good idea.

Karl Schleicher
 


Date: Tue, 13 Dec 2011 15:31:48 -0800
From: glp...@yahoo.com.br
To: seis...@mailman.mines.edu
Subject: [Seisunix] How to create a 3d seismic data


Dear

How to create a 3D seismic data, and view this data cube? If the "Su" does not make, which free software to obtain a 3d data?

Best regards

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

muhammad toqeer

unread,
Dec 14, 2011, 10:44:05 AM12/14/11
to seis...@mailman.mines.edu
Dear Glauco Lira
you can generate and process the 3D seismic data by <<Madagascar>> which is open source seismic data processing utility like SU.
You can find this software and wiki at:
you can use << sfzomig3 >> function to model/migrate 3D data.
Regards
Toqeer



Glauco Lira

unread,
Dec 14, 2011, 1:32:55 PM12/14/11
to Karl Schleicher, seismic unix
Dear Karl Schleicher,

I am making a algoritm for 3d prestack depth migration and i need a cube data and a estimated velocity model to migrate. During the construction phase of the algorithm, I think that the better way is to build a velocity model and obtain a 3d data then to migrate and verify if teh amplitudes are in correct depth. The procedure is to calibrate the algoritm.

I have interest in teapot dome 3D dataset from RMOTC and I not found on the Internet. I'll be grateful if you provide it.

I tried to use "suxmovie". Look:

$ segyread tape=Seis-1shot.sgy vebose=1 endian=0 | segyclean > Seis-1shot.su
$ surange<Seis_Two-Hor-Layer-1shot.su
121 traces:
tracl    1 121 (1 - 121)
tracr    1 121 (1 - 121)
fldr     1
tracf    1 121 (1 - 121)
ep       1
cdp      1
cdpt     1
trid     1
gdel     1 11 (1 - 11)
sdel     1 11 (1 - 11)
scalel   1
scalco   1
sx       2500
sy       2500
gx       0 5000 (0 - 5000)
gy       0 5000 (0 - 5000)
counit   1
ns       1001
dt       4000
$ suxmovie<Seis-1shot.su n1=1001 ntr=121 n2=11 n3=11

But I saw only 3d data in 2d visualisation.

Another manner is using (tip of muhammad toqeer):

$ supscube < Seis-1shot.su > Seis-1shot.ps

But still do not know how to change the view.

Thanks all for contribution.


De: Karl Schleicher <k_schl...@hotmail.com>
Para: glp...@yahoo.com.br; seismic unix <seis...@mailman.mines.edu>
Enviadas: Quarta-feira, 14 de Dezembro de 2011 11:49
Assunto: RE: [Seisunix] How to create a 3d seismic data

Karl Schleicher

unread,
Dec 15, 2011, 6:33:26 PM12/15/11
to seismic unix
Glauco Lira,

To get the Teapot Dome data go to:
http://www.rmotc.doe.gov/datasets.html
In the first paragraph there is a "contact us" link you can click.  Fill out the form and they will send you an account and password you can use to ftp the data.  I think the filtered segy is in a Poststack directory.  When you get it you can run:

segyread bfile=checkstack.bin hfile=checkstack.asc tape=filt_mig.sgy endian=0 conv=1 | segyclean > checkstack.su
< checkstack.su suxmovie n2=188 loop=1 fframe=1 dframe=1 title="filtered migration %g" perc=99

This will make a movie of the data.  If you want a 3d perspective rendering, I do not know how to do that.

Terms of used for the teapot dome is at the bottom of the web page listed above.  It is very generous.

I do not understand what the Seis_1shot.sgy data is supposed to be, so I am lost on your processing scripts.  I never heard of supscube, so I cannnot help you there either.

Karl



Date: Wed, 14 Dec 2011 10:32:55 -0800
From: glp...@yahoo.com.br
Subject: Re: [Seisunix] How to create a 3d seismic data
To: k_schl...@hotmail.com
CC: seis...@mailman.mines.edu

Reginald Beardsley

unread,
Dec 15, 2011, 10:45:30 PM12/15/11
to seis...@mailman.mines.edu

cmovie might do much of what you want.  It's been broken  for some time because of the switch from 8 bit to 24 bit  graphics in X11 between R5 & R6.  However, I was able to get the code  transferred to Mines when Unocal got bought by Chevron.

Might be time to put that up on the CWP site and see if someone will volunteer to fix it.  I know quite a lot about the code as an earlier version was my introduction to C & Unix as a summer intern at ARCO's Plano Research Center (aka PRC) in 1989.  In addition I added the ability to do color overlay of velocity information and fixed a lot of bugs  at Unocal 10-12 years later.  It felt a lot like a complete rewrite at the time.  Rick is very sharp, but  given to being too clever by half.

I am NOT volunteering to fix it.  However, if someone wants to really learn X11  I will be happy to explain what the problem is.  It  requires a very small amount of programming and a lot of reading and research. 

For those who already know X11, it's a get default visual problem.  I *really* don't want to remember what the exact call is.   cmovie  expects an 8 bit default visual which is a rather ancient concept.  It can be fixed by feeding arguments to the X server, but a proper fix would be to modify cmovie.  Among other things this would greatly improve the dynamic range of the display when doing color overlays.

It's really worth the work. It's a very useful program.  I just don't like doing the same thing more than once, and especially dislike doing it 3 times.  Though for suitably large sums I can be persuaded ;-)

Have Fun!
Reg

Note:  There are a variety of versions floating around.  So there may be a working version already, though probably w/o the color overlay.  That was a feature in another pair of Rick's programs, vmerge & vdisplay.

muhammad toqeer

unread,
Dec 16, 2011, 4:26:11 AM12/16/11
to seis...@mailman.mines.edu
Dear Glauco Lira
Suxmovie is to view the movie style visualization of the 3D data; So when you will run the << suxmovie>> command you will see all the inlines and cross lines of the data that are passing rapidly as just seems in the movie, so you will not see the 3D cube.
You can output the ps graph of your data in 3D by running <<supscube>>. Suppose you have a cube of total 900 traces, 30 inlines and 30 cross lines, then you can run <<supscube >> to get the 3D data in ps formate

supscube<data.su ntr=900 n2=30 n3=30 >tt.ps 

you have to tune the other parameters for specific color combinations and axis labels.
Other thing is that you data ::
surange<Seis_Two-Hor-Layer-1shot.su
gives cdp=1 and cdpt=1. Are you sure this is 3D data?
If you have installed  <<Madgascar>> then try to view the data by <<sfgrey3>>.(Offcourse you have to convert the .segy or .su data into .rsf formate)
Let me know if you have success in data visualization.
Regards
Toqeer

Glauco Lira

unread,
Dec 20, 2011, 7:23:22 AM12/20/11
to seis...@mailman.mines.edu
Karl Schleicher,

Thank you, your tips are already of great value

Regards

Glauco Lira.

smart...@gmail.com

unread,
May 21, 2013, 4:37:01 AM5/21/13
to osdeve_mirror_ge...@googlegroups.com, seis...@mailman.mines.edu, mtoq...@gmail.com
HI,

How are you/?
Do you have any way to visualize the 3D cube SU format seismic data/. in any software without converting to segy? because when i convert from SU to SEGY, inlines and crosslines are needed to be remapped which I dont know how to do.

Other point that is it possible to repair the geometry like sorting in inlines and crosslines again after stacking the data?

waiting for help and reply

Regards

Khawar
Reply all
Reply to author
Forward
0 new messages