def points(sp,sa,wp,xm,qq,par):
# ------------------------------------------------------------
ls = 2*par['lo']
par['zs' ] = par['oz'] + par['lo']
par['zs0'] = par['zt']+ls
par['zs1'] = par['zt']
par['zs2'] = par['zt']-ls
par['xs' ] = par['xt']
par['xs0'] = par['xt']
par['xs1'] = par['xt']-ls
par['xs2'] = par['xt']+ls
par['jzs0'] = par['zs0'] / par['dz']
par['jzs1'] = par['zs1'] / par['dz']
par['jzs2'] = par['zs2'] / par['dz']
par['jxs0'] = par['xs0'] / par['dx']
par['jxs1'] = par['xs1'] / par['dx']
par['jxs2'] = par['xs2'] / par['dx']
# ------------------------------------------------------------
# source positions
fdmod.point('ss0',par['xs0'],par['zs0'],par)
fdmod.point('ss1',par['xs1'],par['zs1'],par)
fdmod.point('ss2',par['xs2'],par['zs2'],par)
Flow(sp,'ss0 ss1 ss2','cat axis=2 space=n ${SOURCES[1:3]} | window n1=2')
Plot(sp,'window |' + fdmod.ssplot('plotcol=2',par))
fdmod.point(sa,par['xs'],par['zs'],par)
Plot(sa,'window |' + fdmod.ssplot('',par))
# ------------------------------------------------------------
# source wavelet
fdmod.wavelet(wp+'_',par['fo'],par)
Flow(wp,wp+'_','spray axis=2 n=3 o=0 d=1 | transp')
# ------------------------------------------------------------
# scatterrers
Flow(xm,None,
'''
spike nsp=3 mag=1,1,1
n1=%(nz)d o1=%(oz)g d1=%(dz)g k1=%(jzs0)d,%(jzs1)d,%(jzs2)d l1=%(jzs0)d,%(jzs1)d,%(jzs2)d
n2=%(nx)d o2=%(ox)g d2=%(dx)g k2=%(jxs0)d,%(jxs1)d,%(jxs2)d l2=%(jxs0)d,%(jxs1)d,%(jxs2)d
''' % par)
Plot( xm,'smooth rect1=3 rect2=3 repeat=3 |'
+fdmod.cgrey('pclip=100',par))
Result(xm,[xm,qq],'Overlay')
def points(sp, sa, wp, xm, qq, par):
# ------------------------------------------------------------
ls = 2 * par['lo']
par['zs'] = par['oz'] + par['lo']
par['zs0'] = par['zt'] + ls
par['zs1'] = par['zt']
par['zs2'] = par['zt'] - ls
par['xs'] = par['xt']
par['xs0'] = par['xt']
par['xs1'] = par['xt'] - ls
par['xs2'] = par['xt'] + ls
par['jzs0'] = par['zs0'] / par['dz']
par['jzs1'] = par['zs1'] / par['dz']
par['jzs2'] = par['zs2'] / par['dz']
par['jxs0'] = par['xs0'] / par['dx']
par['jxs1'] = par['xs1'] / par['dx']
par['jxs2'] = par['xs2'] / par['dx']
# ------------------------------------------------------------
# source positions
fdmod.point('ss0', par['xs0'], par['zs0'], par)
fdmod.point('ss1', par['xs1'], par['zs1'], par)
fdmod.point('ss2', par['xs2'], par['zs2'], par)
Flow(sp, 'ss0 ss1 ss2', 'cat axis=2 space=n ${SOURCES[1:3]} | window n1=2')
Plot(sp, 'window |' + fdmod.ssplot('plotcol=2', par))
fdmod.point(sa, par['xs'], par['zs'], par)
Plot(sa, 'window |' + fdmod.ssplot('', par))
# ------------------------------------------------------------
# source wavelet
fdmod.wavelet(wp + '_', par['fo'], par)
Flow(wp, wp + '_', 'spray axis=2 n=3 o=0 d=1 | transp')
# ------------------------------------------------------------
# scatterrers
Flow(
xm, None, '''
spike nsp=3 mag=1,1,1
n1=%(nz)d o1=%(oz)g d1=%(dz)g k1=%(jzs0)d,%(jzs1)d,%(jzs2)d l1=%(jzs0)d,%(jzs1)d,%(jzs2)d
n2=%(nx)d o2=%(ox)g d2=%(dx)g k2=%(jxs0)d,%(jxs1)d,%(jxs2)d l2=%(jxs0)d,%(jxs1)d,%(jxs2)d
''' % par)
Plot(xm,
'smooth rect1=3 rect2=3 repeat=3 |' + fdmod.cgrey('pclip=100', par))
Result(xm, [xm, qq], 'Overlay')
def thrsou(ss,wav,par):
par['xsou0']=par['xsou']-0.25
par['zsou0']=par['zsou']-0.25
par['xsou1']=par['xsou']
par['zsou1']=par['zsou']
par['xsou2']=par['xsou']+0.25
par['zsou2']=par['zsou']+0.25
for j in range(3):
fdmod.point(ss+str(j),par['xsou'+str(j)],par['zsou'+str(j)],par)
Flow(ss,[ss+'0',ss+'1',ss+'2'],'cat axis=2 space=n ${SOURCES[1:3]}')
fdmod.wavelet(wav+'0_',50,par)
Flow(wav+'0',
wav+'0_',
'scale rscale=1.0 | window f1=100 | pad end1=100 | put o1=%(ot)g n2=1' %par)
fdmod.wavelet(wav+'1_',50,par)
Flow(wav+'1',
wav+'1_',
'window' %par)
fdmod.wavelet(wav+'2_',50,par)
Flow(wav+'2',
wav+'2_',
'scale rscale=1.0 | pad beg1=100 | window n1=%(nt)d | put o1=%(ot)g n2=1' %par)
Flow(wav,[wav+'0',wav+'1',wav+'2'],
'''
cat axis=2 space=n ${SOURCES[1:3]} |
transp
''')
def point(sp,sa,wp,xm,qq,par):
# ------------------------------------------------------------
ls = 2*par['lo']
par['zs' ] = par['oz'] + par['lo']
par['zs0'] = par['zt']
par['xs' ] = par['xt']
par['xs0'] = par['xt']
par['jzs0'] = (par['zs0']-par['oz']) / par['dz']
par['jxs0'] = (par['xs0']-par['ox']) / par['dx']
# ------------------------------------------------------------
# source positions
fdmod.point(sp,par['xs0'],par['zs0'],par)
Plot(sp,'window |' + fdmod.ssplot('',par))
fdmod.point(sa,par['xs'],par['zs'],par)
Plot(sa,'window |' + fdmod.ssplot('',par))
# ------------------------------------------------------------
# source wavelet
fdmod.wavelet(wp+'_',par['fo'],par)
Flow(wp,wp+'_','transp')
# ------------------------------------------------------------
# scatterrers
Flow(xm,None,
'''
spike nsp=3 mag=1,1,1
n1=%(nz)d o1=%(oz)g d1=%(dz)g k1=%(jzs0)d l1=%(jzs0)d
n2=%(nx)d o2=%(ox)g d2=%(dx)g k2=%(jxs0)d l2=%(jxs0)d
''' % par)
Plot( xm,'smooth rect1=3 rect2=3 repeat=3 |'
+fdmod.cgrey('pclip=100',par))
Result(xm,[xm,qq],'Overlay')
def point(sp, sa, wp, xm, qq, par):
# ------------------------------------------------------------
ls = 2 * par['lo']
par['zs'] = par['oz'] + par['lo']
par['zs0'] = par['zt']
par['xs'] = par['xt']
par['xs0'] = par['xt']
par['jzs0'] = (par['zs0'] - par['oz']) / par['dz']
par['jxs0'] = (par['xs0'] - par['ox']) / par['dx']
# ------------------------------------------------------------
# source positions
fdmod.point(sp, par['xs0'], par['zs0'], par)
Plot(sp, 'window |' + fdmod.ssplot('', par))
fdmod.point(sa, par['xs'], par['zs'], par)
Plot(sa, 'window |' + fdmod.ssplot('', par))
# ------------------------------------------------------------
# source wavelet
fdmod.wavelet(wp + '_', par['fo'], par)
Flow(wp, wp + '_', 'transp')
# ------------------------------------------------------------
# scatterrers
Flow(
xm, None, '''
spike nsp=3 mag=1,1,1
n1=%(nz)d o1=%(oz)g d1=%(dz)g k1=%(jzs0)d l1=%(jzs0)d
n2=%(nx)d o2=%(ox)g d2=%(dx)g k2=%(jxs0)d l2=%(jxs0)d
''' % par)
Plot(xm,
'smooth rect1=3 rect2=3 repeat=3 |' + fdmod.cgrey('pclip=100', par))
Result(xm, [xm, qq], 'Overlay')
def thrsou(ss, wav, par):
par['xsou0'] = par['xsou'] - 0.25
par['zsou0'] = par['zsou'] - 0.25
par['xsou1'] = par['xsou']
par['zsou1'] = par['zsou']
par['xsou2'] = par['xsou'] + 0.25
par['zsou2'] = par['zsou'] + 0.25
for j in range(3):
fdmod.point(ss + str(j), par['xsou' + str(j)], par['zsou' + str(j)],
par)
Flow(ss, [ss + '0', ss + '1', ss + '2'],
'cat axis=2 space=n ${SOURCES[1:3]}')
fdmod.wavelet(wav + '0_', 50, par)
Flow(
wav + '0', wav + '0_',
'scale rscale=1.0 | window f1=100 | pad end1=100 | put o1=%(ot)g n2=1'
% par)
fdmod.wavelet(wav + '1_', 50, par)
Flow(wav + '1', wav + '1_', 'window' % par)
fdmod.wavelet(wav + '2_', 50, par)
Flow(
wav + '2', wav + '2_',
'scale rscale=1.0 | pad beg1=100 | window n1=%(nt)d | put o1=%(ot)g n2=1'
% par)
Flow(
wav, [wav + '0', wav + '1', wav + '2'], '''
cat axis=2 space=n ${SOURCES[1:3]} |
transp
''')
def reflectivity(refl, sp, sa, wp, xm, qq, par):
par['zs'] = par['lo']
par['zs0'] = par['zt']
par['xs'] = par['ox'] + 0.50 * par['nx'] * par['dx']
par['xs0'] = par['xt']
par['jzs0'] = 0.85 * par['nz']
par['jxs0'] = 0.50 * par['nx']
Flow(
refl + '_mask', refl, '''
window f1=300 |
mask min=-0.01 max=0.01 |
dd type=float |
scale rscale=-1 |
add add=1 |
dd type=integer |
pad beg1=300
''')
Result(refl + '_mask', 'sfdd type=float |' + fdmod.cgrey('pclip=100', par))
par['nspk'] = par['nz'] * par['nx']
for x in ('x1', 'x2'):
if (x == 'x1'): o = sp + '_z'
if (x == 'x2'): o = sp + '_x'
Flow(
o, [refl, refl + '_mask'], '''
math output=%s |
put n1=1 n2=%d |
headerwindow mask=${SOURCES[1]} |
window
''' % (x, par['nspk']))
Flow(sp, [sp + '_x', sp + '_z'],
'''
cat axis=2 space=n
${SOURCES[0:2]} | transp
''',
stdin=0)
Plot(sp, 'window |' + fdmod.ssplot('', par))
fdmod.point(sa, par['xs'], par['zs'], par)
Plot(sa, 'window |' + fdmod.ssplot('', par))
# ------------------------------------------------------------
# source wavelet
Flow(
sp + '_r', [refl, refl + '_mask'], '''
put n1=1 n2=%d |
headerwindow mask=${SOURCES[1]} |
window
''' % (par['nspk']))
Flow(
sp + '_mask', sp + '_r', '''
spray axis=2 n=%(nt)d o=%(ot)g d=%(dt)g
''' % par)
fdmod.wavelet(wp + '_', par['fo'], par)
Flow(
wp, [wp + '_', sp + '_mask'], '''
spray axis=2 n=32914 o=0 d=1 |
transp |
math m=${SOURCES[1]} output="input*m"
''' % par)
# ------------------------------------------------------------
# scatterers
Flow(xm, refl, 'window')
Plot(xm,
'smooth rect1=3 rect2=3 repeat=3 |' + fdmod.cgrey('pclip=98', par))
Result(xm, [xm, qq], 'Overlay')
def segments(sp, sa, wp, xm, qq, par):
par['ff'] = 30
par['zs'] = par['oz'] + par['lo']
par['zs0'] = par['zt']
par['xs'] = par['ox'] + 0.35 * par['nx'] * par['dx']
par['xs0'] = par['xt']
par['jzs0'] = 0.85 * par['nz']
par['jxs0'] = 0.50 * par['nx']
n = 200
m = 26
par['nspk'] = n * m
par['allx'] = ''
par['allz'] = ''
par['allr'] = ''
par['alli'] = ''
ind = 0
for j in range(m):
for i in range(n):
xc = par['jxs0'] + (
i - n / 2) - (j - m / 2 + 0.5) * 7 * par['lo'] / par['dx']
zc = par['jzs0'] + (i - n / 2) * math.tan(math.radians(par['ff']))
par['allx'] += '%d,' % xc
par['allz'] += '%d,' % zc
par['allr'] += '1,'
ind += 1
par['alli'] += '%d,' % ind
# ------------------------------------------------------------
# source positions
Flow(sp + '_x', None,
'spike nsp=%(nspk)s mag=%(allx)s n1=%(nspk)s k1=%(alli)s' % par)
Flow(sp + '_z', None,
'spike nsp=%(nspk)s mag=%(allz)s n1=%(nspk)s k1=%(alli)s' % par)
Flow(sp + '_r', None,
'spike nsp=%(nspk)s mag=%(allr)s n1=%(nspk)s k1=%(alli)s' % par)
Flow(sp, [sp + '_x', sp + '_z'],
'''
cat axis=2 space=n
${SOURCES[0:2]} | transp
''',
stdin=0)
Plot(sp, 'window |' + fdmod.ssplot('', par))
fdmod.point(sa, par['xs'], par['zs'], par)
Plot(sa, 'window |' + fdmod.ssplot('', par))
# ------------------------------------------------------------
# source wavelet
fdmod.wavelet(wp + '_', par['fo'], par)
Flow(wp, wp + '_', 'transp' % par)
# Flow(wp,wp+'_','spray axis=2 n=%(nspk)d o=0 d=1 | transp' % par)
# ------------------------------------------------------------
# scatterers
Flow(
xm, None, '''
spike nsp=%(nspk)s mag=%(allr)s
n1=%(nz)d o1=%(oz)g d1=%(dz)g k1=%(allz)s
n2=%(nx)d o2=%(ox)g d2=%(dx)g k2=%(allx)s
''' % par)
Plot(xm,
'smooth rect1=5 rect2=5 repeat=1 |' + fdmod.cgrey('pclip=98', par))
Result(xm, [xm, qq], 'Overlay')
def run(par):
# ------------------------------------------------------------
fdmod.point('ss1',par['xsou1'],par['oz'],par)
fdmod.point('ss2',par['xsou2'],par['oz'],par)
Plot('ss1','window |' + fdmod.ssplot('plotcol=5',par))
Plot('ss2','window |' + fdmod.ssplot('plotcol=5',par))
# ------------------------------------------------------------
# velocity
Plot( 'vel',fdmod.cgrey('bias=4.8 allpos=y pclip=99 color=j',par))
Result('vel','vel ss1 ss2','Overlay')
# slowness
Flow('slo','vel',
'''
math output=1/input |
transp |
transp plane=23 |
put o2=0 d2=1 label2=y
''')
Result('slo','window | transp |'
+ fdmod.cgrey('allpos=y pclip=95 bias=0.125',par))
# reflectivity
Flow('ref','del',
'''
transp |
transp plane=23 |
put o2=0 d2=1 label2=y
''')
Result('ref','window | transp |' + fdmod.cgrey('pclip=99',par))
# ------------------------------------------------------------
# wavelet
fdmod.wavelet('wav',8,par)
Result('wav','window n1=500 |' + fdmod.waveplot('',par))
# ------------------------------------------------------------
for i in ('1','2'):
Flow('spk'+i,'wav',
'''
pad beg2=%d n2out=%d |
put o2=%g d2=%g
''' % ( (par['xsou'+i]-par['ox'])/par['dx'],par['nx'],par['ox'],par['dx']) )
Result('spk'+i,fdmod.dgrey('pclip=100',par))
# ------------------------------------------------------------
for i in ('1','2'):
# source wavefield (from time to frequency)
zomig.wflds('dds'+i,'spk'+i,par)
# wavefield extrapolation MODELING
spmig.modelPW3('ddr'+i,'slo','dds'+i,'ref',par)
# ------------------------------------------------------------
# shots 1 and 2
Flow('dds0','dds1 dds2','add ${SOURCES[1]}')
Flow('ddr0','ddr1 ddr2','add ${SOURCES[1]}') # both shots
for i in ('1','2','0'):
Result('dds'+i,'window | real | smooth rect1=5 | sfgrey pclip=100')
Result('ddr'+i,'window | real | smooth rect1=5 | sfgrey pclip=100')
# ------------------------------------------------------------
# available dds[0,1,2], ddr[0,1,2], slo
# ------------------------------------------------------------
for i in ('1','2','0'):
# recorded data (from frequency to time)
Flow('ttr'+i,'ddr'+i,
'''
window |
transp |
pad beg1=%(fw)d n1out=%(kw)d |
fft1 inv=y opt=n
''' % par)
Result('ttr'+i,fdmod.dgrey('pclip=100 min1=1 max1=6 screenratio=0.5 screenht=7',par))
# wavefield extrapolation MIGRATION
spmig.imagePW3('ii'+i,'cc'+i,'slo','dds'+i,'ddr'+i,par)
Plot( 'ii'+i,'window min1=%(lox)d max1=%(hix)d | transp |' %par
+ fdmod.cgrey('pclip=99',par))
Result('ii'+i,['ii'+i,'ss1','ss2'],'Overlay')
Flow('ii','ii1 ii2','add ${SOURCES[1]}')
Plot('ii','window min1=%(lox)d max1=%(hix)d | transp |' %par
+ fdmod.cgrey('pclip=100',par))
Result('ii',['ii','ss1','ss2'],'Overlay')
# ------------------------------------------------------------
# datuming
zomig.Cwfone3('wfs','dds0','slo',par) # source wavefield for one shot
zomig.Awfone3('wfr','ddr0','slo',par) # receiver wavefield for two shots
# ------------------------------------------------------------
# data in the time domain
for k in ('s','r'):
Flow('q'+k,'wf'+k,
'''
window min1=%(lox)d max1=%(hix)d j1=2 j3=2 |
transp plane=23 memsize=500 |
transp plane=12 memsize=500 |
pad beg1=%(fw)d n1out=%(kw)d |
fft1 inv=y opt=n |
window max1=%(tcut)g
''' % par)
Plot('q'+k,'window j3=10 |' +
fdmod.dgrey('gainpanel=a pclip=99',par),view=1)
Flow( 'qi','qs qr','add mode=p ${SOURCES[1]}')
Plot('qi','window j3=10 |'
+ fdmod.dgrey('gainpanel=a pclip=100',par),view=1)
# SIC
Flow('kk',['qs','qr'],
'''
sic3d ur=${SOURCES[1]} nbuf=500 verb=y stack=n
oanga=%(oanga)g nanga=%(nanga)d danga=%(danga)g
oangb=%(oangb)g nangb=%(nangb)d dangb=%(dangb)g
nl=%(nl)d dl=%(dl)g
sig=%(sig)g
''' % par)
Plot('kk','transp plane=23 | stack | transp |'
+ fdmod.cgrey('pclip=100',par))
Result('kk',['kk','ss1','ss2'],'Overlay')
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=g', 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)
Plot('u' + i, 'window j3=10 |' + fdmod.wgrey('pclip=99', par), view=1)
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'):
Plot('q' + i,
'window j3=10 |' + fdmod.dgrey('gainpanel=a pclip=100', par),
view=1)
Plot('o' + i,
'window j3=10 | transp |' +
fdmod.egrey('gainpanel=a pclip=100', par),
view=1)
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 run(par):
# ------------------------------------------------------------
fdmod.point('ss1', par['xsou1'], par['oz'], par)
fdmod.point('ss2', par['xsou2'], par['oz'], par)
Plot('ss1', 'window |' + fdmod.ssplot('plotcol=5', par))
Plot('ss2', 'window |' + fdmod.ssplot('plotcol=5', par))
# ------------------------------------------------------------
# velocity
Plot('vel', fdmod.cgrey('bias=4.8 allpos=y pclip=99 color=j', par))
Result('vel', 'vel ss1 ss2', 'Overlay')
# slowness
Flow(
'slo', 'vel', '''
math output=1/input |
transp |
transp plane=23 |
put o2=0 d2=1 label2=y
''')
Result(
'slo',
'window | transp |' + fdmod.cgrey('allpos=y pclip=95 bias=0.125', par))
# reflectivity
Flow(
'ref', 'del', '''
transp |
transp plane=23 |
put o2=0 d2=1 label2=y
''')
Result('ref', 'window | transp |' + fdmod.cgrey('pclip=99', par))
# ------------------------------------------------------------
# wavelet
fdmod.wavelet('wav', 8, par)
Result('wav', 'window n1=500 |' + fdmod.waveplot('', par))
# ------------------------------------------------------------
for i in ('1', '2'):
Flow(
'spk' + i, 'wav', '''
pad beg2=%d n2out=%d |
put o2=%g d2=%g
''' % ((par['xsou' + i] - par['ox']) / par['dx'], par['nx'],
par['ox'], par['dx']))
Result('spk' + i, fdmod.dgrey('pclip=100', par))
# ------------------------------------------------------------
for i in ('1', '2'):
# source wavefield (from time to frequency)
zomig.wflds('dds' + i, 'spk' + i, par)
# wavefield extrapolation MODELING
spmig.modelPW3('ddr' + i, 'slo', 'dds' + i, 'ref', par)
# ------------------------------------------------------------
# shots 1 and 2
Flow('dds0', 'dds1 dds2', 'add ${SOURCES[1]}')
Flow('ddr0', 'ddr1 ddr2', 'add ${SOURCES[1]}') # both shots
for i in ('1', '2', '0'):
Result('dds' + i, 'window | real | smooth rect1=5 | sfgrey pclip=100')
Result('ddr' + i, 'window | real | smooth rect1=5 | sfgrey pclip=100')
# ------------------------------------------------------------
# available dds[0,1,2], ddr[0,1,2], slo
# ------------------------------------------------------------
for i in ('1', '2', '0'):
# recorded data (from frequency to time)
Flow(
'ttr' + i, 'ddr' + i, '''
window |
transp |
pad beg1=%(fw)d n1out=%(kw)d |
fft1 inv=y opt=n
''' % par)
Result(
'ttr' + i,
fdmod.dgrey('pclip=100 min1=1 max1=6 screenratio=0.5 screenht=7',
par))
# wavefield extrapolation MIGRATION
spmig.imagePW3('ii' + i, 'cc' + i, 'slo', 'dds' + i, 'ddr' + i, par)
Plot(
'ii' + i, 'window min1=%(lox)d max1=%(hix)d | transp |' % par +
fdmod.cgrey('pclip=99', par))
Result('ii' + i, ['ii' + i, 'ss1', 'ss2'], 'Overlay')
Flow('ii', 'ii1 ii2', 'add ${SOURCES[1]}')
Plot(
'ii', 'window min1=%(lox)d max1=%(hix)d | transp |' % par +
fdmod.cgrey('pclip=100', par))
Result('ii', ['ii', 'ss1', 'ss2'], 'Overlay')
# ------------------------------------------------------------
# datuming
zomig.Cwfone3('wfs', 'dds0', 'slo', par) # source wavefield for one shot
zomig.Awfone3('wfr', 'ddr0', 'slo',
par) # receiver wavefield for two shots
# ------------------------------------------------------------
# data in the time domain
for k in ('s', 'r'):
Flow(
'q' + k, 'wf' + k, '''
window min1=%(lox)d max1=%(hix)d j1=2 j3=2 |
transp plane=23 memsize=500 |
transp plane=12 memsize=500 |
pad beg1=%(fw)d n1out=%(kw)d |
fft1 inv=y opt=n |
window max1=%(tcut)g
''' % par)
Plot('q' + k,
'window j3=10 |' + fdmod.dgrey('gainpanel=a pclip=99', par),
view=1)
Flow('qi', 'qs qr', 'add mode=p ${SOURCES[1]}')
Plot('qi',
'window j3=10 |' + fdmod.dgrey('gainpanel=a pclip=100', par),
view=1)
# SIC
Flow(
'kk', ['qs', 'qr'], '''
sic3d ur=${SOURCES[1]} nbuf=500 verb=y stack=n
oanga=%(oanga)g nanga=%(nanga)d danga=%(danga)g
oangb=%(oangb)g nangb=%(nangb)d dangb=%(dangb)g
nl=%(nl)d dl=%(dl)g
sig=%(sig)g
''' % par)
Plot('kk',
'transp plane=23 | stack | transp |' + fdmod.cgrey('pclip=100', par))
Result('kk', ['kk', 'ss1', 'ss2'], 'Overlay')
def reflectivity(refl,sp,sa,wp,xm,qq,par):
par['zs' ] = par['lo']
par['zs0'] = par['zt']
par['xs' ] = par['ox']+0.50*par['nx']*par['dx']
par['xs0'] = par['xt']
par['jzs0'] = 0.85*par['nz']
par['jxs0'] = 0.50*par['nx']
Flow(refl+'_mask',refl,
'''
window f1=300 |
mask min=-0.01 max=0.01 |
dd type=float |
scale rscale=-1 |
add add=1 |
dd type=integer |
pad beg1=300
''')
Result(refl+'_mask',
'sfdd type=float |' + fdmod.cgrey('pclip=100',par))
par['nspk']=par['nz']*par['nx']
for x in ('x1','x2'):
if(x=='x1'): o=sp+'_z'
if(x=='x2'): o=sp+'_x'
Flow(o,[refl,refl+'_mask'],
'''
math output=%s |
put n1=1 n2=%d |
headerwindow mask=${SOURCES[1]} |
window
''' % (x,par['nspk']))
Flow(sp,[sp+'_x',sp+'_z'],
'''
cat axis=2 space=n
${SOURCES[0:2]} | transp
''', stdin=0)
Plot(sp,'window |' + fdmod.ssplot('',par))
fdmod.point(sa,par['xs'],par['zs'],par)
Plot(sa,'window |' + fdmod.ssplot('',par))
# ------------------------------------------------------------
# source wavelet
Flow(sp+'_r',[refl,refl+'_mask'],
'''
put n1=1 n2=%d |
headerwindow mask=${SOURCES[1]} |
window
''' % (par['nspk']) )
Flow(sp+'_mask',sp+'_r',
'''
spray axis=2 n=%(nt)d o=%(ot)g d=%(dt)g
''' % par)
fdmod.wavelet(wp+'_',par['fo'],par)
Flow(wp,[wp+'_',sp+'_mask'],
'''
spray axis=2 n=32914 o=0 d=1 |
transp |
math m=${SOURCES[1]} output="input*m"
''' % par)
# ------------------------------------------------------------
# scatterers
Flow(xm,refl,'window')
Plot( xm,'smooth rect1=3 rect2=3 repeat=3 |'
+fdmod.cgrey('pclip=98',par))
Result(xm,[xm,qq],'Overlay')
def segments(sp,sa,wp,xm,qq,par):
par['ff']=30
par['zs' ] = par['oz'] + par['lo']
par['zs0'] = par['zt']
par['xs' ] = par['ox']+0.35*par['nx']*par['dx']
par['xs0'] = par['xt']
par['jzs0'] = 0.85*par['nz']
par['jxs0'] = 0.50*par['nx']
n = 200
m = 26
par['nspk'] = n*m
par['allx'] = ''
par['allz'] = ''
par['allr'] = ''
par['alli'] = ''
ind = 0
for j in range(m):
for i in range(n):
xc = par['jxs0'] + (i-n/2) - (j-m/2+0.5) * 7*par['lo']/ par['dx']
zc = par['jzs0'] + (i-n/2) * math.tan(math.radians(par['ff']))
par['allx'] += '%d,'%xc
par['allz'] += '%d,'%zc
par['allr'] += '1,'
ind +=1
par['alli'] += '%d,'%ind
# ------------------------------------------------------------
# source positions
Flow(sp+'_x',None,'spike nsp=%(nspk)s mag=%(allx)s n1=%(nspk)s k1=%(alli)s' %par)
Flow(sp+'_z',None,'spike nsp=%(nspk)s mag=%(allz)s n1=%(nspk)s k1=%(alli)s' %par)
Flow(sp+'_r',None,'spike nsp=%(nspk)s mag=%(allr)s n1=%(nspk)s k1=%(alli)s' %par)
Flow(sp,[sp+'_x',sp+'_z'],
'''
cat axis=2 space=n
${SOURCES[0:2]} | transp
''', stdin=0)
Plot(sp,'window |' + fdmod.ssplot('',par))
fdmod.point(sa,par['xs'],par['zs'],par)
Plot(sa,'window |' + fdmod.ssplot('',par))
# ------------------------------------------------------------
# source wavelet
fdmod.wavelet(wp+'_',par['fo'],par)
Flow(wp,wp+'_','transp' % par)
# Flow(wp,wp+'_','spray axis=2 n=%(nspk)d o=0 d=1 | transp' % par)
# ------------------------------------------------------------
# scatterers
Flow(xm,None,
'''
spike nsp=%(nspk)s mag=%(allr)s
n1=%(nz)d o1=%(oz)g d1=%(dz)g k1=%(allz)s
n2=%(nx)d o2=%(ox)g d2=%(dx)g k2=%(allx)s
''' % par)
Plot( xm,'smooth rect1=5 rect2=5 repeat=1 |'
+fdmod.cgrey('pclip=98',par))
Result(xm,[xm,qq],'Overlay')
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=g',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))
Plot('wo',fdmod.wgrey('',par),view=1)
Plot('wd',fdmod.wgrey('',par),view=1)
# source data and wavefield
fdmod.awefd1(
'ds','ws',
'wav','velo','dens','sx','rr','',par)
Plot('ws','window j3=20 |' + fdmod.wgrey('',par),view=1)
# 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')
Plot('wr','window j3=20 |' + fdmod.wgrey('',par),view=1)
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)
Plot('u'+i,'window j3=10 |' + fdmod.wgrey('pclip=99',par),view=1)
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'):
Plot('q'+i,'window j3=10 |'
+ fdmod.dgrey('gainpanel=a pclip=100',par),view=1)
Plot('o'+i,'window j3=10 | transp |'
+ fdmod.egrey('gainpanel=a pclip=100',par),view=1)
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 onesou(ss, qq, wav, par):
fdmod.point(ss, par['xsou'], par['zsou'], par)
fdmod.wavelet(wav + '_', 50, par)
Flow(wav, wav + '_', 'transp')
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)
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)
def wavelet(wav, par):
fdmod.wavelet(wav + '_', par['fo'], par)
Flow(wav, wav + '_', 'transp')
def onesou(ss,qq,wav,par):
fdmod.point(ss,par['xsou'],par['zsou'],par)
fdmod.wavelet(wav+'_',50,par)
Flow(wav,wav+'_','transp')
def wavelet(wav,par):
fdmod.wavelet(wav+'_',par['fo'],par)
Flow(wav,wav+'_','transp')
