App-SeismicUnixGui
view release on metacpan or search on metacpan
lib/App/SeismicUnixGui/big_streams/iApply_mute.pm view on Meta::CPAN
$suxwigb->clear();
$suxwigb->defaults( 'iSurf4', 'top_middle' );
$suxwigb->key('tracl');
$suxwigb->title( $iApply_mute->{_file_in} );
#$suxwigb-> pmin($CFG->{sutaup}{1}{pmin});
#$suxwigb-> dp($dp);
$suxwigb->box_width(400);
$suxwigb->box_height(500);
#$suxwigb-> percent(99.9);
$suxwigb->clip(1);
$suxwigb->windowtitle( quotemeta('Muted Taup Data') );
$suxwigb[1] = $suxwigb->Step();
=head2
DEFINE FLOW(S)
=cut
lib/App/SeismicUnixGui/big_streams/iSelect_tr_Sumute.pm view on Meta::CPAN
=cut
$suxwigb->clear();
$suxwigb->defaults( 'iSurf4', 'top_middle' );
$suxwigb->key('tracl');
$suxwigb->title( $iSelect_tr_Sumute->{_file_in} );
# $suxwigb-> pmin($CFG->{sutaup}{1}{pmin});
# $suxwigb-> dp($dp);
# $suxwigb-> picks($sutaup_outbound_pickfile[1]);
#$suxwigb-> percent(99.9);
$suxwigb->clip(1);
if ( $iSelect_tr_Sumute->{_number_of_tries} > 0 ) {
$iSelect_tr_Sumute->{_TX_outbound} =
'itemp_'
. $iSelect_tr_Sumute->{_purpose}
. '_picks_'
. $iSelect_tr_Sumute->{_file_in} . '_'
. $iSelect_tr_Sumute->{_gather_type}
. $iSelect_tr_Sumute->{_gather_num};
lib/App/SeismicUnixGui/big_streams/iSunmo.pm view on Meta::CPAN
$suximage->cmap( quotemeta('hsv2') );
$suximage->first_x( quotemeta(1) ); #d2=1,f2=1
$suximage->dx( quotemeta(1) );
$suximage->loclip( quotemeta( $suximage_col_bar_min[1] ) );
$suximage->hiclip( quotemeta( $suximage_col_bar_max[1] ) );
$suximage->verbose( quotemeta($off) );
#$suximage -> dx_major_divisions(quotemeta($iSunmo->{_velocity_increment}*10));
#$suximage -> dy_minor_divisions(quotemeta($number_minor_time_divisions);
#$suximage -> dy_major_divisions(quotemeta($time_inc_major);
#$suximage -> percent4clip(95.0);
#$suximage -> first_tick_number_x($iSunmo->{_first_velocity});
# $suximage -> picks(\$iSunmo->{_Tvel_outbound});
#print("Writing picks to $iSunmo->{_Tvel_outbound}\n\n");
$suximage[1] = $suximage->Step();
=head2
DEFINE FLOW(S)
in interactive mode:
lib/App/SeismicUnixGui/big_streams/iSuvelan.pm view on Meta::CPAN
$suximage->verbose( quotemeta($off) );
#print ("d2num is ($iSuvelan->{_velocity_increment})\n\n");
$suximage->dx_major_divisions( quotemeta($dx_major_divisions) );
$suximage->dy_minor_divisions(
quotemeta($number_minor_time_divisions) );
$suximage->dy_major_divisions( quotemeta($time_inc_major) );
$suximage->tend_s( quotemeta( $iSuvelan->{_tmax_s} ) );
#$suximage -> percent4clip(quotemeta(95.0);
$suximage->first_tick_number_x(
quotemeta( $iSuvelan->{_first_velocity} ) );
$suximage->picks( $iSuvelan->{_Tvel_outbound} );
print("iSuvelan, calcNdisplay: Writing picks to $iSuvelan->{_Tvel_outbound}\n\n");
=head2 conditions
when number_of_tries is >=2
there should be a pre-exisiting digitized
lib/App/SeismicUnixGui/configs/plot/supswigb.config view on Meta::CPAN
labelfont = Helvetica
labelsize = 24
loclip|wclip =
n1 =
n2 =
nbpi = 72
npair =
num_minor_ticks_betw_time_ticks|n1tic =
num_minor_ticks_betw_distance_ticks|n2tic =
orientation|style = 'seismic'
percent|perc =
pick|mpicks =
plotfile =
shading|va =
title = 'nu'
tend_s|x1end =
ticwidth =
tstart_s|x1beg =
titlefont = Helvetica-Bold
titlecolor =
titlesize = 24
lib/App/SeismicUnixGui/configs/plot/supswigp.config view on Meta::CPAN
interp =
labelcolor =
labelfont =
loclip|wclip =
n1 =
n2 =
npair =
num_minor_ticks_betw_time_ticks|n1tic =
num_minor_ticks_betw_distance_ticks|n2tic =
orientation|style = 'seismic'
percent|perc =
pick|mpicks =
plotfile =
shading|va =
title = 'nu'
tend_s|x1end =
tstart_s|x1beg =
titlefont =
titlecolor =
tmpdir =
trace_inc|d2|dx|dp =
lib/App/SeismicUnixGui/configs/plot/suximage.config view on Meta::CPAN
lheight =
lwidth = 16
lx = 3
ly =
picks|mpicks = '/dev/tty'
n1 =
num_minor_ticks_betw_time_ticks|n1tic = 1
n2 =
num_minor_ticks_betw_distance_ticks|n2tic = 1
npair =
percent4clip|perc = 100.0
plotfile = 'plotfile.ps'
orientation|style = 'seismic'
title = 'suximage'
titlecolor = 'red'
titlefont = 'Rom22'
tmpdir =
tstart_s|x1beg =
tend_s|x1end =
units = 'unit'
verbose = 1
lib/App/SeismicUnixGui/configs/plot/suxwigb.config view on Meta::CPAN
interp =
labelcolor =
labelfont =
loclip|wclip =
n1 =
n2 =
npair =
num_minor_ticks_betw_time_ticks|n1tic =
num_minor_ticks_betw_distance_ticks|n2tic =
orientation|style = 'seismic'
percent|perc =
pick|mpicks =
plotfile =
shading|va =
title = 'nu'
tend_s|x1end =
tstart_s|x1beg =
titlefont =
titlecolor =
tmpdir =
trace_inc|d2|dx|dp =
lib/App/SeismicUnixGui/developer/Stripped/NMO_Vel_Stk/suvelan_nsel.su.main.velocity_analysis view on Meta::CPAN
achieve maximum amplitude of unity. The constant is just the inverse of the
total number of crosscorrelations included in the sum. The selection is
made using a parabolic approximation of the differential moveout and
imposing a threshold for those differential moveouts.
That threshold is the parameter tau in this program, which varies between 0
to 1. A value of tau=0, means conventional crosscorrelation sum is applied
implying that all crosscorrelations are included in the sum. In contrast,
a value of tau=1 (not recomended) means that only the crosscorrelation
formed by the trace pair involving the shortest and longest offset is
included in the sum. Intermediate values will produce percentages of the
crosscorrelations included in the sum that will be shown in the screen
before computing the velocity spectra. Typical values for tau are between
0.2 and 0.6, producing approximated percentages of crosscorrelations summed
between 60% and 20%. The higher the value of tau the lower the percentage
and higher the increase in the resolving power of velocity spectra.
Keeping the percentage of crosscorrelations included in the sum between 20%
and 60% will increase resolution and avoid the precense of artifacts in
the results. In data contaminated by random noise or statics distortions
is recomended to mantaing the percentage of crosscorrelations included in
the sum above 25%. After computing the velocity spectra one might want to
adjust the level and number of contours before velocity picking.
Credits: CWP: Valmore Celis, Sept 2002
Based on the original code: suvelan.c
Colorado School of Mines: Dave Hale c. 1989
References:
lib/App/SeismicUnixGui/developer/Stripped/NMO_Vel_Stk/suvelan_usel.su.main.velocity_analysis view on Meta::CPAN
moveout, thus increasing the resolving power in the velocity spectra
compared to that achieved by conventional methods. The selection is made
using a parabolic approximation of the differential moveout and imposing a
threshold for those differential moveouts.
That threshold is the parameter tau in this program, which varies between
0 to 1. A value of tau=0, means conventional crosscorrelation sum is
applied implying that all crosscorrelations are included in the sum. In
contrast, a value of tau=1 (not recomended) means that only the
crosscorrelation formed by the trace pair involving the shortest and longest
offset is included in the sum. Intermediate values will produce percentages
of the crosscorrelations included in the sum that will be shown in the
screen before computing the velocity spectra. Typical values for tau are
between 0.2 and 0.6, producing approximated percentages of crosscorrelations
summed between 60% and 20%. The higher the value of tau the lower the
percentage and higher the increase in the resolving power of velocity
spectra.
Keeping the percentage of crosscorrelations included in the sum between 20%
and 60% will increase resolution and avoid the precense of artifacts in the
results. In data contaminated by random noise or statics distortions is
recomended to mantaing the percentage of crosscorrelations included in the
sum above 25%. After computing the velocity spectra one might want to
adjust the level and number of contours before velocity picking.
Credits: CWP: Valmore Celis, Sept 2002
Based on the original code: suvelan_.c
Colorado School of Mines: Dave Hale c. 1989
lib/App/SeismicUnixGui/developer/Stripped/filter/succfilt.su.main.filters view on Meta::CPAN
SUCCFILT - FX domain Correlation Coefficient FILTER
sucff < stdin > stdout [optional parameters]
Optional parameters:
cch=1.0 Correlation coefficient high pass value
ccl=0.3 Correlation coefficient low pass value
key=ep ensemble identifier
padd=25 FFT padding in percentage
Notes:
This program uses "get_gather" and "put_gather" so requires that
the data be sorted into ensembles designated by "key", with the ntr
field set to the number of traces in each respective ensemble.
Example:
susort ep offset < data.su > datasorted.su
suputgthr dir=Data verbose=1 < datasorted.su
sugetgthr dir=Data verbose=1 > dataupdated.su
lib/App/SeismicUnixGui/developer/Stripped/filter/suphidecon.su.main.decon_shaping view on Meta::CPAN
SUPHIDECON - PHase Inversion Deconvolution
suphidecon < stdin > stdout
Required parameters:
none
Optional parameters:
... time range used for wavelet extraction:
tm=-0.1 Pre zero time (maximum phase component )
tp=+0.4 Post zero time (minimum phase component + multiples)
percpad=50 percentage padding for nt prior to cepstrum calculation
pnoise=0.001 Pre-withening (assumed noise to prevent division by zero)
Notes:
The wavelet is separated from the reflectivity and noise based on
their different 'smoothness' in the pseudo cepstrum domain.
The extracted wavelet also includes multiples.
The wavelet is reconstructed in frequency domain, end removed ",
from the trace. (Method by Lichman and Northwood, 1996.)
lib/App/SeismicUnixGui/developer/Stripped/header/suedit.su.main.headers view on Meta::CPAN
suedit diskfile (open for possible header modification if writable)
suedit <diskfile (open read only)
The following commands are recognized:
number read in that trace and print nonzero header words
<CR> go to trace one step away (step is initially -1)
+ read in next trace (step is set to +1)
- read in previous trace (step is set to -1)
dN advance N traces (step is set to N)
% print some percentiles of the trace data
r print some ranks (rank[j] = jth smallest datum)
p [n1 [n2]] tab plot sample n1 to n2 on current trace
g [tr1 tr2] ximage plot the trace [traces tr1 to tr2]
w [tr1 tr2] xwigb plot the trace [traces tr1 to tr2]
f [tr1 tr2] ximage plot the amplitude spectra of the trace
u [tr1 tr2] apply user pipeline to specified traces
! key=val change a value in a field (e.g. ! tracr=101)
? print help file
q quit
lib/App/SeismicUnixGui/developer/Stripped/header/suresstat.su.main.statics view on Meta::CPAN
niter=5 number of iterations
nshot=240 largest shot number (fldr=1 to nshot)
nr=335 largest receiver number (tracf=1 to nr)
nc=574 maximum number of cmp's (for array allocation)
sfold=96 maximum shot gather fold
rfold=96 maximum receiver gather fold
cfold=48 maximum cmp gather fold
sub=0 subtract super trace 1 from super trace 2 (=1)
sub=0 strongly biases static to a value of 0
mode=0 use global maximum in cross-correllation window
=1 choose the peak perc=percent smaller than the global max.
perc=10. percent of global max (used only for mode=1)
verbose=0 print diagnostic output (verbose=1)
Notes:
Estimates surface-consistent source and receiver statics, meaning that
there is one static correction value estimated for each shot and receiver
position.
The method employed here is based on the method of Ronen and Claerbout:
Geophysics 50, 2759-2767 (1985).
The output files are binary files containing the source and receiver
statics, as a function of shot number (trace header fldr) and
receiver station number (trace header tracf).
The code builds a supertrace1 and supertrace2, which are subsequently
cross-correllated. The program then picks the time lag associated with
the largest peak in the cross-correllation according to two possible
criteria set by the parameter "mode". If mode=0, the maximum of the
cross-correllation window is chosen. If mode=1, the program will pick
a peak which is up to perc=percent smaller than the global maximum, but
closer to zero lag than the global maximum. (Choosing mode=0 is
recommended.)
The geometry can be irregular: the program simply computes a static
correction for each shot record (fldr=1 to fldr=nshot), with any missing
shots being assigned a static of 0. A static correction for each
receiver station (tracf=1 to tracf=nr) is calculated, with missing
receivers again assigned a static of 0. ",
The ntracesces parameter must be equal to the number of prestack traces.
lib/App/SeismicUnixGui/developer/Stripped/header/suxedit.su.main.headers view on Meta::CPAN
suxedit diskfile (open for possible header modification if writable)
suxedit <diskfile (open read only)
The following commands are recognized:
number read in that trace and print nonzero header words
<CR> go to trace one step away (step is initially -1)
+ read in next trace (step is set to +1)
- read in previous trace (step is set to -1)
dN advance N traces (step is set to N)
% print some percentiles of the trace data
r print some ranks (rank[j] = jth smallest datum)
p [n1 [n2]] tab plot sample n1 to n2 on current trace
g [tr1 tr2] ["opts"] wiggle plot (graph) the trace
[traces tr1 to tr2]
f wiggle plot the Fourier transform of the trace
! key=val change a value in a field (e.g. ! tracr=101)
? print help file
q quit
NOTE: sample numbers are 1-based (first sample is 1).
lib/App/SeismicUnixGui/developer/Stripped/model/addrvl3d.par.main view on Meta::CPAN
=2 random vels are always higher than v_orig
cvel=2000 layer filled with constant velocity cvel
(overides vlsd,add,how params)
---->Smoothing parameters (0 = no smoothing)
r1=0.0 1st dimension operator length in samples
r2=0.0 2nd dimension operator length in samples
r3=0.0 3rd dimension operator length in samples
slowness=0 =1 smoothing on slowness; =0 smoothing on velocity
nrvl=n1/10 number of const velocity layers to add
pdv=10. percentage velocity deviation (max) from input model
Notes:
1. Smoothing radii usually fall in the range of [0,20].
2. Smoothing radii can be used to set aspect ratio of random velocity
anomalies in the new layer. For example (r1=5,r2=0,r3=0) will
result in vertical vel streaks that mimick vertical fracturing.
3. Smoothing on slowness works better to preserve traveltimes relative
to the unsmoothed case.
4. Default case is a random velocity (+/-30%) near surface layer whose
thickness is 20% of the total 2D model thickness.
lib/App/SeismicUnixGui/developer/Stripped/model/suremac2d.su.main.synthetics_waveforms_testpatterns view on Meta::CPAN
4. A 20 grid points wide border at the sides and the bottom of
the modeling grid is used for sponge boundary conditions
(default: iabso=1).
Source and receiver lines should be placed some (e.g. 10) grid
points away from the absorbing boundaries in order to reduce
reflections due to obliquely incident wavefronts.
5. Dominant frequency is about fmax/2 (sflag=2), absolute maximum
is delayed by 3/fmax from beginning of wavelet.
6. If opflag!=1 the source should be not a spike in space; the
parameter w determines at which distance (in grid points) from
the source's center the Gaussian weight decays to 10 percent
of its maximum. w=2 may be a reasonable choice; however, the
waveform will be distorted.
7. Horizontal and vertical receiver line sections are written to
separate files. Each file can hold more than one line.
8. Parameter vmaxu may be enlarged if the modeling run becomes
unstable. This happens if the largest eigenvalue of the modeling
operator L is larger than estimated from the largest velocity.
In particular if using the variable density acoustic wave
equation the eigenvalues depend also on the density and it is
impossible to estimated the largest eigenvalue analytically.
lib/App/SeismicUnixGui/developer/Stripped/model/suremel2dan.su.main.synthetics_waveforms_testpatterns view on Meta::CPAN
the modeling grid is used for sponge boundary conditions
(default: iabso=1).
Source and receiver lines should be placed some (e.g. 10) grid
points away from the absorbing boundaries in order to reduce
reflections due to obliquely incident wavefronts.
5. Dominant frequency is about fmax/2 (sflag=2), absolute maximum
is delayed by 3/fmax from beginning of wavelet.
6. If source is not single force (i.e. pressure or shear source)
it should be not a spike in space; the parameter w determines
at which distance (in grid points) from the source's center
the Gaussian weight decays to 10 percent of its maximum.
w=2 may be a reasonable choice; however, the waveform will be
distorted.
7. Horizontal and vertical receiver line sections are written to
separate files. Each file can hold more than one line.
8. Parameter vmaxu may be enlarged if the modeling run becomes
unstable. This happens if the largest eigenvalue of the modeling
operator L is larger than estimated from the largest velocity.
9. Bessel coefficients can be precomputed (prec=1) and stored on
disk to save CPU time when several shots need to be run.
In this case computation of Bessel coefficients can be skipped
lib/App/SeismicUnixGui/developer/Stripped/par/mrafxzwt.par.main view on Meta::CPAN
Required Parameters:
n1= size of first (fast) dimension
n2= size of second (slow) dimension
Optional Parameters:
p1= maximum integer such that 2^p1 <= n1
p2= maximum integer such that 2^p2 <= n2
order=6 order of Daubechies wavelet used (even, 4<=order<=20)
mralevel1=3 maximum multi-resolution analysis level in dimension 1
mralevel2=3 maximum multi-resolution analysis level in dimension 2
trunc=0.0 truncation level (percentage) of the reconstruction
verbose=0 =1 to print some useful information
reconfile= reconstructed data file to write
reconmrafile= reconstructed data file in MRA domain to write
dfile= difference between infile and reconfile to write
dmrafile= difference between mrafile and reconmrafile to write
dconly=0 =1 keep only dc component of MRA
verbose=0 =1 to print some useful information
if (n1 or n2 is not integer powers of 2) specify the following:
nc1=n1/2 center of trimmed image in the 1st dimension
nc2=n2/2 center of trimmed image in the 2nd dimension
lib/App/SeismicUnixGui/developer/Stripped/par/randvel3d.par.main view on Meta::CPAN
=2 random vels are always higher than v_orig
cvel=2000 layer filled with constant velocity cvel
(overides vlsd,add,how params)
---->Smoothing parameters (0 = no smoothing)
r1=0.0 1st dimension operator length in samples
r2=0.0 2nd dimension operator length in samples
r3=0.0 3rd dimension operator length in samples
slowness=0 =1 smoothing on slowness; =0 smoothing on velocity
nrvl=n1/10 number of const velocity layers to add
pdv=10. percentage velocity deviation (max) from input model
Notes:
1. Smoothing radii usually fall in the range of [0,20].
2. Smoothing radii can be used to set aspect ratio of random velocity
anomalies in the new layer. For example (r1=5,r2=0,r3=0) will
result in vertical vel streaks that mimick vertical fracturing.
3. Smoothing on slowness works better to preserve traveltimes relative
to the unsmoothed case.
4. Default case is a random velocity (+/-30%) near surface layer whose
thickness is 20% of the total 2D model thickness.
lib/App/SeismicUnixGui/developer/Stripped/plot/pscube.psplot.main view on Meta::CPAN
front= name of file containing front panel
side= name of file containing side panel
top= name of file containing top panel
faces=0 =1 to read faces from stdin (data format 2)
d1=1.0 sampling interval in 1st dimension
f1=0.0 first sample in 1st dimension
d2=1.0 sampling interval in 2nd dimension
f2=0.0 first sample in 2nd dimension
d3=1.0 sampling interval in 3rd dimension
f3=0.0 first sample in 3rd dimension
perc=100.0 percentile used to determine clip
clip=(perc percentile) clip used to determine bclip and wclip
bperc=perc percentile for determining black clip value
wperc=100.0-perc percentile for determining white clip value
bclip=clip data values outside of [bclip,wclip] are clipped
wclip=-clip data values outside of [bclip,wclip] are clipped
brgb=0.0,0.0,0.0 red, green, blue values corresponding to black
wrgb=1.0,1.0,1.0 red, green, blue values corresponding to white
bhls=0.0,0.0,0.0 hue, lightness, saturation corresponding to black
whls=0.0,1.0,0.0 hue, lightness, saturation corresponding to white
bps=12 bits per sample for color plots, either 12 or 24
d1s=1.0 factor by which to scale d1 before imaging
d2s=1.0 factor by which to scale d2 before imaging
d3s=1.0 factor by which to scale d3 before imaging
lib/App/SeismicUnixGui/developer/Stripped/plot/psimage.psplot.main view on Meta::CPAN
Required Parameters:
n1 number of samples in 1st (fast) dimension
Optional Parameters:
d1=1.0 sampling interval in 1st dimension
f1=0.0 first sample in 1st dimension
n2=all number of samples in 2nd (slow) dimension
d2=1.0 sampling interval in 2nd dimension
f2=0.0 first sample in 2nd dimension
perc=100.0 percentile used to determine clip
clip=(perc percentile) clip used to determine bclip and wclip
bperc=perc percentile for determining black clip value
wperc=100.0-perc percentile for determining white clip value
bclip=clip data values outside of [bclip,wclip] are clipped
wclip=-clip data values outside of [bclip,wclip] are clipped
bclip and wclip will be set to be inside
[lbeg,lend] if lbeg and/or lend are supplied
threecolor=1 supply 3 color values instead of only two,
using not only black and white, but f.e. red,
green and blue
brgb=0.0,0.0,0.0 red, green, blue values corresponding to black
grgb=1.0,1.0,1.0 red, green, blue values corresponding to grey
wrgb=1.0,1.0,1.0 red, green, blue values corresponding to white
lib/App/SeismicUnixGui/developer/Stripped/plot/psmovie.psplot.main view on Meta::CPAN
Required Parameters:
n1 number of samples in 1st (fast) dimension
Optional Parameters:
d1=1.0 sampling interval in 1st dimension
f1=0.0 first sample in 1st dimension
n2=all number of samples in 2nd (slow) dimension
d2=1.0 sampling interval in 2nd dimension
f2=0.0 first sample in 2nd dimension
perc=100.0 percentile used to determine clip
clip=(perc percentile) clip used to determine bclip and wclip
bperc=perc percentile for determining black clip value
wperc=100.0-perc percentile for determining white clip value
bclip=clip data values outside of [bclip,wclip] are clipped
wclip=-clip data values outside of [bclip,wclip] are clipped
d1s=1.0 factor by which to scale d1 before imaging
d2s=1.0 factor by which to scale d2 before imaging
verbose=1 =1 for info printed on stderr (0 for no info)
xbox=1.0 offset in inches of left side of axes box
ybox=1.5 offset in inches of bottom side of axes box
wbox=6.0 width in inches of axes box
hbox=8.0 height in inches of axes box
x1beg=x1min value at which axis 1 begins
lib/App/SeismicUnixGui/developer/Stripped/plot/pswigb.psplot.main view on Meta::CPAN
n1 number of samples in 1st (fast) dimension
Optional Parameters:
d1=1.0 sampling interval in 1st dimension
f1=0.0 first sample in 1st dimension
n2=all number of samples in 2nd (slow) dimension
d2=1.0 sampling interval in 2nd dimension
f2=0.0 first sample in 2nd dimension
x2=f2,f2+d2,... array of sampled values in 2nd dimension
bias=0.0 data value corresponding to location along axis 2
perc=100.0 percentile for determining clip
clip=(perc percentile) data values < bias+clip and > bias-clip are clipped
xcur=1.0 wiggle excursion in traces corresponding to clip
wt=1 =0 for no wiggle-trace; =1 for wiggle-trace
va=1 =0 for no variable-area; =1 for variable-area fill
=2 for variable area, solid/grey fill
SHADING: 2<= va <=5 va=2 lightgrey, va=5 black",
nbpi=72 number of bits per inch at which to rasterize
verbose=1 =1 for info printed on stderr (0 for no info)
xbox=1.5 offset in inches of left side of axes box
ybox=1.5 offset in inches of bottom side of axes box
wbox=6.0 width in inches of axes box
lib/App/SeismicUnixGui/developer/Stripped/plot/pswigp.psplot.main view on Meta::CPAN
n1 number of samples in 1st (fast) dimension
Optional Parameters:
d1=1.0 sampling interval in 1st dimension
f1=0.0 first sample in 1st dimension
n2=all number of samples in 2nd (slow) dimension
d2=1.0 sampling interval in 2nd dimension
f2=0.0 first sample in 2nd dimension
x2=f2,f2+d2,... array of sampled values in 2nd dimension
bias=0.0 data value corresponding to location along axis 2
perc=100.0 percentile for determining clip
clip=(perc percentile) data values < bias+clip and > bias-clip are clipped
xcur=1.0 wiggle excursion in traces corresponding to clip
fill=1 =0 for no fill;
>0 for pos. fill;
<0 for neg. fill
=2 for pos. fill solid, neg. fill grey
=-2for neg. fill solid, pos. fill grey
SHADING: 2<=abs(fill)<=5 2=lightgrey 5=black
linewidth=1.0 linewidth in points (0.0 for thinest visible line)
tracecolor=black color of traces; should contrast with background
backcolor=none color of background; none means no background
lib/App/SeismicUnixGui/developer/Stripped/plot/xpicker.xplot.main view on Meta::CPAN
Optional Parameters:
mpicks=pick_file name of output (input) pick file
d1=1.0 sampling interval in 1st dimension
f1=d1 first sample in 1st dimension
n2=all number of samples in 2nd (slow) dimension
d2=1.0 sampling interval in 2nd dimension
f2=d2 first sample in 2nd dimension
x2=f2,f2+d2,... array of sampled values in 2nd dimension
bias=0.0 data value corresponding to location along axis 2
perc=100.0 percentile for determining clip
clip=(perc percentile) data values < bias+clip and > bias-clip are clipped
xcur=1.0 wiggle excursion in traces corresponding to clip
wt=1 =0 for no wiggle-trace; =1 for wiggle-trace
va=1 =0 for no variable-area; =1 for variable-area fill
=2 for variable area, solid/grey fill
SHADING: 2<=va<=5 va=2 light grey, va=5 black
verbose=1 =1 for info printed on stderr (0 for no info)
xbox=50 x in pixels of upper left corner of window
ybox=50 y in pixels of upper left corner of window
wbox=550 width in pixels of window
hbox=700 height in pixels of window
lib/App/SeismicUnixGui/developer/Stripped/shapeNcut/suresstat.su.main.statics view on Meta::CPAN
niter=5 number of iterations
nshot=240 largest shot number (fldr=1 to nshot)
nr=335 largest receiver number (tracf=1 to nr)
nc=574 maximum number of cmp's (for array allocation)
sfold=96 maximum shot gather fold
rfold=96 maximum receiver gather fold
cfold=48 maximum cmp gather fold
sub=0 subtract super trace 1 from super trace 2 (=1)
sub=0 strongly biases static to a value of 0
mode=0 use global maximum in cross-correllation window
=1 choose the peak perc=percent smaller than the global max.
perc=10. percent of global max (used only for mode=1)
verbose=0 print diagnostic output (verbose=1)
Notes:
Estimates surface-consistent source and receiver statics, meaning that
there is one static correction value estimated for each shot and receiver
position.
The method employed here is based on the method of Ronen and Claerbout:
Geophysics 50, 2759-2767 (1985).
The output files are binary files containing the source and receiver
statics, as a function of shot number (trace header fldr) and
receiver station number (trace header tracf).
The code builds a supertrace1 and supertrace2, which are subsequently
cross-correllated. The program then picks the time lag associated with
the largest peak in the cross-correllation according to two possible
criteria set by the parameter "mode". If mode=0, the maximum of the
cross-correllation window is chosen. If mode=1, the program will pick
a peak which is up to perc=percent smaller than the global maximum, but
closer to zero lag than the global maximum. (Choosing mode=0 is
recommended.)
The geometry can be irregular: the program simply computes a static
correction for each shot record (fldr=1 to fldr=nshot), with any missing
shots being assigned a static of 0. A static correction for each
receiver station (tracf=1 to tracf=nr) is calculated, with missing
receivers again assigned a static of 0. ",
The ntracesces parameter must be equal to the number of prestack traces.
lib/App/SeismicUnixGui/developer/Stripped/transform/suradon.su.main.transforms view on Meta::CPAN
pmax=400 maximum moveout in ms on reference offset
dp=16 moveout increment in ms on reference offset
pmula=80 moveout in ms on reference offset where multiples begin
at maximum time
pmulb=200 moveout in ms on reference offset where multiples begin
at zero time
depthref=500. Reference depth for Foster/Mosher hyperbolic transform
nwin=1 number of windows to use through the mute zone
f1=60. High-end frequency before taper off
f2=80. High-end frequency
prewhite=0.1 Prewhitening factor in percent.
cdpkey=cdp name of header word for defining ensemble
offkey=offset name of header word with spatial information
nxmax=240 maximum number of input traces per ensemble
ltaper=7 taper (integer) for mute tapering function
Optimizing Parameters:
The following parameters are occasionally used to avoid spatial aliasing
problems on the linear tau-p transform. Not recommended for other
transforms...
ninterp=0 number of traces to interpolate between each input trace
lib/App/SeismicUnixGui/sunix/NMO_Vel_Stk/sunmo.pm view on Meta::CPAN
smute=1.5 samples with NMO stretch exceeding smute are zeroed
lmute=25 length (in samples) of linear ramp for stretch mute
sscale=1 =1 to divide output samples by NMO stretch factor
invert=0 =1 to perform (approximate) inverse NMO
upward=0 =1 to scan upward to find first sample to kill
voutfile= if set, interplolated velocity function v[cdp][t] is
output to named file.
Notes:
For constant-velocity NMO, specify only one vnmo=constant and omit tnmo.
NMO interpolation error is less than 1 percent for frequencies less than
60 percent of the Nyquist frequency.
Exact inverse NMO is impossible, particularly for early times at large
offsets and for frequencies near Nyquist with large interpolation errors.
The "offset" header field must be set.
Use suazimuth to set offset header field when sx,sy,gx,gy are all
nonzero.
For NMO with a velocity function of time only, specify the arrays
vnmo=v1,v2,... tnmo=t1,t2,...
lib/App/SeismicUnixGui/sunix/NMO_Vel_Stk/suvelan_nsel.pm view on Meta::CPAN
That threshold is the parameter tau in this program, which varies between 0
to 1. A value of tau=0, means conventional crosscorrelation sum is applied
implying that all crosscorrelations are included in the sum. In contrast,
a value of tau=1 (not recomended) means that only the crosscorrelation
formed by the trace pair involving the shortest and longest offset is
included in the sum. Intermediate values will produce percentages of the
crosscorrelations included in the sum that will be shown in the screen
before computing the velocity spectra. Typical values for tau are between
0.2 and 0.6, producing approximated percentages of crosscorrelations summed
between 60 0x0p+0nd 20%. The higher the value of tau the lower the percentage
and higher the increase in the resolving power of velocity spectra.
Keeping the percentage of crosscorrelations included in the sum between 20%
and 60% will increase resolution and avoid the precense of artifacts in
the results. In data contaminated by random noise or statics distortions
is recomended to mantaing the percentage of crosscorrelations included in
the sum above 25%. After computing the velocity spectra one might want to
adjust the level and number of contours before velocity picking.
Credits: CWP: Valmore Celis, Sept 2002
( run in 0.532 second using v1.01-cache-2.11-cpan-05162d3a2b1 )