sconstruct

国外免费地震资料处理软件包

from rsfproj import *

# Fetch the velocity model from the data server
###############################################
vstr = 'sigsbee2a_stratigraphy.sgy'
Fetch(vstr,'sigsbee')
Flow('zvstr tzvstr ./shead ./bshead',vstr,
     '''
     segyread
     tape=$SOURCE
     tfile=${TARGETS[1]}
     hfile=${TARGETS[2]}
     bfile=${TARGETS[3]}
     ''',stdin=0)

Flow('vel','zvstr',
     '''
     put d1=0.025 o2=10.025 d2=0.025 |
     window f1=200 |
     scale dscale=0.001
     ''')
Plot('vel',
     '''
     grey wanttitle=n color=j allpos=y 
     label1="Depth (kft)" label2="Lateral (kft)"
     screenratio=0.3125 screenht=4 labelsz=4
     ''')

# Compute traveltimes by solving the eikonal equation
#####################################################

eikonal = 'eikonal yshot=40 zshot=5 br1=5 br2=5'

Flow('eiko','vel',eikonal)
Plot('eiko',
     '''
     contour nc=30 dc=0.2 plotcol=7
     wantaxis=n wanttitle=n
     screenratio=0.3125 screenht=4
     ''')

Result('eiko','vel eiko','Overlay')

samples = range(2,50)

errs = []
spaces = []

for sample in samples:

    # Subsample velocity and traveltime
    ###################################

    # Solve the eikonal equation
    ############################

    # Compute error
    ###############

    pass

# Concatenate results
#####################

# Plot error versus spacing
###########################
    
End()