bamp

NAME
SYNOPSIS
DESCRIPTION
FILES
SEE ALSO
BUGS
COPYRIGHT
AUTHOR

NAME

bamp - BSU measure inelastic amplitude decay as a function of frequency

SYNOPSIS

bamp [ -h | infile | emin | emax | fmin | fmax | delf | bw | tmax ]

DESCRIPTION

Basic Seismic Utilities (BSU) measures inelastic amplitude decay (corrected for spherical divergence). Assumes that the data are from a VSP which has not had any gain processing beyond recorded gain recovery. Program creates an output file, bamp.his, for use in a joint inversion for the two inelastic wave equation coefficients (stiffness and damping) under the Kelvin constitutive model. Requires plplot shared library if postscript plots to be produced. PLPLOT is a GNU function/subroutine library which can generate PostScript plots (http://plplot.sourceforge.net).

Options

-h

Online help giving details on command line arguments

infile

Input file name

emin

Minimum geophone elevation to include in the analysis zone.

emax

Maximum geophone elevation to include in the analysis zone.

fmin

Minimum frequency for the scan.

fmax

Maximum frequency for the scan.

delf

Frequency increment for the scan. One should not try to exceed the resolution implied by the recording aperture. Thus, delf should not be less than 1/(npts*fsamin), where npts is the number of samples in the trace, fsamin is the sample interval. Or put another way, be sure to record enough data to permit adaquate resolution in the frequency domain.

bw

Bandwidth of the filters. See comment above on resolution.

tmax

Maximum time to scan. Note: If the data have been aligned on the first arrivals, tmax will mark the time after the first arrival. See programs bpic and bshf.

NOTE:
If invoked with no options, will prompt user for input parameters.

EXAMPLE:
bamp bshftwav.seg 840.45 843.60 6.0 100.0 2.0 2.0 .5

Input file bshftwav.seg is processed for the subsurface interval between elevations 840.45 m (bottom) and 843.6 m (top). The minimum frequency scanned will be 6 Hz, the maximum frequency scanned will be 100 hz. The frequencies will be scanned every 2 Hz (thus center frequencies scanned are 6, 8, 10, ... to 100 Hz). The bandwidth of each filter is 2.0 hz and the input data window extends to 0.5 seconds.

FILES

bamp.ps

Summary Postscript plot of amplitude decay (nepers/m) as a function of frequency.

bampqc.ps

QC Postscript plot of amplitude decay as a function of range for each frequency.

bamp.his

Text file of frequency (hz), amplitude decay (nepers/m), standard deviation of decay.

bampdb.his

Text file of frequency (hz), amplitude decay (decibels/m), standard deviation of decay.

bampxxxx.lst

Listing file for bamp. The characters, xxxx, will be the first 4 char of the input file name (for example, would be bshf if the last process was bshf, to align the data on first arrivals).

standard output

list of frequency, decay (dB/m should be <0), decay (nepers/m, should be >0)

SEE ALSO

bhelp(1), bvas(1), bpic(1), bshf(1)

BUGS

No known bugs. In cases of noisy data, or a thin zone of analysis, one may obtain decay values suggesting an increase in amplitude with propagation distance (dB/m>0). The code fixes such cases to zero decay in the history file. Increasing the zone thickness (emax-emin) or limiting the scan of frequencies to those with the best signal to noise will reduce the number of zero fixed values. When examining bampqc.ps, good results correspond to negative decay values (sloping down with increasing frequency). Results which slope upward with increasing range are in conflict with the basic model, and are zeroed in the bamp.his file as mentioned above.

COPYRIGHT

Copyright © 2024 by Paul Michaels

This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

AUTHOR

P. Michaels, PE. <paulmichaels@.boisestate.edu>