App-SeismicUnixGui
view release on metacpan or search on metacpan
lib/App/SeismicUnixGui/developer/Stripped/header/suresstat.su.main.statics view on Meta::CPAN
suresstat <stdin [optional parameters]
Required parameters:
ssol= output file source statics
rsol= output file receiver statics
ntraces= number of traces in input data set (must be correct!)
Optional parameters:
ntpick=50 maximum static shift (samples)
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.
lib/App/SeismicUnixGui/developer/Stripped/header/suresstat.su.main.statics view on Meta::CPAN
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.
The ntpick parameter sets the maximum allowable shift desired (in
samples NOT time).
The niter parameter sets the number of iterations desired.
The nshot parameter must be equal to the maximum fldr number in
the data. Note that this number might be different from the actual
number of shot records in the data (i.e., the maximum ep number).
For getting the correct maximum fldr number, you may use the surange
command.
The nr parameter must be equal to the largest number of receivers
per shot in the whole data.
The nc parameter must be equal to the number of prestack traces in
the data.
The sfold parameter must be equal to the nr parameter.
lib/App/SeismicUnixGui/developer/Stripped/model/sufdmod1.su.main.synthetics_waveforms_testpatterns view on Meta::CPAN
}
else for (iz=0; iz<nz; iz++) rd[iz]=2500;
for (iz=0; iz<nz-1; iz++) rd1_5[iz]=(rd[iz]+rd[iz+1])/2;
rd1_5[nz-1]=rd[nz-1];
/* time step computation
vmax=0;
for (iz=0; iz<nz; iz++) if (rv[iz]>vmax) vmax=rv[iz];if (verbose) warn( "vmax= %f ", vmax);
dt=dz/1.414/vmax/2; if (verbose) warn( "time step dt= %f ", dt);
/* maximum number of iterations
if (nt==0) nt=1+tmax/dt;
if (verbose) warn( "number of time steps nt= %d ", nt);
if (NINT( (float) nt/((float)td))+1>SU_NFLTS) err("too many time steps. Increase td ?");
/* source parameter computation
alpha=2*9.8696*freq*freq;
/* time shift to get a t0 centered source
if ((styp==0) || (styp == 2)) epst0=fabs(source (0, styp, dt, dz, 0, alpha) / 1e4);
lib/App/SeismicUnixGui/developer/Stripped/model/susynlvfti.su.main.synthetics_waveforms_testpatterns view on Meta::CPAN
For transversely isotropic media:
angxs=0.0 angle of symmetry axis with the vertical (degrees)
define the media using either
a=1.0 corresponding to the ratio of elastic coef.(c1111/c3333)
f=0.4 corresponding to the ratio of elastic coef. (c1133/c3333)
l=0.3 corresponding to the ratio of elastic coef. (c1313/c3333)
Alternately use Tompson\'s parameters:
delta=0 Thomsen's 1986 defined parameter
epsilon=0 Thomsen's 1986 defined parameter
ntries=40 number of iterations in Snell's law and offset searches
epsx=.001 lateral offset tolerance
epst=.0001 reflection time tolerance
nitmax=12 max number of iterations in travel time integrations
Notes:
Offsets are signed - may be positive or negative. Receiver locations
are computed by adding the signed offset to the source location.
Specify either midpoint sampling or shotpoint sampling, but not both.
If neither is specified, the default is the midpoint sampling above.
More than one ref (reflector) may be specified. When obliquity factors
lib/App/SeismicUnixGui/developer/Stripped/shapeNcut/suresstat.su.main.statics view on Meta::CPAN
suresstat <stdin [optional parameters]
Required parameters:
ssol= output file source statics
rsol= output file receiver statics
ntraces= number of traces in input data set (must be correct!)
Optional parameters:
ntpick=50 maximum static shift (samples)
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.
lib/App/SeismicUnixGui/developer/Stripped/shapeNcut/suresstat.su.main.statics view on Meta::CPAN
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.
The ntpick parameter sets the maximum allowable shift desired (in
samples NOT time).
The niter parameter sets the number of iterations desired.
The nshot parameter must be equal to the maximum fldr number in
the data. Note that this number might be different from the actual
number of shot records in the data (i.e., the maximum ep number).
For getting the correct maximum fldr number, you may use the surange
command.
The nr parameter must be equal to the largest number of receivers
per shot in the whole data.
The nc parameter must be equal to the number of prestack traces in
the data.
The sfold parameter must be equal to the nr parameter.
lib/App/SeismicUnixGui/developer/Stripped/statsMath/suharlan.su.main.noise view on Meta::CPAN
SUHARLAN - signal-noise separation by the invertible linear
transformation method of Harlan, 1984
suharlan <infile >outfile [optional parameters]
Required Parameters:
<none>
Optional Parameters:
FLAGS:
niter=1 number of requested iterations
anenv=1 =1 for positive analytic envelopes
=0 for no analytic envelopes (not recommended)
scl=0 =1 to scale output traces (not recommended)
plot=3 =0 for no plots. =1 for 1-D plots only
=2 for 2-D plots only. =3 for all plots
norm=1 =0 not to normalize reliability values
verbose=1 =0 not to print processing information
rgt=2 =1 for uniform random generator
=2 for gaussian random generator
sts=1 =0 for no smoothing (not recommended)
lib/App/SeismicUnixGui/developer/Stripped/statsMath/suharlan.su.main.noise view on Meta::CPAN
np=100 number of p-values for taup transform
prewhite=0.01 prewhitening value (suggested between 0.1 and 0,01)
offref=2000 reference offset for p-values (m)
depthref=500 reference depth for Foster/Mosher taup (if gopt=4)
pmula=pmax1 maximum p-value preserved in the data (ms)
pmulb=pmax1 minimum p-value muted on the data (ms)
ninterp=0 number of traces to interpolate in input data
Extraction Parameters:
nintlh=50 number of intervals (bins) in histograms
sditer=5 number of steepest descent iterations to compute ps
c=0.04 maximum noise allowed in a sample of signal(%)
rel1=0.5 reliability value for first pass of the extraction
rel2=0.75 reliability value for second pass of the extraction
Smoothing Parameters: ",
r1=10 number of points for damped lsq vertical smoothing
r2=2 number of points for damped lsq horizontal smoothing
Output Files:
lib/App/SeismicUnixGui/geopsy/dinver.pm view on Meta::CPAN
Optimization options:
-target <TARGET> Set targets from file TARGET. It can be a .target
or a .dinver file. A .dinver file contains both
parameters and targets. Provide it to both options
'-param' and '-target'.
-param <PARAM> Set parameters from file PARAM. It can be a .param
or a .dinver file. A .dinver file contains both
parameters and targets. Provide it to both options
'-param' and '-target'.
-itmax <ITMAX> Number of iterations started by option -run
(default=50).
-ns0 <NS0> Number of initial models (default=50).
-ns <NS> Number of models generated (default=50).
-nr <NR> Number of best cells (default=50).
-seed <SEED> Set random seed to SEED. This option is for debug
only.
-o <REPORT> Ouput report (default='run.report'). To suppress
output to .report file, set option '-o' to "".
-f Force overwrite if output report already exists.
-resume If the output report file already exists, it is
imported into the parameter space before starting
the inversion. Better to set NS0 to zero in this
case. This way, it is possible to continue an
existing inversion adding new iterations.
Importance sampling options:
-dof <DOF> Set degrees of freedom for conversion from misfit
to a posteriori probability function.
-param <PARAM> Set parameters from file PARAM. It can be a .param
or a .dinver file. A .dinver file contains both
parameters and targets. Provide it to both options
'-param' and '-target'.
-ns <NS> Number of models generated (default=5000).
-seed <SEED> Set random seed to SEED. This option is for debug
lib/App/SeismicUnixGui/geopsy/dinver.pm view on Meta::CPAN
$Dinver->{_note} = $Dinver->{_note} . ' -param ' . $Dinver->{_param};
$Dinver->{_Step} = $Dinver->{_Step} . ' -param ' . $Dinver->{_param};
}
####
=pod
=head2 Subroutine itmax
Number of iterations for inversion to run
0 iterations = Pure Monte Carlo
=cut
sub itmax {
my ( $sub, $itmax ) = @_;
$Dinver->{_itmax} = $itmax if defined($itmax);
$Dinver->{_note} = $Dinver->{_note} . ' -itmax ' . $Dinver->{_itmax};
$Dinver->{_Step} = $Dinver->{_Step} . ' -itmax ' . $Dinver->{_itmax};
}
lib/App/SeismicUnixGui/geopsy/dinver.pm view on Meta::CPAN
$Dinver->{_note} = $Dinver->{_note} . $Dinver->{_force};
$Dinver->{_Step} = $Dinver->{_Step} . $Dinver->{_force};
}
=pod
=head2 Subroutine resume
If output file already exists, it is improted before starting inversion
Safer to set ns0 -> 0 for this case
Adds new iterations to existing inversion
=cut
sub resume {
my ( $sub, $resume ) = @_;
$Dinver->{_resume} = $resume if defined($resume);
$Dinver->{_note} = $Dinver->{_note} . ' -resume ' . $Dinver->{_resume};
$Dinver->{_Step} = $Dinver->{_Step} . ' -resume ' . $Dinver->{_resume};
}
lib/App/SeismicUnixGui/sunix/header/suresstat.pm view on Meta::CPAN
rsol= output file receiver statics
ntraces= number of traces in input data set (must be correct!)
Optional parameters:
ntpick=50 maximum static shift (samples)
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
lib/App/SeismicUnixGui/sunix/header/suresstat.pm view on Meta::CPAN
receivers again assigned a static of 0. ",
The ntracesces parameter must be equal to the number of prestack traces.
The ntpick parameter sets the maximum allowable shift desired (in
samples NOT time).
The niter parameter sets the number of iterations desired.
The nshot parameter must be equal to the maximum fldr number in
the data. Note that this number might be different from the actual
number of shot records in the data (i.e., the maximum ep number).
For getting the correct maximum fldr number, you may use the surange
command.
lib/App/SeismicUnixGui/sunix/header/suresstat_old.pm view on Meta::CPAN
Notes:
Package name is the same as the file name
Moose is a package that allows an object-oriented
syntax to organizing your programs
=========================================================
MODIFIED BY: Martial Morrison
DATE: July 14 2016
DESCRIPTION: Updated for Seismic Unix v44
Version 1.1
Added new subroutines: fn, imax, input_file (fn), max_fldr_cdp_tracf (imax), iterations (niter), max_sample_shift (ntpick)
-These parameters were added to suresstat in SU v44
Removed outdated subroutines: cfold, rfold, ntraces, nr, nc, nshot, number_of_traces
-These parameters are no longer part of suresstat
=cut
=head2 Notes from Seismic Unix
SURESSTAT - Surface consistent source and receiver statics calculation
suresstat fn= [optional parameters]
Required parameters:
fn= seismic file
ssol= output file source statics
rsol= output file receiver statics
Optional parameters:
ntpick=50 maximum static shift (samples)
niter=5 number of iterations
imax=100000 largest shot (fldr),reciver(tracf) or cmp(cdp) number
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
lib/App/SeismicUnixGui/sunix/header/suresstat_old.pm view on Meta::CPAN
=cut
my $suresstat = {
_fn => '',
_input_file => '',
_imax => '',
_max_fldr_cdp_tracf => '',
_ntpick => '',
_max_sample_shift => '',
_niter => '',
_iterations => '',
_subtract => '',
_mode => '',
_perc => '',
_ssol => '',
_source_statics_file_output => '',
_rsol => '',
_receiver_statics_file_output => '',
_note => '',
_Step => '',
_verbose => ''
lib/App/SeismicUnixGui/sunix/header/suresstat_old.pm view on Meta::CPAN
=cut
sub clear {
$suresstat->{_fn} = '';
$suresstat->{_input_file} = '';
$suresstat->{_imax} = '';
$suresstat->{_max_fldr_cdp_tracf} = '';
$suresstat->{_ntpick} = '';
$suresstat->{_max_sample_shift} = '';
$suresstat->{_niter} = '';
$suresstat->{_iterations} = '';
$suresstat->{_subtract} = '';
$suresstat->{_mode} = '';
$suresstat->{_perc} = '';
$suresstat->{_ssol} = '';
$suresstat->{_source_statics_file_output} = '';
$suresstat->{_rsol} = '';
$suresstat->{_receiver_statics_file_output} = '';
$suresstat->{_note} = '';
$suresstat->{_Step} = '';
$suresstat->{_verbose} = '';
lib/App/SeismicUnixGui/sunix/header/suresstat_old.pm view on Meta::CPAN
if ($niter) {
$suresstat->{_niter} = $niter;
$suresstat->{_Step} =
$suresstat->{_Step} . ' niter=' . $suresstat->{_niter};
$suresstat->{_note} =
$suresstat->{_note} . ' niter=' . $suresstat->{_niter};
}
}
=head2 subroutine iterations
=cut
sub iterations {
my ( $variable, $iterations ) = @_;
if ($iterations) {
$suresstat->{_iterations} = $iterations;
$suresstat->{_Step} =
$suresstat->{_Step} . ' niter=' . $suresstat->{_iterations};
$suresstat->{_note} =
$suresstat->{_note} . ' iterations=' . $suresstat->{_iterations};
}
}
=head2 subroutine note
=cut
sub note {
my ( $variable, $note ) = @_;
$suresstat->{_note} = $note;
lib/App/SeismicUnixGui/sunix/model/susynlvfti.pm view on Meta::CPAN
f=0.4 corresponding to the ratio of elastic coef. (c1133/c3333)
l=0.3 corresponding to the ratio of elastic coef. (c1313/c3333)
Alternately use Tompson\'s parameters:
delta=0 Thomsen's 1986 defined parameter
epsilon=0 Thomsen's 1986 defined parameter
ntries=40 number of iterations in Snell's law and offset searches
epsx=.001 lateral offset tolerance
epst=.0001 reflection time tolerance
nitmax=12 max number of iterations in travel time integrations
Notes:
Offsets are signed - may be positive or negative. Receiver locations
are computed by adding the signed offset to the source location.
( run in 1.020 second using v1.01-cache-2.11-cpan-96521ef73a4 )