def new(par):
fdmod.point3('ssold0', par['ox'] + 3 * par['nx'] / 10 * par['dx'],
par['oz'], 1, par)
fdmod.point3('ssold1', par['ox'] + 5 * par['nx'] / 10 * par['dx'],
par['oz'], 1, par)
fdmod.point3('ssold2', par['ox'] + 6 * par['nx'] / 10 * par['dx'],
par['oz'], 1, par)
Flow('ssold', ['ssold0', 'ssold1', 'ssold2'],
'cat axis=2 space=n ${SOURCES[0:3]}',
stdin=0)
Flow('wavold', 'wav', 'window')
fdmod.lmodel('doold', 'woold', 'ddold', 'wdold', 'wavold', 'velo', 'refl',
'ssold', 'rr', 'jsnap=100 nbz=100 nbx=100 ', par)
Result('woold', fdmod.wgrey('', par))
Result('wdold', fdmod.wgrey('', par))
fdmod.point('ssnew0', par['ox'] + 3 * par['nx'] / 10 * par['dx'],
par['oz'], par)
fdmod.point('ssnew1', par['ox'] + 5 * par['nx'] / 10 * par['dx'],
par['oz'], par)
fdmod.point('ssnew2', par['ox'] + 6 * par['nx'] / 10 * par['dx'],
par['oz'], par)
Flow('ssnew', ['ssnew0', 'ssnew1', 'ssnew2'],
'cat axis=2 space=n ${SOURCES[0:3]}',
stdin=0)
Flow('wavnew', 'wav', 'window squeeze=n')
fdmod.lwefd1('donew', 'wonew', 'ddnew', 'wdnew', 'wavnew', 'velo', 'dens',
'refl', 'ssnew', 'rr', 'jsnap=100 nb=100 ', par)
Result('wonew', fdmod.wgrey('', par))
Result('wdnew', fdmod.wgrey('', par))
def run(par):
# experiments
fdmod.horizontal('rr', 0, par)
Plot('rr', 'window j2=10|' + fdmod.rrplot('', par))
Plot('ss', 'window |' + fdmod.ssplot('', par))
Plot('sx', 'window |' + fdmod.ssplot('plotcol=5', par))
# wavelet
fdmod.wavelet('wav_', par['frq'], par)
Flow('wav', 'wav_', 'transp')
Result('wav', 'window n2=200 |' + fdmod.waveplot('', par))
# velocity
Flow(
'vbck', None, '''
math n1=%(nz)d o1=%(oz)g d1=%(dz)g output="2000" |
spray axis=2 n=%(nx)d o=%(ox)g d=%(dx)g
''' % par)
Flow(
'vprt', None, '''
spike nsp=1 mag=1
n1=%(nz)d o1=%(oz)g d1=%(dz)g k1=%(kz)d l1=%(mz)d
n2=%(nx)d o2=%(ox)g d2=%(dx)g k2=%(kx)d l2=%(mx)d |
smooth rect1=25 rect2=25 repeat=3 |
scale axis=123 |
scale rscale=%(vprt)g
''' % par)
Flow('velo', 'vbck vprt', 'add ${SOURCES[1]}')
Plot('velo', fdmod.cgrey('allpos=y bias=1200 pclip=100 color=F', par))
Result('velo', ['velo', 'ss', 'sx'], 'Overlay')
# density
Flow('dens', 'velo', 'math output=1')
# reflector
Plot(
'refl', 'refl velo', '''
depth2time velocity=${SOURCES[1]} dt=%(dt)g nt=%(nt)d |
scale rscale=-1 |
ricker1 frequency=%(frq)g |
time2depth velocity=${SOURCES[1]} dz=%(dz)g nz=%(nz)d |
''' % par + fdmod.cgrey('pclip=100', par))
Result('refl', ['refl', 'ss', 'sx'], 'Overlay')
# mask
Flow(
'mask', None, '''
spike nsp=1 mag=1
n1=%(nx)d o1=%(ox)g d1=%(dx)g k1=101 l1=900
n2=%(nt)d o2=%(ot)g d2=%(dt)g |
smooth rect1=100 |
scale axis=123
''' % par)
Result('mask', 'transp |' + fdmod.dgrey('allpos=y pclip=100', par))
# F-D modeling (born)
fdmod.lwefd1('do', 'wo', 'dd', 'wd', 'wav', 'velo', 'dens', 'refl', 'ss',
'rr', 'jsnap=100', par)
Result(
'do', 'transp | window min1=0.25 |' +
fdmod.dgrey('min1=0.25 pclip=100', par))
Result(
'dd', 'transp | window min1=0.25 |' +
fdmod.dgrey('min1=0.25 pclip=100', par))
Result('wo', fdmod.wgrey('', par))
Result('wd', fdmod.wgrey('', par))
# source data and wavefield
fdmod.awefd1('ds', 'ws', 'wav', 'velo', 'dens', 'sx', 'rr', '', par)
Result('ws', 'window j3=20 |' + fdmod.wgrey('', par))
# receiver wavefield
Flow('du', 'dd mask',
'add mode=p ${SOURCES[1]} | reverse which=2 opt=i verb=y')
fdmod.awefd('dx', 'wx', 'du', 'velo', 'dens', 'rr', 'rr', '', par)
Flow('dr', 'dx', 'reverse which=2 opt=i verb=y')
Flow('wr', 'wx', 'reverse which=4 opt=i verb=y')
Result('wr', 'window j3=20 |' + fdmod.wgrey('', par))
for i in range(0, par['nt'] / 100, 1):
fdmod.wframe('wo' + '-' + str(i), 'wo', i, 'pclip=99.9', par)
fdmod.wframe('wd' + '-' + str(i), 'wd', i, 'pclip=100', par)
for i in range(0, par['nt'] / 100, 1):
fdmod.wframe('wx' + '-' + str(i), 'wx', i * 25, 'pclip=99.9', par)
# ------------------------------------------------------------
minx = 500
maxx = 1500
minz = par['oz'] + par['nz'] * par['dz'] / 2
numz = par['nz'] / 2
mint = 0.1
numt = 150
maxt = mint + numt * par['dt'] * par['jsnap']
# wavefield
for i in ('s', 'r'):
Flow(
'u' + i, 'w' + i, '''
window min1=%(zmin)g max1=%(zmax)g min2=%(xmin)g max2=%(xmax)g |
scale axis=123
''' % par)
Result('u' + i, 'window j3=10 |' + fdmod.wgrey('pclip=99', par))
for k in range(0, par['nt'] / par['jsnap'], 25):
fdmod.wframe('u' + i + '-' + str(k / 25), 'u' + i, k, 'pclip=99',
par)
# windowed wavefields
Flow(
'p' + i, 'u' + i, '''
window min1=%g n1=%g min2=%g max2=%g min3=%g n3=%g
''' % (minz, numz, minx, maxx, mint, numt))
Flow('q' + i, 'p' + i, 'transp plane=13 memsize=500')
Flow('o' + i, 'q' + i, 'transp plane=23 memsize=500')
Flow('qi', 'qs qr', 'add mode=p ${SOURCES[1]}')
Flow('oi', 'os or', 'add mode=p ${SOURCES[1]}')
for i in ('s', 'r', 'i'):
Result('q' + i,
'window j3=10 |' + fdmod.dgrey('gainpanel=a pclip=100', par))
Result(
'o' + i, 'window j3=10 | transp |' +
fdmod.egrey('gainpanel=a pclip=100', par))
Flow(['q' + i + 'plt', 'q' + i + 'bar'], 'q' + i,
'byte bar=${TARGETS[1]} gainpanel=a pclip=100')
for k in range(10):
Result(
'q' + i + 'plt' + str(k), ['q' + i + 'plt', 'q' + i + 'bar'],
'window n3=1 f3=%d |' % (10 * k) + fdmod.dgrey(
'''
bar=${SOURCES[1]} min1=%g max1=%g min2=%g max2=%g
labelsz=8 labelfat=3 screenratio=1.5
''' % (mint, maxt, minx, maxx), par))
# cut along the reflectors
Flow(
'cut', 'refl', '''
window min1=%g n1=%g min2=%g max2=%g |
spray axis=3 n=%d o=%g d=%g |
transp plane=13 memsize=500
''' % (minz, numz, minx, maxx, numt, 0.1, par['dt'] * 4))
for i in ('s', 'r'):
Flow(
'c' + i, ['q' + i, 'cut'], '''
add mode=p ${SOURCES[1]} |
transp plane=23 |
stack
''')
Result('c' + i, fdmod.dgrey('', par))
Flow('f' + i, 'q' + i, 'window n3=1 min3=300')
Result('f' + i, fdmod.dgrey('', par))
# ------------------------------------------------------------
# conventional IC
Flow('ii', ['ps', 'pr'], 'ic ur=${SOURCES[1]} version=0 nbuf=500 verb=y')
Plot('ii', fdmod.cgrey('pclip=99.9', par))
Result('ii', ['ii', 'ss', 'sx'], 'Overlay')
def test(vp, vs, ro, epsilon, delta, ss, rr, par):
# ------------------------------------------------------------
# source/receiver coordinates
Plot(
rr, 'window n1=2 | dd type=complex | window j2=10 | ' +
fdmod.cgraph('wantscalebar=y symbol=o plotcol=1', par))
Plot(
ss, 'window n1=2 | dd type=complex | window | ' +
fdmod.cgraph('wantscalebar=y symbol=x plotcol=2', par))
# ------------------------------------------------------------
# acoustic source
fdmod.wavelet('wava0', par['frq'], par)
Flow('wava', 'wava0', 'transp')
Result(
'wava', 'transp | window n1=500 |' +
fdmod.waveplot('title="Acoustic source"', par))
# ------------------------------------------------------------
# elastic source
fdmod.wavelet('hor0', par['frq'], par)
fdmod.wavelet('ver0', par['frq'], par)
Flow('hor', 'hor0', 'math output=input*%(hscale)g' % par)
Flow('ver', 'ver0', 'math output=input*%(vscale)g' % par)
Flow('wave0', 'ver hor', 'cat axis=2 space=n ${SOURCES[1:2]}')
Flow(
'wave', 'wave0', '''
transp plane=12 |
transp plane=23 |
transp plane=12
''')
Plot(
'ver', 'wave', 'window n2=1 f2=0 | window n1=500 |' +
fdmod.waveplot('title="Elastic vertical source"', par))
Plot(
'hor', 'wave', 'window n2=1 f2=1 | window n1=500 |' +
fdmod.waveplot('title="Elastic horizontal source"', par))
Result('wave', 'hor ver', 'Movie')
# ------------------------------------------------------------
Plot(vp, fdmod.cgrey('wantscalebar=y allpos=y bias=1.0 pclip=100', par))
Plot(vs, fdmod.cgrey('wantscalebar=y allpos=y bias=1.0 pclip=100', par))
Plot(ro, fdmod.cgrey('wantscalebar=y allpos=y bias=100000 pclip=100', par))
Plot(epsilon,
fdmod.cgrey('wantscalebar=y allpos=y pclip=100', par))
Plot(delta,
fdmod.cgrey('wantscalebar=y allpos=y pclip=100', par))
Result(vp, [vp, ss, rr], 'Overlay')
Result(vs, [vs, ss, rr], 'Overlay')
Result(ro, [ro, ss, rr], 'Overlay')
Result(epsilon, [epsilon, ss, rr], 'Overlay')
Result(delta, [delta, ss, rr], 'Overlay')
fdmod.anisotropic('cc', 'vp', 'vs', 'ro', 'epsilon', 'delta', par)
# ------------------------------------------------------------
# acoustic modeling
amodel('da', 'wa', 'wava', vp, ro, ss, rr, '', par)
Flow(
'waw', 'wa', '''
window min1=%g max1=%g min2=%g max2=%g |
scale axis=123
''' % (par['zmin'], par['zmax'], par['xmin'], par['xmax']))
Result('wa', fdmod.wgrey('pclip=99 title="Acoustic wavefield"', par))
Result(
'da',
'transp | window f1=%(kt)d | put o1=%(ot)g | pad end1=%(kt)d |' % par +
fdmod.dgrey('pclip=99 title="Acoustic data" grid=y', par))
# elastic modeling
emodel('de', 'we', 'wave', 'cc', ro, ss, rr, 'ssou=%(ssou)s opot=n' % par,
par)
for i in range(2):
Flow(
'we' + str(i + 1), 'we', '''
window n3=1 f3=%d |
window min1=%g max1=%g min2=%g max2=%g |
scale axis=123
''' % (i, par['zmin'], par['zmax'], par['xmin'], par['xmax']))
Result('we' + str(i + 1),
fdmod.wgrey('title=u%s pclip=99' % str(i + 1), par))
Result(
'de' + str(i + 1), 'de', '''
window n2=1 f2=%d |
transp |
window f1=%d | put o1=%g | pad end1=%d |
''' % (i, par['kt'], par['ot'], par['kt']) +
fdmod.dgrey('title=u%s pclip=99 grid=y' % str(i + 1), par))
Flow('weall', 'we1 we2', 'cat axis=1 space=n ${SOURCES[1]}')
Result(
'weall', '''
grey title="Elastic Wavefields" wantaxis=y screenratio=%f screenht=8
gainpanel=a pclip=99
grid1=y grid2=y g1num=0.25 g2num=0.25
''' % (2 * par['ratio']))
Flow('wall', 'waw we1 we2', 'cat axis=1 space=n ${SOURCES[1:3]}')
Result(
'wall', '''
grey title="" wantaxis=y screenratio=%f screenht=10
gainpanel=a pclip=99
grid1=y grid2=y g1num=0.1 g2num=0.1
''' % (3 * par['ratio']))
# wavefield movie frames
for j in range(0, par['nt'] / par['jsnap'], 1):
fdmod.wframe('wa-' + str(j), 'wa', j, 'pclip=99.9', par)
fdmod.wframe('we1-' + str(j), 'we1', j, 'pclip=99.9', par)
fdmod.wframe('we2-' + str(j), 'we2', j, 'pclip=99.9', par)
