Gpprofile

Synopsis

gpprofile [OPTIONS] [FILE]

Description

Prints profile for specified parameters, from a layered model given through stdin or FILE

Format for layered models:

 Line 1    <number of layers n including half-space for first model>
 Line 2    <thickness (m)> <Vp (m/s)> <Vs (m/s)> <Density (kg/m3)>[ <Qp> <Qs>]
 ....
 Line n+1    0 <Vp (m/s)> <Vs (m/s)> <Density (kg/m3)>[ <Qp> <Qs>]
 Line n+3  <number of layers including half-space for second model, if wished>
 ....

Quality factors are not mandatory. Any number of models can be given as input.

Gpprofile options

• -vp Export Vp profiles
• -vs Export Vs profiles (default)
• -rho Export density profiles
• -nu Export Poisson's ratio profiles
• -imp Export impedance profiles
• -minmax Export minimum and maximum selected profiles from input models
• -d, -depth <depth> Maximum depth for minimum and maximum profiles (default=100)
• -n <count> Number of samples for minimum and maximum profiles (default=100)
• -average-at <depth> Returns the average profile over <depth> meters (travel time average)

The travel time average velocity at depth $D$, $V_{s,D}$, is computed as:

$V_{s,D}={\frac {1}{D}}({\frac {D-\sum _{i=1}^{j}D_{i}}{V_{j+1}}}+\sum _{i=1}^{j}{\frac {H_{i}}{V_{i}}})$

with $D_{i}$ the depth of the layer $i$, $V_{i}$, the shear wave velocity of the layer $i$ and with $j$ so that $\sum _{i=1}^{j}D_{i}\leq D$ and $\sum _{i=1}^{j+1}D_{i}>D$

Generic options

• -help <SECTION> Show help about options. SECTION may be empty or: all, html, latex, examples, generic, gpprofile
• -version Show version information
• -app-version Show short version information
• -reportbug Start bug report dialog, information about bug is passed through stdin. This option is used internally to report bugs if option -nobugreport is not specified.
• -reportint Start bug report dialog, information about interruption is passed through stdin. This option is used internally to report interruptions if option -nobugreport is not specified.
• -nobugreport Do not generate bug report in case of error

Examples

The function gpprofile is particulary useful to compute the values of $V_{S,30}$, crucial for earthquake engineering:

gpprofile -average-at 30 file

It can even be computed for all the models estimated by dinver as realistic (here with a misfit lower than 0.1) during an inversion process:

gpdcreport file.report -m 0.1 | gpprofile -average-at 30 > out

Then, the distribution of this parameter, written in the file 'out', can be studied (average, standard deviation...).

Copyright

Copyright (c) 2002-2008 Marc Wathelet. All rights reserved. See http://www.geopsy.org for more information.