Marmousi Model

Model

The Marmousi model contains 158 horizontally layered horizons. A series of normal faults and resulting tilted blocks complicates the model towards its center. The model sits under approximately 32 m of water and is 9.2 km in length and 3 km in depth.

The velocity model found in the Madagascar repository, marmvel.hh can easily be displayed. This grid contains 751 data points in the Z direction and 2301 data points in the x direction. Table 2 displays the proper header configuration.

n1=751	d1=.004	o1=0	label= Depth	unit1=km
n2=2301	d2=.004	o2=0	label2=Position	unit2=km

Table 2. Header information for Marmousi velocity models

The script found at marmousi/model/SConstruct was written to obtain the Marmousi model datasets, append the headers as necessary and display the data. This file is presented in Table 3 while the velocity model image is displayed in Figure 1.

from rsf.proj import *
# Fetch Files from repository
raw=['marmvel.hh','marmsmooth.HH']
for file in raw:
    Fetch(file,"marm")
    if file is 'marmvel.hh':
        d=.004
        fileOut='marmvel'
        t='Velocity Model'
    if file is 'marmsmooth.HH':
        d=.024
        fileOut='marmsmooth'
        t='Smoothed Velocity Model'
# Convert Files to RSF and update headers
    Flow(fileOut,file,'''dd form=native | 
        scale rscale=.001 | put
        label1=Depth label2=Position unit1=km unit2=km
        d1=%f d2=%f''' % (d,d))
# Plotting Section
    Result(fileOut,'''window $SOURCE  | 
        grey color=I gainpanel=a allpos=y scalebar=y 
        title=%s barlabel=Velocity barunit=km/s screenratio=.326
        screenht=3 wheretitle=t labelsz=4 titlesz=6 ''' % t)
End()

Table 3. SConstruct script generating the Marmousi velocity model images

Typing Command 1 within the Marmousimodel directory runs the script.

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

marmvel
Figure 1. Velocity model

Marmousi Model