def simulateddata(sdat, rdat, slo, sd, ss, par):
# single reflector
par['kxref'] = 50
par['lxref'] = par['nximg'] - par['kxref']
par['kzref'] = par['nzimg'] - par['nzpad'] / par['jzimg'] - 10
Flow(
'ref', None, '''
spike nsp=1 mag=1
n1=%(nximg)d d1=%(dximg)g o1=%(oximg)g k1=%(kxref)d l1=%(lxref)d
n2=%(nzimg)d d2=%(dzimg)g o2=%(ozimg)g k2=%(kzref)d |
smooth rect1=25 repeat=3 |
spray axis=2 n=1 o=0 d=1 |
put label1=x label2=y label3=z
''' % par)
Plot('ref',
'window | transp | smooth rect1=3 |' + fdmod.cgrey('pclip=100', par))
Result('ref', ['ref', ss], 'Overlay')
# source wavelet
par['frq'] = 5
par['kt'] = 50
Flow(
'wvl', None, '''
spike nsp=1 mag=1
n1=%(nt)d d1=%(dt)g o1=0 k1=%(kt)d |
ricker1 frequency=%(frq)g |
scale axis=123 |
fft1 |
window squeeze=n n1=%(nw)d min1=%(ow)g
''' % par)
# source data on the surface
for iexp in range(par['ns']):
etag = "-e%03d" % iexp
xsou = par['os'] + (par['fs'] + iexp * par['js']) * par['ds']
isou = (xsou - par['oximg']) / par['dximg']
Flow(
'spk' + etag, 'wvl', '''
pad beg2=%d n2out=%d |
put label1=w label2=x label3=y o2=%g d2=%g |
transp memsize=250 plane=12 |
transp memsize=250 plane=23
''' % (
isou,
par['nximg'],
par['oximg'],
par['dximg'],
))
Flow('spk', map(lambda x: 'spk-e%03d' % x, range(par['ns'])),
'cat space=n axis=4 ${SOURCES[1:%d]}' % par['ns'])
# datumed source data
zomig.Cdtone3(sdat, 'spk', sd, par)
Result(sdat, fdmod.fgrey('', par))
# datumed receiver data
spmig.modelPW3(rdat, slo, sdat, 'ref', par)
Result(rdat, fdmod.fgrey('', par))
def simulateddata(sdat,rdat,slo,sd,ss,par):
# single reflector
par['kxref']=50
par['lxref']=par['nximg']-par['kxref']
par['kzref']=par['nzimg']-par['nzpad']/par['jzimg']-10
Flow('ref',None,
'''
spike nsp=1 mag=1
n1=%(nximg)d d1=%(dximg)g o1=%(oximg)g k1=%(kxref)d l1=%(lxref)d
n2=%(nzimg)d d2=%(dzimg)g o2=%(ozimg)g k2=%(kzref)d |
smooth rect1=25 repeat=3 |
spray axis=2 n=1 o=0 d=1 |
put label1=x label2=y label3=z
''' % par )
Plot('ref','window | transp | smooth rect1=3 |'
+ fdmod.cgrey('pclip=100',par))
Result('ref',['ref',ss],'Overlay')
# source wavelet
par['frq']=5
par['kt']=50
Flow('wvl',None,
'''
spike nsp=1 mag=1
n1=%(nt)d d1=%(dt)g o1=0 k1=%(kt)d |
ricker1 frequency=%(frq)g |
scale axis=123 |
fft1 |
window squeeze=n n1=%(nw)d min1=%(ow)g
''' % par)
# source data on the surface
for iexp in range(par['ns']):
etag = "-e%03d" % iexp
xsou = par['os'] + (par['fs'] + iexp * par['js']) * par['ds']
isou = (xsou-par['oximg'])/par['dximg']
Flow('spk'+etag,'wvl',
'''
pad beg2=%d n2out=%d |
put label1=w label2=x label3=y o2=%g d2=%g |
transp memsize=250 plane=12 |
transp memsize=250 plane=23
''' % (isou,
par['nximg'],
par['oximg'],
par['dximg'],) )
Flow('spk',
map(lambda x: 'spk-e%03d' % x,range(par['ns'])),
'cat space=n axis=4 ${SOURCES[1:%d]}'%par['ns'])
# datumed source data
zomig.Cdtone3(sdat,'spk',sd,par)
Result(sdat,fdmod.fgrey('',par))
# datumed receiver data
spmig.modelPW3(rdat,slo,sdat,'ref',par)
Result(rdat,fdmod.fgrey('',par))
def originaldata(sdat,rdat,sd,par):
Flow( 'wav',None,'spike nsp=1 mag=1 k1=1 n1=%(nt)d d1=%(dt)g o1=%(ot)g' % par)
Result('wav','window n1=500 |' + fdmod.waveplot('',par))
sigsbee.getdata('data',par)
sigsbee.makeshots('sall','data',par)
Flow('swin',
'sall',
'''
window n3=%(ns)d f3=%(fs)d j3=%(js)d |
bandpass flo=2 fhi=10
''' % par)
for file in (['swin','sall']):
for iexp in range(par['ns']):
etag = "-e%03d" % iexp
Plot(file+etag,file,
'window n3=1 f3=%d | put o2=%g |' % (iexp,par['jos']+iexp*par['jds'])
+ fdmod.dgrey('pclip=99',par))
Result(file,map(lambda x: file+'-e%03d' % x,range(par['ns'])),'Movie')
encode.shot2grid('sdata-t','rdata-t','wav','swin',par)
Flow('S-sou-t','sdata-t','window | transp' %par)
Flow('S-rec-t','rdata-t','window | transp' %par)
encode.time2freq('S-sou-t','S-sou-w',par)
encode.time2freq('S-rec-t','S-rec-w',par)
# datum through water
spmig.datum3('S-dfs',
'S-dfr',
sd,
'S-sou-w',
'S-rec-w',par)
# window datumed data (all shots)
Flow(sdat,'S-dfs','window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)
Flow(rdat,'S-dfr','window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)
Result(sdat,fdmod.fgrey('pclip=99',par))
Result(rdat,fdmod.fgrey('pclip=99',par))
def originaldata(sdat, rdat, sd, par):
Flow('wav', None,
'spike nsp=1 mag=1 k1=1 n1=%(nt)d d1=%(dt)g o1=%(ot)g' % par)
Result('wav', 'window n1=500 |' + fdmod.waveplot('', par))
sigsbee.getdata('data', par)
sigsbee.makeshots('sall', 'data', par)
Flow(
'swin', 'sall', '''
window n3=%(ns)d f3=%(fs)d j3=%(js)d |
bandpass flo=2 fhi=10
''' % par)
for file in (['swin', 'sall']):
for iexp in range(par['ns']):
etag = "-e%03d" % iexp
Plot(
file + etag, file, 'window n3=1 f3=%d | put o2=%g |' %
(iexp, par['jos'] + iexp * par['jds']) +
fdmod.dgrey('pclip=99', par))
Result(file, map(lambda x: file + '-e%03d' % x, range(par['ns'])),
'Movie')
encode.shot2grid('sdata-t', 'rdata-t', 'wav', 'swin', par)
Flow('S-sou-t', 'sdata-t', 'window | transp' % par)
Flow('S-rec-t', 'rdata-t', 'window | transp' % par)
encode.time2freq('S-sou-t', 'S-sou-w', par)
encode.time2freq('S-rec-t', 'S-rec-w', par)
# datum through water
spmig.datum3('S-dfs', 'S-dfr', sd, 'S-sou-w', 'S-rec-w', par)
# window datumed data (all shots)
Flow(sdat, 'S-dfs', 'window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)
Flow(rdat, 'S-dfr', 'window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)
Result(sdat, fdmod.fgrey('pclip=99', par))
Result(rdat, fdmod.fgrey('pclip=99', par))
