Marmousi2 model

OBC div data

The divergence operator was applied to the multicomponent OBC Pluto dataset. These files are obc_div_v_1.segy and obc_div_v_2.segy. Taking the divergence separates out the acoustic component of the data.

The file marmousi2 div SConstruct contains a list of rules that tell Madagascar to gather the div data files, append the header and produce plots of the data. This script is reproduced in table 7 and a plot of shot 50 is shown in figure

from rsf.proj import *
# Fetch Files from repository
Fetch("obc_div_v_1.segy","marm2")
Fetch("obc_div_v_2.segy","marm2")

# Convert Files to RSF and update header
Flow('obc_div_v_1','obc_div_v_1.segy', '''segyread tape=$SOURCE |
        put n2=1381 n3=320 o1=0 o2=0 o3=3000
        d2=12.32 d3=25 label1=Depth Z label2=Distance X label3=Shot-Cord
        unit1=s unit2=m unit2=m''')

Flow('obc_div_v_2','obc_div_v_2.segy', '''segyread tape=$SOURCE | 
        put n2=1381 n3=160 o1=0 o2=0 o3=11025
        d2=12.32 d3=25 label1=Depth Z label2=Distance X label2=Shot-Cord
        unit1=s unit2=m unit3=m''')
# Concatinate Datasets
Flow('div',['obc_div_v_1','obc_div_v_2'],'cat ${SOURCES[0:2]} axis=3',stdin=0)

# Plot Data
Result('movie','div','''window $SOURCE  
               min3=4250 max3=4250 n3=1 |
               grey color=I gainpanel=a 
               title=OBC Div Shot 50''')

End()

Table 7. SCons script generating images of the Marmousi2 Vx data

movie
Figure 4. Marmousi2 shot 50 of div data.

Marmousi2 model