App-SeismicUnixGui
view release on metacpan or search on metacpan
lib/App/SeismicUnixGui/developer/Stripped/datum/sudatumk2dr.su.main.datuming view on Meta::CPAN
fields tr.sx and tr.gx. Offset is computed automatically.
Author: Trino Salinas, 05/01/96, Colorado School of Mines
This code is based on sukzmig2d.c written by Zhenyue Liu, 03/01/95.
Subroutines from Dave Hale's modeling library were adapted in
this code to define topography using cubic splines.
This code implements a Kirchhoff extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is a far field approximation of the Berryhill's original
formula (Berryhill, 1979). This equation is the result of a stationary
phase analysis to get an analog asymptotic expansion for the two-and-one
half dimensional extrapolation formula (Bleistein, 1984).
The extrapolation formula permits the downward continuation of upgoing
waves and upward continuation of downgoing waves. For upward conti-
nuation of upgoing waves and downward continuation of downgoing waves,
the conjugate transpose of the equation is used (Bevc, 1993).
References :
Berryhill, J.R., 1979, Wave equation datuming: Geophysics,
44, 1329--1344.
_______________, 1984, Wave equation datuming before stack
lib/App/SeismicUnixGui/developer/Stripped/datum/sudatumk2ds.su.main.datuming view on Meta::CPAN
fields tr.sx and tr.gx. Offset is computed automatically.
Author: Trino Salinas, 05/01/96, Colorado School of Mines
This code is based on sukzmig2d.c written by Zhenyue Liu, 03/01/95.
Subroutines from Dave Hale's modeling library were adapted in
this code to define topography using cubic splines.
This code implements a Kirchhoff extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is a far field approximation of the Berryhill's original
formula (Berryhill, 1979). This equation is the result of a stationary
phase analysis to get an analog asymptotic expansion for the two-and-one
half dimensional extrapolation formula (Bleistein, 1984).
The extrapolation formula permits the downward continuation of upgoing
waves and upward continuation of downgoing waves. For upward conti-
nuation of upgoing waves and downward continuation of downgoing waves,
the conjugate transpose of the equation is used (Bevc, 1993).
References :
Berryhill, J.R., 1979, Wave equation datuming: Geophysics,
44, 1329--1344.
_______________, 1984, Wave equation datuming before stack
lib/App/SeismicUnixGui/developer/Stripped/datum/sukdmdcr.su.main.datuming view on Meta::CPAN
SUKDMDCR - 2.5D datuming of receivers for prestack, common source
data using constant-background data mapping formula.
(See selfdoc for specific survey requirements.)
sukdmdcr infile= outfile= [parameters]
Required file parameters:
infile=stdin file for input seismic traces
outfile=stdout file for output
ttfile= file for input traveltime tables
Required parameters describing the traveltime tables:
lib/App/SeismicUnixGui/developer/Stripped/datum/sukdmdcr.su.main.datuming view on Meta::CPAN
input file. The input file recfile should be a single column ascii file
with the depth of the recording surface at every surface location (first
source to last offset), with spacing equal to dxgo.
The same holds for the datuming surface, using datsurf, zdat, and datfile.
Assumptions and limitations:
1. This code implements a 2.5D extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is an adaptation of Bleistein & Jaramillo's 2.5D data
mapping formula for receiver extrapolation. This is the result of a
stationary phase analysis of the data mapping equation in the case of
a constant source location (shot gather).
Credits:
Authors: Steven D. Sheaffer (CWP), 11/97
References: Sheaffer, S., 1999, "2.5D Downward Continuation of the Seismic
lib/App/SeismicUnixGui/developer/Stripped/datum/sukdmdcs.su.main.datuming view on Meta::CPAN
SUKDMDCS - 2.5D datuming of sources for prestack common receiver
data, using constant-background data-mapping formula.
(See selfdoc for specific survey geometry requirements.)
sukdmdcs infile= outfile= [parameters]
Required parameters:
infile=stdin file for input seismic traces
outfile=stdout file for output
ttfile=file for input traveltime tables
Required parameters describing the traveltime tables:
fzt= first depth sample in traveltime table
lib/App/SeismicUnixGui/developer/Stripped/datum/sukdmdcs.su.main.datuming view on Meta::CPAN
input file. The input file recfile should be a single column ascii file
with the depth of the recording surface at every surface location (first
source to last offset), with spacing equal to dxgo.
The same holds for the datuming surface, using datsurf, zdat, and datfile.
Assumptions and limitations:
1. This code implements a 2.5D extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is an adaptation of Bleistein & Jaramillo's 2.5D data
mapping formula for receiver extrapolation. This is the result of a
stationary phase analysis of the data mapping equation in the case of
a constant input receiver location (receiver gather).
Credits:
Authors: Steven D. Sheaffer (CWP), 11/97
References: Sheaffer, S., 1999, "2.5D Downward Continuation of the Seismic
lib/App/SeismicUnixGui/developer/Stripped/header/surandhw.su.main.headers view on Meta::CPAN
Optional parameters:
key=tstat header key word to set
a=0 =1 flag to add original value to final key
noise=gauss noise probability distribution
=flat for uniform; default Gaussian
seed=from_clock random number seed (integer)
min=0 minimum random number
max=1 maximum radnom number
NOTES:
The value of header word key is computed using the formula:
val(key) = a * val(key) + rand
Example:
surandhw <indata key=tstat a=0 min=0 max=10 > outdata
lib/App/SeismicUnixGui/developer/Stripped/header/sushw.su.main.headers view on Meta::CPAN
b=0,... increment(s) within group
c=0,... group increment(s)
d=0,... trace number shift(s)
j=ULONG_MAX,ULONG_MAX,... number of elements in group
Notes:
Fields that are getparred must have the same number of entries as key
words being set. Any field that is not getparred is set to the default
value(s) above. Explicitly setting j=0 will set j to ULONG_MAX.
The value of each header word key is computed using the formula:
i = itr + d
val(key) = a + b * (i % j) + c * (int(i / j))
where itr is the trace number (first trace has itr=0, NOT 1)
Examples:
1. set every dt field to 4ms
sushw <indata key=dt a=4000 |...
2. set the sx field of the first 32 traces to 6400, the second 32 traces
to 6300, decrementing by -100 for each 32 trace groups
...| sushw key=sx a=6400 c=-100 j=32 |...
lib/App/SeismicUnixGui/developer/Stripped/inversion/suinvco3d.3D.Suinvco3d view on Meta::CPAN
yt1=0.0 y-coordinate of last input tables
fmax=0.25/dt Maximum frequency set for operator antialiasing
ang=180 Maximum dip angle allowed in the image
apet=45 aperture open angle for summation
alias=0 =1 to set the anti-aliasing filter
verbose=1 =1 to print some useful information
Notes:
The information needed in the computation of the weighting factor
in Kirchhoff inversion formula includes traveltime, amplitude,
and Beylkin determinant at each grid point for each source/receiver
point. For a 3-D nonzero common-offset inversion, the Beylkin
determinant is computed from a 3x3 matrix with each element
containing a sum of quantities from the source and the receiver.
The nine elements in the Beylkin matrix for each source/receiver
can be determined by eight quantities. These quantities can be
computed by the dynamic ray tracer. Tables of traveltime, amplitude,
and Beylkin matrix elements from each source and receiver are
pre-computed and stored in files.
lib/App/SeismicUnixGui/developer/Stripped/inversion/suinvco3d.3D.Suinvco3d view on Meta::CPAN
interpolation is used to determine the weighting factor at each
output grid point, and weighted diffraction summation is then
applied. For each midpoint, the traveltimes and weight factors are
calculated in the horizontal range of (xm-nxb*dx-z*tan(apet),
xm+nxb*dx+z*tan(apet)).
Offsets are signed - may be positive or negative. ",
This algorithm is based on the inversion formulas in chaper 5 of
_Mathematics of Multimensional Seismic Migration, Imaging and Inversion_
(Springer-Verlag, 2000), by Bleistein, N., Cohen, J.K.
and Stockwell, Jr., J. W.
lib/App/SeismicUnixGui/developer/Stripped/inversion/suinvzco3d.su.main.migration_inversion view on Meta::CPAN
nyo=101 number of output traces in y-direction
fzo=0.0 z-coordinate of first point in output trace
dzo=15.0 vertical spacing of output trace
nzo=101 number of points in output trace ",
fmax=0.25/dt Maximum frequency set for operator antialiasing
ang=180 Maximum dip angle allowed in the image
verbose=1 =1 to print some useful information
Notes:
This algorithm is based on formula (50) in Geophysics Vol. 51,
1552-1558, by Cohen, J., Hagin, F., and Bleistein, N.
Traveltime and amplitude are calculated by ray tracing.
Interpolation is used to calculate traveltime and amplitude. ",
For each midpoint, traveltime and amplitude are calculated in the
horizontal range of (xm-nxb*dx, xm+nxb*dx). Velocity is changed by
linear interpolation in two upper trianglar corners whose width is
nxc*dx and height is nzc*dz. ",
Eikonal equation will fail to solve if there is a polar turned ray.
lib/App/SeismicUnixGui/developer/Stripped/model/normray.tri.main view on Meta::CPAN
MODIFIED: Boyi Ou, Colorado School of Mines, 4/14/95
Notes:
This code can shoot rays from specified interface by users, normally you
need to use gbmodel2 to generate interface parameters for this code, both
code have a parameter named nrays, it should be same. If you just want to
shoot rays from one specified location, you need to specify xs1,zs1,
otherwise, leave them alone. If you want to shoot rays from surface, you need
to define surface equal to 1. The rays from one location will be
approximately symmetric with direction Normal_direction - ashift.(if nangle is
odd, it is symmetric, even, almost symmetric. The formula for the first take
off angle is: angle=normal_direction-nangle/2*dangle-ashift. If you only want to
see caustics, you specify caustic=1, if you want to see rays which does not
reach surface, you specify nonsurface=1.
/
lib/App/SeismicUnixGui/developer/Stripped/model/susynvxz.su.main.synthetics_waveforms_testpatterns view on Meta::CPAN
fxm=0.0 first midpoint (m)
fpeak=0.2/dt peak frequency of symmetric Ricker wavelet (Hz)
ref="1:1,2;4,2" reflector(s): "amplitude:x1,z1;x2,z2;x3,z3;...
smooth=0 =1 for smooth (piecewise cubic spline) reflectors
ls=0 =1 for line source; default is point source
tmin=10.0*dt minimum time of interest (sec)
ndpfz=5 number of diffractors per Fresnel zone
verbose=0 =1 to print some useful information
Notes:
This algorithm is based on formula (58) in Geo. Pros. 34, 686-703,
by N. Bleistein.
Offsets are signed - may be positive or negative. ",
Traveltime and amplitude are calculated by finite differences which
is done only in part of midpoints; in the skiped midpoint, interpolation
is used to calculate traveltime and amplitude. ",
More than one ref (reflector) may be specified.
Note that reflectors are encoded as quoted strings, with an optional
reflector amplitude: preceding the x,z coordinates of each reflector.
lib/App/SeismicUnixGui/developer/Stripped/model/susynvxzcs.su.main.synthetics_waveforms_testpatterns view on Meta::CPAN
ref="1:1,2;4,2" reflector(s): "amplitude:x1,z1;x2,z2;x3,z3;...
smooth=0 =1 for smooth (piecewise cubic spline) reflectors
ls=0 =1 for line source; =0 for point source
tmin=10.0*dt minimum time of interest (sec)
ndpfz=5 number of diffractors per Fresnel zone
cable=1 roll reciever spread with shot
=0 static reciever spread
verbose=0 =1 to print some useful information
Notes:
This algorithm is based on formula (58) in Geo. Pros. 34, 686-703,
by N. Bleistein.
Traveltime and amplitude are calculated by finite difference which
is done only in one of every NXD receivers; in skipped receivers,
interpolation is used to calculate traveltime and amplitude. ",
For each receiver, traveltime and amplitude are calculated in the
horizontal range of (xg-nxb*dx, xg+nxb*dx). Velocity is changed by
constant extropolation in two upper trianglar corners whose width is
nxc*dx and height is nzc*dz.
lib/App/SeismicUnixGui/developer/Stripped/statsMath/suocext.su.main.interp_extrap view on Meta::CPAN
The tdmo and vdmo arrays specify a velocity function of time that is
used to implement a first-order correction for depth-variable velocity.
The times in tdmo must be monotonically increasing.
For each offset, the minimum time at which a non-zero sample exists is
used to determine a mute time. Output samples for times earlier than this",
mute time will be zeroed. Computation time may be significantly reduced
if the input traces are zeroed (muted) for early times at large offsets.
A term for better amplitude reconstruction was added to Hale's formulation.
Credits: Carlos E. Theodoro (modification of Hale's SUDMOFK program)
Technical Reference:
C. Theodoro & K. Larner, 1998
Extrapolation of seismic data to small offsets (CWP-276).
Dip-Moveout Processing - SEG Course Notes
Dave Hale, 1988
lib/App/SeismicUnixGui/script/L_SU.pl view on Meta::CPAN
Moves SaveAs to L_SU Menu and removes Save button
V 0.3.2 has 4 flow panels
V 0.3.3 has dragNdrop deactivated to stabilize version
V 0.3.4 has classifies sunix programs using tabbed notebooks Sept. 12, 2018
V0.3.7 removed all ticks from strings in GUIS using control module
From now on users can write words with gaps and commas and L_SU will accept these
value and formulate the correct Seismic Unix sytnax.
V 0.3.8 Standardized format with PerlTidy, tidyviewer .perltidyrc Aug., 2019
V 0.3.9 Introduce Moose attributes to record real-time GUI history
V 0.4.5 Include PDL packages to handle interactive modeling and reading fortran-generated
files
V 0.5.0 new color_listbox class handles occupancy and vacancies among the listboxes March 2021
lib/App/SeismicUnixGui/sunix/datum/sudatumk2dr.pm view on Meta::CPAN
fields tr.sx and tr.gx. Offset is computed automatically.
Author: Trino Salinas, 05/01/96, Colorado School of Mines
This code is based on sukzmig2d.c written by Zhenyue Liu, 03/01/95.
Subroutines from Dave Hale's modeling library were adapted in
this code to define topography using cubic splines.
This code implements a Kirchhoff extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is a far field approximation of the Berryhill's original
formula (Berryhill, 1979). This equation is the result of a stationary
phase analysis to get an analog asymptotic expansion for the two-and-one
half dimensional extrapolation formula (Bleistein, 1984).
The extrapolation formula permits the downward continuation of upgoing
waves and upward continuation of downgoing waves. For upward conti-
nuation of upgoing waves and downward continuation of downgoing waves,
the conjugate transpose of the equation is used (Bevc, 1993).
References :
Berryhill, J.R., 1979, Wave equation datuming: Geophysics,
44, 1329--1344.
_______________, 1984, Wave equation datuming before stack
lib/App/SeismicUnixGui/sunix/datum/sudatumk2ds.pm view on Meta::CPAN
This code is based on sukzmig2d.c written by Zhenyue Liu, 03/01/95.
Subroutines from Dave Hale's modeling library were adapted in
this code to define topography using cubic splines.
This code implements a Kirchhoff extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is a far field approximation of the Berryhill's original
formula (Berryhill, 1979). This equation is the result of a stationary
phase analysis to get an analog asymptotic expansion for the two-and-one
half dimensional extrapolation formula (Bleistein, 1984).
The extrapolation formula permits the downward continuation of upgoing
waves and upward continuation of downgoing waves. For upward conti-
nuation of upgoing waves and downward continuation of downgoing waves,
the conjugate transpose of the equation is used (Bevc, 1993).
References :
lib/App/SeismicUnixGui/sunix/datum/sukdmdcr.pm view on Meta::CPAN
=head3 NOTES
=head4 Examples
=head2 SYNOPSIS
=head3 SEISMIC UNIX NOTES
SUKDMDCR - 2.5D datuming of receivers for prestack, common source
data using constant-background data mapping formula.
(See selfdoc for specific survey requirements.)
sukdmdcr infile= outfile= [parameters]
Required file parameters:
lib/App/SeismicUnixGui/sunix/datum/sukdmdcr.pm view on Meta::CPAN
Assumptions and limitations:
1. This code implements a 2.5D extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is an adaptation of Bleistein & Jaramillo's 2.5D data
mapping formula for receiver extrapolation. This is the result of a
stationary phase analysis of the data mapping equation in the case of
a constant source location (shot gather).
Credits:
lib/App/SeismicUnixGui/sunix/datum/sukdmdcs.pm view on Meta::CPAN
fields tr.sx and tr.gx. Offset is computed automatically.
Author: Trino Salinas, 05/01/96, Colorado School of Mines
This code is based on sukzmig2d.c written by Zhenyue Liu, 03/01/95.
Subroutines from Dave Hale's modeling library were adapted in
this code to define topography using cubic splines.
This code implements a Kirchhoff extraplolation operator that allows to
transfer data from one reference surface to another. The formula used in
this application is a far field approximation of the Berryhill's original
formula (Berryhill, 1979). This equation is the result of a stationary
phase analysis to get an analog asymptotic expansion for the two-and-one
half dimensional extrapolation formula (Bleistein, 1984).
The extrapolation formula permits the downward continuation of upgoing
waves and upward continuation of downgoing waves. For upward conti-
nuation of upgoing waves and downward continuation of downgoing waves,
the conjugate transpose of the equation is used (Bevc, 1993).
References :
Berryhill, J.R., 1979, Wave equation datuming: Geophysics,
44, 1329--1344.
_______________, 1984, Wave equation datuming before stack
lib/App/SeismicUnixGui/sunix/header/surandhw.pm view on Meta::CPAN
seed=from_clock random number seed (integer)
min=0 minimum random number
max=1 maximum radnom number
NOTES:
The value of header word key is computed using the formula:
val(key) = a * val(key) + rand
Example:
surandhw <indata key=tstat a=0 min=0 max=10 > outdata
lib/App/SeismicUnixGui/sunix/header/sushw.pm view on Meta::CPAN
b=0,... increment(s) within group
c=0,... group increment(s)
d=0,... trace number shift(s)
j=ULONG_MAX,ULONG_MAX,... number of elements in group
Notes:
Fields that are getparred must have the same number of entries as key
words being set. Any field that is not getparred is set to the default
value(s) above. Explicitly setting j=0 will set j to ULONG_MAX.
The value of each header word key is computed using the formula:
i = itr + d
val(key) = a + b * (i % j) + c * (int(i / j))
where itr is the trace number (first trace has itr=0, NOT 1)
Examples:
1. set every dt field to 4ms
sushw <indata key=dt a=4000 |...
2. set the sx field of the first 32 traces to 6400, the second 32 traces
to 6300, decrementing by -100 for each 32 trace groups
...| sushw key=sx a=6400 c=-100 j=32 |...
lib/App/SeismicUnixGui/sunix/inversion/suinvco3d.pm view on Meta::CPAN
verbose=1 =1 to print some useful information
Notes:
The information needed in the computation of the weighting factor
in Kirchhoff inversion formula includes traveltime, amplitude,
and Beylkin determinant at each grid point for each source/receiver
point. For a 3-D nonzero common-offset inversion, the Beylkin
determinant is computed from a 3x3 matrix with each element
containing a sum of quantities from the source and the receiver.
The nine elements in the Beylkin matrix for each source/receiver
lib/App/SeismicUnixGui/sunix/inversion/suinvco3d.pm view on Meta::CPAN
Offsets are signed - may be positive or negative. ",
This algorithm is based on the inversion formulas in chaper 5 of
_Mathematics of Multimensional Seismic Migration, Imaging and Inversion_
(Springer-Verlag, 2000), by Bleistein, N., Cohen, J.K.
and Stockwell, Jr., J. W.
=head2 User's notes (Juan Lorenzo)
lib/App/SeismicUnixGui/sunix/inversion/suinvzco3d.pm view on Meta::CPAN
ang=180 Maximum dip angle allowed in the image
verbose=1 =1 to print some useful information
Notes:
This algorithm is based on formula (50) in Geophysics Vol. 51,
1552-1558, by Cohen, J., Hagin, F., and Bleistein, N.
Traveltime and amplitude are calculated by ray tracing.
Interpolation is used to calculate traveltime and amplitude. ",
For each midpoint, traveltime and amplitude are calculated in the
lib/App/SeismicUnixGui/sunix/model/normray.pm view on Meta::CPAN
code have a parameter named nrays, it should be same. If you just want to
shoot rays from one specified location, you need to specify xs1,zs1,
otherwise, leave them alone. If you want to shoot rays from surface, you need
to define surface equal to 1. The rays from one location will be
approximately symmetric with direction Normal_direction - ashift.(if nangle is
odd, it is symmetric, even, almost symmetric. The formula for the first take
off angle is: angle=normal_direction-nangle/2*dangle-ashift. If you only want to
see caustics, you specify caustic=1, if you want to see rays which does not
reach surface, you specify nonsurface=1.
/
=head2 User's notes (Juan Lorenzo)
lib/App/SeismicUnixGui/sunix/model/susynvxz.pm view on Meta::CPAN
tmin=10.0*dt minimum time of interest (sec)
ndpfz=5 number of diffractors per Fresnel zone
verbose=0 =1 to print some useful information
Notes:
This algorithm is based on formula (58) in Geo. Pros. 34, 686-703,
by N. Bleistein.
Offsets are signed - may be positive or negative. ",
Traveltime and amplitude are calculated by finite differences which
is done only in part of midpoints; in the skiped midpoint, interpolation
lib/App/SeismicUnixGui/sunix/model/susynvxzcs.pm view on Meta::CPAN
cable=1 roll reciever spread with shot
=0 static reciever spread
verbose=0 =1 to print some useful information
Notes:
This algorithm is based on formula (58) in Geo. Pros. 34, 686-703,
by N. Bleistein.
Traveltime and amplitude are calculated by finite difference which
is done only in one of every NXD receivers; in skipped receivers,
interpolation is used to calculate traveltime and amplitude. ",
( run in 0.442 second using v1.01-cache-2.11-cpan-3cd7ad12f66 )