?? 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()
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -