EW0210 SIOSEIS scripts - real-time brute stack




#! /bin/csh
#
#   Remember to run script "init_script" (sioseis-seisnet daemon)
#
if( $#argv < 2 ) then
    echo "Usage: stack line-name plot-direction(ltr/rtl)"
    exit 1
endif
set LINENO = $1
set DIR = $2

#   This should be run in /export/home2/scratch/ew0210/brute so the
#  sioseis tmp files do not conflict with other sioseis scripts

/home/henkart/bin/sioseis << eof
  
procs segddin prout weight geom header gather nmo mute stack
     diskoe filter avenor plot end

segddin
   ftr 241 ltr 480
   fcset 1 lcset 1
   format seisnet
   stack latest
   logpath /export/home3/seisnet/0210/logs/$LINENO.log
   end
end
avenor
   hold 300
   addwb yes sets .2 .7   1.2 1.7  2.2 2.7  8.2 8.7 end
end
gains
  addwb yes tgp 0 .8 1.5 .8 3 4 7 12 16 24 end
end
weight
    fno 0 lno 999999 twp 236 0 end
end
geom
  type 9 # Use realtime GPS
  fs 1 ls 999999  # all shot have the same parameters (preset)
  gxp 480 -180   # RESET the closest group only.
  ggx -12.5         # Used to extrapolate gxp!
# dfls 37.5   # ignored with type 9 
  dbrps 6.25 smear 6.25  
  rpadd 1000 end
end

mute
   fno 1 lno 999999
   addwb yes xtp 200 -.2  end
end
  
diskoa   # Write out the filtered stack file disk file
  opath /export/home3/seisnet/0210/stacked_lines/line$LINENO.stack-filter end
end
diskob 
# write every 50th shot to a "circular" file 
# remember that segdin limited the traces read (ftr/ltr/trinc)
    fno 1 lno 999999 noinc 50 rewind 1 
    opath /export/home3/seisnet/0210/shots/latest.shot end 
end
diskod 
# write out the gather
    fno 1 lno 999999 noinc 50 rewind 1 
    opath /export/home3/seisnet/0210/shots/latest.gather end 
end
diskoe   # Write out disk file
  opath /export/home3/seisnet/0210/stacked_lines/line$LINENO.stack end
end

prout
 fno 0 lno 99999 ftr 479 ltr 479 noinc 10 end
end

header
   fno 0 lno 9999999 ftr 1 ltr 9999 
   r50 r54 / 750.  # convert water depth to water time for addwb (mute)
   end
end

gather
#   maxtrs 90 maxrps 500 end
   maxtrs 50 maxrps 250 end   # half the streamer
end

nmo
# real time nmo, replace interpolation by RP to WB depth in Meters.
# If water depth changes by > 500 m, use previous value. Water-depth
# velocity functions derived from ESP5, interpolation by iso-velocity layering
 vtrkwb 500 stretc 1.50

fno 10 lno 10
vtp 1500 0.013 
 1700 0.884 
 2171 1.312 
 2486 1.687 
 3092 2.208 
 5733 8.588 
 6556 13.088 end
   
fno 500 lno 500
vtp 1500 0.667 
 1630 1.549 
 1977 1.978 
 2238 2.353 
 2775 2.873 
 5673 8.391 
 6788 13.391 end
   
fno 1000 lno 1000
vtp 1500 1.333 
 1564 2.216 
 1761 2.644 
 1932 3.019 
 2089 3.332 
 5413 7.470 
 7181 13.470 end
   
fno 2000 lno 2000
vtp 1500 2.667 
 1521 3.549 
 1593 3.978 
 1664 4.353 
 1734 4.665 
 4410 7.769 
 7122 14.269 end
   
fno 2500 lno 2500
vtp 1500 3.333 
 1502 3.627 
 1521 3.885 
 1546 4.135 
 1555 4.260 
 3312 6.329 
 7244 13.829 end
   
fno 3000 lno 3000
vtp 1500 4.000 
 1501 4.294 
 1504 4.408 
 1505 4.483 
 1564 4.858 
 2900 6.551 
 7064 14.051 end
end

filter
    pass 5 60 ftype 0 dbdrop 48 end
end

 agc 
   winlen 1. center .1 end 
 end 

 plot 
   dir $DIR
  scalar -1
  tlines 0.5 1 nibs 7224 ann gmtint anninc 5 
  def 0.05 trpin 125 wiggle 0 
   nsecs 16 vscale 1.25 clip .03
  opath /export/home3/seisnet/0210/stack_plots/$LINENO.atlantek
  end 
end

end
eof