I am currently processing some seismic data for an experiment and I am having issue with using the program sufbpickw from http://www.cwp.mines.edu/cwpcodes/:
"Potash SU: SU style codes developed by Balazs Nemeth"
I am trying to run this program to do first break picking for my experimental survey. I can use the example file, and the program works just fine using the commands below:#visualize the seismic image before first arrival pickingsugain wagc=0.1 agc=1 < 71799offset_catted.su | suximage
#first arrival picking#visualize the first arrival picking resultsugain wagc=0.1 agc=1 < name_of_output.su | suximageShannon Chollett, M.S. Candidate
However, when i try to use my file, dry100211AVGshift.su, the program won't run and it gives me this error:sufbpickw: ealloc.c: malloc failed (Cannot allocate memory)
My file is created in a process like the one shown below:#Convert ascii to binary$infile=$file.'.txt';$num=$num1/$stacks;system("a2b <$infile >$file.bin n1=1");#Add SEG-Y headerssystem("suaddhead ns=$num <$file.bin | sushw key=dt a=13 > $file.su");#Display Datasystem("sufilter f=100,200,9000,13000 <$file.su | sugain agc=1 wagc=.001 | suximage");I am wondering why my file gives the error and what I can do to make it work? Also, is there another first break program in the seismic unix package that I could use instead of this one?Thank you for your time and effort,
Craft & Hawkins Department of Petroleum Engineering
Louisiana State University
Patrick F. Taylor Hall, Room 1413
Baton Rouge, LA 70803
Cell: 208-794-9296
1) If you get:
surange < 71799offset_catted.su
surange: fgettr.c: on trace #1, number of samples in header (41025)
differs from number for first trace (8)
then you need to put the file through: suoldtonew
suoldtonew < 71799offset_catted.su > 71799offset_catted_new.su
suoldtonew: converted 72 traces to XDR format
2) If surange looks like this, then the data are ok
surange < 71799offset_catted_new.su
72 traces:
fldr 1000
tracf 1 72 (25 - 72)
offset 10 223 (223 - 10)
ns 2048
dt 500
grnlof 0 1 (0 - 1)
3) you can view them with
suxwigb < 71799offset_catted_new.su perc=99
4) Now with the autopicker:
./sufbpickw < 71799offset_catted_new.su window=.01 | suxwigb perc=99
you can see the first arrivals, and you could use sumax to get the
values
./sufbpickw < 71799offset_catted_new.su window=.01 | sumax mode=max verbose=1
5) now try it for the dry
On my system I have to do:
suuoldtonew < dry100211AVGshift.su > dry100211AVGshift_new.su
suxwigb < dry100211AVGshift_new.su perc=99
sufbpickw < dry100211AVGshift_new.su window=.0005 | suxwigb perc=99
You have choose the window so that it is realistic for the data. The
default window is .01 sec. The dry data do not go to .01 sec.
-John
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.
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.