sconstruct

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

from rsfproj import *

Flow('spike',None,
     'spike n1=500 n2=50 k1=100,300 nsp=2 o2=1 d2=1 | ricker1 frequency=5')
Flow('offset','spike','window n1=1 | noise type=n rep=y seed=2005 range=3')
Flow('offset2','offset','transp | headersort head=$SOURCE')
Flow('velocity','spike','window n2=1 | math output="3+0.5*x1" ')

Flow('data','spike offset2 velocity',
     'inmo offset=${SOURCES[1]} velocity=${SOURCES[2]} half=n')
Plot('reg','data',
     '''
     wiggle label2="Trace Number" unit2= title=Regular
     yreverse=y transp=y poly=y
     ''')

Plot('ireg','data offset2',
     '''
     wiggle xpos=${SOURCES[1]} xmax=3 xmin=-3
     label2=Offset title=Irregular
     yreverse=y transp=y poly=y
     ''')

Result('dwiggle','reg ireg','SideBySideAniso')

End()