bamp - BSU measure inelastic amplitude decay as a function of frequency
bamp [ -h | infile | emin | emax | fmin | fmax | delf | bw | tmax ]
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.
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)
bhelp(1), bvas(1), bpic(1), bshf(1)
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 © 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.
P. Michaels, PE. <[email protected]>