SYNOPSIS
PERL PROGRAM NAME:
AUTHOR: Juan Lorenzo (Perl module only)
DATE:
DESCRIPTION:
Version:
USE
NOTES
Examples
SYNOPSIS
SEISMIC UNIX NOTES SUGASSMAN - Model reflectivity change with rock/fluid properties
sugassman [optional parameters] > data.su
Optional parameters:
nt=500 number of time samples
ntr=200 number of traces
dt=0.004 time sampling interval in seconds
mode=0 model isolated gassmann refl coefficient
=1 embed gassmann RC in random RC series
=2 R0 parameter sensitivity output
p=0.15 parameter sensitivity test range (if mode=2)
.... Environment variables ...
temp=140 Temperature in degrees C
pres=20 Pressure in megaPascals
.... Caprock variables ....
v1=37900 caprock P-wave speed (m/s)
r1=44300 caprock mass density (g/cc)
.... Reservoir fluid variables ....
g=0.56 Gas specific gravity 0.56 (methane)-1.8 (condensate)
api=50 Gas specific gravity 10 (heavy)-50 (ultra light)
s=35 Brine salinity in ppm/(1000 000
so=.7 Oil saturation (0-1)
sg=.2 Gas saturation (0-1)
.... Reservoir rock frame variables ....
kmin=37900 Bulk modulus (MPa) of rock frame mineral(s) [default=quartz]
mumin=44300 Shear modulus (MPa) of rock frame mineral(s) [default=quartz]
rmin=2.67 Mass density (g/cc) of rock frame mineral(s) [default=quartz]
phi=.24 Rock frame porosity (0-1)
a=1 Fitting parameters: Mdry/Mmineral ~ 1/(a + b phi^c)
b=15 ... where M is P-wave modulus and defaults are for
c=1 ... Glenn sandstone [see Liner (2nd ed, table 26.2)]
h=20 Reservoir thickness (m)
Notes:
Creates a reflection coefficient series based on Gassmann
theory of velocity and density for porous elastic media
Credits: UHouston: Chris Liner 9/23/2009
trace header fields set: User's notes (Juan Lorenzo) untested
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
sub note
collects switches and assembles bash instructions by adding the program name
sub clear
sub a
sub api
sub b
sub c
sub dt
sub g
sub h
sub kmin
sub mode
sub mumin
sub nt
sub ntr
sub p
sub phi
sub pres
sub r1
sub rmin
sub s
sub sg
sub so
sub temp
sub v1
sub get_max_index
max index = number of input variables -1