SYNOPSIS
PERL PROGRAM NAME:
AUTHOR: Juan Lorenzo (Perl module only)
DATE:
DESCRIPTION:
Version:
USE
NOTES
Examples
SYNOPSIS
SEISMIC UNIX NOTES
UNIF2ANISO - generate a 2-D UNIFormly sampled profile of elastic
constants from a layered model.
unif2aniso < infile [Parameters]
Required Parameters:
none
Optional Parameters:
ninf=5 number of interfaces
nx=100 number of x samples (2nd dimension)
nz=100 number of z samples (1st dimension)
dx=10 x sampling interval
dz=10 z sampling interval
npmax=201 maximum number of points on interfaces
fx=0.0 first x sample
fz=0.0 first z sample
x0=0.0,0.0,..., distance x at which vp00 and vs00 are specified
z0=0.0,0.0,..., depth z at which vp00 and vs00 are specified
vp00=1500,2000,..., P-velocity at each x0,z0 (m/sec)
vs00=866,1155..., S-velocity at each x0,z0 (m/sec)
rho00=1000,1100,..., density at each x0,z0 (kg/m^3)
q00=110,120,130,.., attenuation Q, at each x0,z0 (kg/m^3)
eps00=0,0,0..., Thomsen or Sayers epsilon
delta00=0,0,0..., Thomsen or Sayers delta
gamma00=0,0,0..., Thomsen or Sayers gamma
dqdx=0.0,0.0,..., x-derivative of Q (d q/dx)
dqdz=0.0,0.0,..., z-derivative of Q (d q/dz)
drdx=0.0,0.0,..., x-derivative of density (d rho/dx)
drdz=0.0,0.0,..., z-derivative of density (d rho/dz)
dvpdx=0.0,0.0,..., x-derivative of P-velocity (dvp/dx)
dvpdz=0.0,0.0,..., z-derivative of P-velocity (dvs/dz)
dvsdx=0.0,0.0,..., x-derivative of S-velocity (dvs/dx)
dvsdz=0.0,0.0,..., z-derivative of S-velocity (dvs/dz)
dedx=0.0,0.0,..., x-derivative of epsilon (de/dx)
dedz=0.0,0.0,..., z-derivative of epsilon with depth z (de/dz)
dddx=0.0,0.0,..., x-derivative of delta (dd/dx)
dddz=0.0,0.0,..., z-derivative of delta (dd/dz)
dgdz=0.0,0.0,..., x-derivative of gamma (dg/dz)
dgdx=0.0,0.0,..., z-derivative of gamma (dg/dx)
phi00=0,0,..., rotation angle(s) in each layer
...output filenames
c11_file=c11_file output filename for c11 values
c13_file=c13_file output filename for c13 values
c15_file=c15_file output filename for c15 values
c33_file=c33_file output filename for c33 values
c35_file=c35_file output filename for c35 values
c44_file=c44_file output filename for c44 values
c55_file=c55_file output filename for c55 values
c66_file=c66_file output filename for c66 values
rho_file=rho_file output filename for density values
q_file=q_file output filename for Q values
paramtype=1 =1 Thomsen parameters, =0 Sayers parameters(see below)
method=linear for linear interpolation of interface
=mono for monotonic cubic interpolation of interface
=akima for Akima's cubic interpolation of interface
=spline for cubic spline interpolation of interface
tfile= =testfilename if set, a sample input dataset is
output to "testfilename".
prevents completion of demos -JL 7.14.21
Notes:
The input file is an ASCII file containing x z values representing a
piecewise continuous velocity model with a flat surface on top.
The surface and each successive boundary between media is represented
by a list of selected x z pairs written column form. The first and
last x values must be the same for all boundaries. Use the entry
1.0 -99999 to separate the entries for successive boundaries. No
boundary may cross another. Note that the choice of the method of
interpolation may cause boundaries to cross that do not appear to
cross in the input data file.
The number of interfaces is specified by the parameter "ninf". This
number does not include the top surface of the model. The input data
format is the same as a CSHOT model file with all comments removed.
The algorithm works by transforming the P-wavespeed , S-wavespeed,
density and the Thomsen or Sayers parameters epsilon, delta, and gamma
into elastic stiffness coefficients. Furthermore, the user can specify
rotations, phi, to the elasticity tensor in each layer.
Common ranges of Thomsen parameters are
epsilon: 0.0 -> 0.5
delta: -0.2 -> 0.4
gamma: 0.0 -> 0.4
If only P-wave, S-wave velocities and density is given as input,
the model is, by definition, isotropic.
If files containing Thomsen/Sayers parameters are given, the model
will be assumed to have VTI symmetry.
Example using test input file generating feature:
unif2aniso tfile=testfilename produces a 5 interface demonstration model
unif2aniso < testfilename
ximage < c11_file n1=100 n2=100
ximage < c13_file n1=100 n2=100
ximage < c15_file n1=100 n2=100
ximage < c33_file n1=100 n2=100
ximage < c35_file n1=100 n2=100
ximage < c44_file n1=100 n2=100
ximage < c55_file n1=100 n2=100
ximage < c66_file n1=100 n2=100
ximage < rho_file n1=100 n2=100
ximage < q_file n1=100 n2=100
Credits:
CWP: John Stockwell, April 2005.
CWP: based on program unif2 by Zhenyue Liu, 1994 CHANGES and their DATES
Import packages
instantiation of packages
Encapsulated hash of private variables
sub Step
collects switches and assembles bash instructions by adding the program name
unif2aniso writes anisotropy files locally wherever the executable is run
sub note
collects switches and assembles bash instructions by adding the program name
sub clear
sub c11_file
sub c13_file
sub c15_file
sub c33_file
sub c35_file
sub c44_file
sub c55_file
sub c66_file
sub dddx
sub dddz
sub dedx
sub dedz
sub delta00
sub dgdx
sub dgdz
sub dqdx
sub dqdz
sub drdx
sub drdz
sub dvpdx
sub dvpdz
sub dvsdx
sub dvsdz
sub dx
sub dz
sub eps00
sub fx
sub fz
sub gamma00
sub method
sub n1
sub ninf
sub npmax
sub nx
sub nz
sub paramtype
sub phi00
sub q00
sub q_file
sub rho00
sub rho_file
sub tfile
sub vp00
sub vs00
sub x0
sub z0
sub get_max_index
max index = number of input variables -1