Amoco model

Model

This model is a 2D subset of a 3D model, however the model does not vary perpendicularly to this line. The velocity model is 22 km across and 4 km in depth.

The velmodel.hh file did not need to be updated appreciably in this example. However, the appropriate header settings are found in table 2. Datums were spread every 12.5 meters to produce a 22km by 4 km grid.

	n1=321	d1=0.0125	o1=0	label1=Depth	unit1=km
	n2=1761	d2=0.0125	o2=0	label2=Position	unit2=km

Table 2. Amoco unsmoothed velocity model header information

A python SConstruct script that fetches the data sets, appends the header slightly and plots the velocity model can be found in Table 3. An image of the velocity profile is found in Figure 1.

from rsfproj import *
import fdmod

# Fetch Files from repository
Fetch("velsmooth.HH","Amoco")
Fetch("velmodel.hh","Amoco")
# Convert Files to RSF
Flow('velsmooth','velsmooth.HH', '''dd form=native | put
    unit1=km unit2=km label1=Depth Z label2=Position title=Velocity Model''')

Flow('velmodel','velmodel.hh', '''dd form=native | scale rscale=.001 | put
    unit1=km unit2=km label1=Depth Z label2=Position title=Velocity Model  
    ''')

Result('velmodel','''window $SOURCE  max1=4 | grey  scalebar=y
                     color=I gainpanel=a allpos=y screenratio=.18 screenht=2
                     scalebar=y barlabel=km/s title=Amoco Velocity Model 
                     labelsz=4 titlesz=6 wheretitle=t''')
End()

Table 3. Scons script that generates RSF formatted Amoco velocity model

Typing Command 1 within the amocomodel directory runs the script.

$$\texttt{bash-3.1\$ scons view}$$ (1)

velmodel
Figure 1. Amoco velocity model.

Amoco model