A later stack script
Changes:
1) stacking 120 of 480 traces (2km of 6km streamer)
2) A new velocity function with a low velocity zone.
3) New geom option type 9 - binning by realtime lat/long
4) New gains parameter tadd used with subwb yes
#!/bin/csh -f
if( $#argv != 2 ) then
echo "Usage: stack line-number shot_spacing (in meters)"
exit 1
endif
set LINENO = $1
set DX = $2
mkdir /export/home/EW0008/Realtimestack/Line$LINENO
chmod a+rw /export/home/EW0008/Realtimestack/Line$LINENO
#./sioseis << eof
/home/henkart/bin/sioseis-fast << eof
#./sioseis-fast << eof
procs segdin geom prout gather nmo mute stack
diskoe filter gains plot end
segdin
ffilen 99999
fcset 1 lcset 1
ftr 361 ltr 480 ! drop long ranges (1-360 ~= 2km to 4km)
stime 3 secs 4
offline yes # eject after the rewind after EOT
iunit 4 end
end
tredit
limits -2000 0 kill outside lhdr 10 # kill 2km to 6km
fno 0 lno 999999 end
end
prout
fno 0 lno 999999 ftr 480 ltr 480 noinc 10 end
end
weight
fno 0 lno 9999999 twp 51 0 end
end
geom
type 9 # The new dfls from the realtime nav block
# type 2 # increment the shot location based on the shot number
fs 1 ls 999999 # all shot have the same parameters (preset)
ggx -12.5 # Used to extrapolate gxp!
# dfls $DX dbrps 6.25 smear 6.25
dfls 37.5 dbrps 6.25 smear 6.25
gxp 480 -184.25 ggx -12.5 dfls 37.5 dbrps 6.25 smear 6.25
bgp 1 6 2 18 3 30 4 42 5 66 6 90 7 114 8 138 9 162 10 186 11 210
12 234 13 258 14 282 15 306 16 330 17 354 18 378 19 402 20 426
21 438 22 459 23 462 24 474
rpadd 1000 end
end
mute
fno 1 lno 999999
addwb yes xtp 2000 -.03 2500 .6 6000 .7 end
end
plot
nsecs 4
scalar 2.E-04
tlines 0.1 0.5 1 nibs 7224
def 0.150 clip .1 trpin 75 wiggle -40 vscale 5
ann gmtint anninc 5
opath /export/home/EW0008/Realtimestack/Line$LINENO/siopltfil.Line$LINENO
end
end
diskoe # Write out disk file
opath /export/home/EW0008/Realtimestack/Line$LINENO/stack.line$LINENO end
end
gather
maxtrs 90 maxrps 500 end
# maxtrs 60 maxrps 500 end
end
filter
# ftype 0 dbdrop 48 pass 20 120 end
ftype 0 dbdrop 48 pass 40 220 end
end
gains
type 3 alpha 4 subwb yes tadd 3 end
end
nmo
addwb yes stretc 1.
vtp 1510 0 1530 .3 1550 .4 1510 .5 1560 .6 1580 .9 1650 1.5 end
end
end
eof
Return to EW0008 example.
Return to SIOSEIS examples.
Go to the list of seismic processes. Go to SIOSEIS introduction.