Skip to content
Snippets Groups Projects
Commit d4d85049 authored by JanThorbecke's avatar JanThorbecke
Browse files

demo Marchenko

parent 43e64a59
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
marchenko/demo/README 0 → 100644
+ 9
− 0
View file @ d4d85049
Description of files:
1) shots.scr create the shots
2) model.scr computes the model
3) direct_wave.scr crate the direct wave to be removed from the shots
4) remove_direct.scr remove the direct wave from the shots and scale them
5) first_arrival.scr computes the first arrival
6) marchenko.scr perform the Marchenko scheme
7) referenceShot.scr creates the reference Green's function
marchenko/demo/direct.scr 0 → 100755
+ 37
− 0
View file @ d4d85049
#!/bin/bash
export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
dx=2.5
dt=0.0005
makemod sizex=12000 sizez=4000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
orig=-6000,-1000 file_base=noContrast.su
export OMP_NUM_THREADS=8
makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3 scale=0
$HOME/bin/fdelmodc \
file_cp=noContrast_cp.su ischeme=1 iorder=4 \
file_den=noContrast_ro.su \
file_src=wavefw.su \
file_rcv=direct.su \
src_type=7 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=1 \
rec_type_p=1 \
rec_int_vz=2 \
rec_delay=0.3 \
dtrcv=0.004 \
verbose=2 \
tmod=4.394 \
dxrcv=10.0 \
xrcv1=-6000 xrcv2=6000 \
zrcv1=0 zrcv2=0 \
xsrc=0 zsrc=0 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
marchenko/demo/first_arrival.scr 0 → 100755
+ 92
− 0
View file @ d4d85049
#!/bin/bash
#PBS -N fdelmod
#PBS -q verylong
#PBS -l nodes=1
#PBS -k eo
#PBS -j eo
export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
dx=2.5
dt=0.0005
makemod sizex=6000 sizez=2000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
orig=-3000,0 file_base=synclDown.su verbose=2 \
intt=def x=-3000,500,3000 z=195,195,195 poly=1 cp=1950 ro=3700 \
intt=def x=-3000,3000 z=600,600 poly=0 cp=2050 ro=1750 \
intt=def x=-3000,3000 z=680,680 poly=0 cp=2150 ro=2220 \
intt=def x=-3000,3000 z=780,780 poly=0 cp=2230 ro=1700 \
intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=520,580,680,840,680,600,500 poly=2 cp=2400 ro=2800 \
# intt=def x=-3000,0,3000 z=1110,1110,1110 poly=0 cp=2300 ro=1950 \
# intt=def x=-3000,3000 z=1180,1180 poly=0 cp=2480 ro=1820 \
# intt=def x=-3000,0,3000 z=1290,1290,1370 poly=0 cp=2600 ro=2000 \
# intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2720 ro=2050 \
# intt=def x=-3000,3000 z=1480,1480 poly=0 cp=2800 ro=1850
makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
#smooth file_in=synclDown_cp.su power=-1.0 ntsm=29 nxsm=29 niter=15 file_out=syncls_cp.su
#smooth file_in=synclDown_ro.su power=-1.0 ntsm=29 nxsm=29 niter=15 file_out=syncls_ro.su
dxshot=10
ishot=300
nshots=301
export OMP_NUM_THREADS=1
mkdir -p shots
mkdir -p jobs
while (( ishot < nshots ))
do
(( xsrc = -3000 + ${ishot}*${dxshot} ))
# (( xsrc = -1100 + ${ishot}*${dxshot} ))
echo xsrc=$xsrc
file_rcv=shots/shotsmonPz1100_${xsrc}.su
cat << EOF > jobs/pbs_$ishot.job
#!/bin/bash
#
#PBS -q medium
#PBS -N mod_${xsrc}
#PBS -j eo
#PBS -m n
#PBS -l nodes=1
#PBS -V
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
export OMP_NUM_THREADS=1
$HOME/bin64/fdelmodc \
file_cp=synclDown_cp.su ischeme=1 iorder=4 \
file_den=synclDown_ro.su \
file_src=wave.su \
file_rcv=$file_rcv \
src_type=1 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_vz=2 \
rec_delay=0.1 \
dtrcv=0.004 \
verbose=2 \
tmod=2.100 \
dxrcv=10.0 \
xrcv1=-3000 xrcv2=3000 \
zrcv1=0 zrcv2=0 \
xsrc=$xsrc zsrc=1100 \
ntaper=300 \
left=4 right=4 top=4 bottom=4
EOF
qsub jobs/pbs_$ishot.job
(( ishot = $ishot + 1))
done
marchenko/demo/marchenko.scr 0 → 100755
+ 39
− 0
View file @ d4d85049
#!/bin/bash -x
#PBS -N fdelmod
#PBS -q verylong
#PBS -l nodes=1
#PBS -k eo
#PBS -j eo
export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
which makewave
which makemod
which fmute
which syn2d
export OMP_NUM_THREADS=8
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
tmpdir=/tmp/shotI
mkdir -p $tmpdir
#for dt=0.004 with modeling at 0.0005
scale=1.0
w1=1
smooth=3
#smooth=0
fmute file_shot=shots/shotsmonPz1100_0_rp.su file_out=p0plus.su above=-1 shift=-10 verbose=1 check=1 hw=4
suwind itmax=1023 < p0plus.su | \
suwind key=gx min=-3000000 max=3000000 | \
sushw key=fldr a=1 > p0plussx.su
#~/bin/marchenko file_shot=../shots/refl_rp.su file_tinv=p0plussx.su nshots=601 file_green=pgreen_1.su verbose=1 tap=0 ntap=0 reci=0 niter=1 hw=8 shift=7 smooth=$smooth w=1 file_gplus=Gplus0.su file_gmin=Gmin0.su file_f1plus=f1plus0_1.su file_f1min=f1min0_1.su file_pplus=Pplus0_1.su
#
#~/bin/marchenko file_shot=../shots/refl_rp.su file_tinv=p0plussx.su nshots=601 file_green=pgreen_4.su verbose=1 tap=0 ntap=0 reci=0 niter=4 hw=8 shift=7 smooth=$smooth w=1 file_gplus=Gplus0.su file_gmin=Gmin0.su file_f1plus=f1plus0_4.su file_f1min=f1min0_4.su file_pplus=Pplus0_4.su
#for backpropagating pplus in marchenko scheme must be written to file
~/bin/marchenko file_shot=../shots/refl_rp.su file_tinv=p0plussx.su nshots=601 file_green=pgreen.su verbose=1 tap=0 ntap=10 niter=15 hw=8 shift=7 smooth=$smooth file_gplus=Gplus0.su file_gmin=Gmin0.su file_f1plus=f1plus0.su file_f1min=f1min0.su file_pplus=Pplus0.su
exit;
marchenko/demo/model.scr 0 → 100755
+ 212
− 0
View file @ d4d85049
#!/bin/bash
export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
dx=2.5
dt=0.0005
makemod sizex=6000 sizez=2000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
orig=-3000,0 file_base=syncl.su verbose=2 \
intt=def x=-3000,500,3000 z=195,195,195 poly=1 cp=1950 ro=3700 \
intt=def x=-3000,3000 z=600,600 poly=0 cp=2050 ro=1750 \
intt=def x=-3000,3000 z=680,680 poly=0 cp=2150 ro=2220 \
intt=def x=-3000,3000 z=780,780 poly=0 cp=2230 ro=1700 \
intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=520,580,680,840,680,600,500 poly=2 cp=2400 ro=2800 \
intt=def x=-3000,0,3000 z=1110,1110,1110 poly=0 cp=2300 ro=1950 \
intt=def x=-3000,3000 z=1180,1180 poly=0 cp=2480 ro=1820 \
intt=def x=-3000,0,3000 z=1290,1290,1370 poly=0 cp=2600 ro=2000 \
intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2720 ro=2050 \
intt=def x=-3000,3000 z=1480,1480 poly=0 cp=2800 ro=1850
#intt=diffr x=-2000,-1000,0,1000,2000 z=1800,1800,1800,1800,1800 cp=0,0,0,0,0 ro=5000,5000,5000,5000,5000
makemod sizex=6000 sizez=2000 dx=10 dz=5 cp0=1900 ro0=1200 \
orig=-3000,0 file_base=syncl_migr.su verbose=2 \
intt=def x=-3000,500,3000 z=195,195,195 poly=1 cp=1950 ro=3700 \
intt=def x=-3000,3000 z=600,600 poly=0 cp=2050 ro=1750 \
intt=def x=-3000,3000 z=680,680 poly=0 cp=2150 ro=2220 \
intt=def x=-3000,3000 z=780,780 poly=0 cp=2230 ro=1700 \
intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=520,580,680,840,680,600,500 poly=2 cp=2400 ro=2800 \
intt=def x=-3000,0,3000 z=1110,1110,1110 poly=0 cp=2300 ro=1950 \
intt=def x=-3000,3000 z=1180,1180 poly=0 cp=2480 ro=1820 \
intt=def x=-3000,0,3000 z=1290,1290,1370 poly=0 cp=2600 ro=2000 \
intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2720 ro=2050 \
intt=def x=-3000,3000 z=1480,1480 poly=0 cp=2800 ro=1850
exit
#example FD modeling with model defined above
makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3 scale=0
export OMP_NUM_THREADS=1
makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
~/bin/fdelmodc \
file_cp=syncl_cp.su ischeme=1 iorder=4 \
file_den=syncl_ro.su \
file_src=wave.su \
file_rcv=shot_fd.su \
src_type=1 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0010 \
rec_delay=0.1 \
verbose=2 \
tmod=4.195 \
dxrcv=10.0 \
xrcv1=-2250 xrcv2=2250 \
zrcv1=0 zrcv2=0 \
xsrc=0 zsrc=0 \
ntaper=400 \
tsnap1=3.1 tsnap2=2.5 dtsnap=0.1 \
left=4 right=4 top=4 bottom=4
exit
makemod sizex=6000 sizez=2000 dx=$dx dz=$dx cp0=1900 ro0=1200 \
orig=-3000,-1000 file_base=hom.su
~/bin/fdelmodc \
file_cp=hom_cp.su ischeme=1 iorder=4 \
file_den=hom_ro.su \
file_src=wave.su \
file_rcv=shot_hom_fd.su \
src_type=1 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0010 \
rec_delay=0.1 \
verbose=2 \
tmod=4.195 \
dxrcv=10.0 \
xrcv1=-2250 xrcv2=2250 \
zrcv1=0 zrcv2=0 \
xsrc=0 zsrc=0 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
sudiff shot_fd_rp.su shot_hom_fd_rp.su > shot_rp.su
migr file_shot=shot_rp.su file_vel=syncl_migr_cp.su imc=0 file_image=image.su verbose=1
exit
makemod sizex=6000 sizez=4000 dx=$dx dz=$dx cp0=$cp cs0=$cs ro0=$rho \
orig=-3000,-1000 file_base=synclTop.su \
intt=def x=-3000,0,3000 z=200,200,200 poly=0 cp=1800 ro=5000 \
intt=def x=-3000,-2000,-1000,-400,0,200,900,1800,3000 z=520,520,560,670,950,790,600,520,500 poly=2 cp=2300 ro=1800 \
~/bin/fdelmodc \
file_cp=synclTop_cp.su ischeme=1 iorder=4 \
file_den=synclTop_ro.su \
file_src=wave.su \
file_rcv=p0.su \
src_type=7 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0010 \
rec_delay=0.1 \
verbose=2 \
tmod=2.100 \
dxrcv=10.0 \
xrcv1=-2250 xrcv2=2250 \
zrcv1=0 zrcv2=0 \
xsrc=0 zsrc=1100 \
ntaper=300 \
left=4 right=4 top=4 bottom=4 &
~/bin/fdelmodc \
file_cp=synclTop_cp.su ischeme=1 iorder=4 \
file_den=synclTop_ro.su \
file_src=wave.su \
file_rcv=shot_top.su \
src_type=7 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0010 \
rec_delay=0.1 \
verbose=2 \
tmod=4.195 \
dxrcv=10.0 \
xrcv1=-2250 xrcv2=2250 \
zrcv1=0 zrcv2=0 \
xsrc=0 zsrc=0 \
ntaper=300 \
left=4 right=4 top=4 bottom=4 &
makemod sizex=6000 sizez=4000 dx=$dx dz=$dx cp0=2300 ro0=1800 \
orig=-3000,-1000 file_base=synclBot.su \
intt=def x=-3000,0,3000 z=1310,1310,1310 poly=0 cp=2450 ro=1950 \
intt=def x=-3000,3000 z=1380,1380 poly=0 cp=2460 ro=1820 \
intt=def x=-3000,0,3000 z=1490,1490,1570 poly=0 cp=2470 ro=2100 \
intt=def x=-3000,3000 z=1580,1580 poly=0 cp=2480 ro=2000 \
intt=def x=-3000,3000 z=1680,1680 poly=0 cp=2490 ro=1850
~/bin/fdelmodc \
file_cp=synclBot_cp.su ischeme=1 iorder=4 \
file_den=synclBot_ro.su \
file_src=wave.su \
file_rcv=pRef.su \
src_type=7 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0010 \
rec_delay=0.1 \
verbose=2 \
tmod=2.100 \
dxrcv=10.0 \
xrcv1=-2250 xrcv2=2250 \
zrcv1=1100 zrcv2=1100 \
xsrc=0 zsrc=1100 \
ntaper=300 \
left=4 right=4 top=4 bottom=4 &
makemod sizex=6000 sizez=4000 dx=$dx dz=$dx cp0=2300 ro0=1800 \
orig=-3000,-1000 file_base=synclBotHom.su
~/bin/fdelmodc \
file_cp=synclBotHom_cp.su ischeme=1 iorder=4 \
file_den=synclBotHom_ro.su \
file_src=wave.su \
file_rcv=pRefHom.su \
src_type=7 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.0010 \
rec_delay=0.1 \
verbose=2 \
tmod=2.100 \
dxrcv=10.0 \
xrcv1=-2250 xrcv2=2250 \
zrcv1=1100 zrcv2=1100 \
xsrc=0 zsrc=1100 \
ntaper=300 \
left=4 right=4 top=4 bottom=4 &
wait
sudiff shot_fd_rp.su shot_hom_fd_rp.su > shot_rp.su
sudiff shot_top_rp.su shot_hom_fd_rp.su > shotTop_rp.su
sudiff pRef_rp.su pRefHom_rp.su > pref_rp.su
marchenko/demo/referenceShot.scr 0 → 100755
+ 44
− 0
View file @ d4d85049
#!/bin/bash
#PBS -N fdelmod
#PBS -q verylong
#PBS -l nodes=1
#PBS -k eo
#PBS -j eo
export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop/Redatum
#makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3
dx=2.5
dt=0.0005
makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
export OMP_NUM_THREADS=8
$HOME/bin/fdelmodc \
file_cp=../syncl_cp.su ischeme=1 iorder=4 \
file_den=../syncl_ro.su \
file_src=wave.su \
file_rcv=virtual_shot_fd_P_zsrc1100.su \
src_type=1 \
src_orient=1 \
src_injectionrate=1 \
rec_type_ud=1 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=0.004 \
rec_delay=0.1 \
verbose=2 \
tmod=2.144 \
dxrcv=10.0 \
xrcv1=-3000 xrcv2=3000 \
zrcv1=0 zrcv2=0 \
xsrc=0 zsrc=1100 \
file_snap=backpropref.su tsnap1=0.1 dtsnap=0.010 tsnap2=2.100 dxsnap=10 dzsnap=10 zsnap1=0 zsnap2=2000 xsnap1=-2250 xsnap2=2250 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
marchenko/demo/remove_direct.scr 0 → 100755
+ 38
− 0
View file @ d4d85049
#!/bin/bash
#PBS -N fdelmod
#PBS -q verylong
#PBS -l nodes=1
#PBS -k eo
#PBS -j eo
export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
dxshot=10
ishot=0
nshots=601
rm shots/refl_rp.su
while (( ishot < nshots ))
do
(( xsrc = -3000 + ${ishot}*${dxshot} ))
(( iishot = ${ishot}*${dxshot}/10 ))
(( tr1 = 601 - ${iishot} ))
(( tr2 = ${tr1} + 600 ))
echo xsrc=$xsrc tr1=$tr1 tr2=$tr2
suwind < direct_rp.su key=tracl min=$tr1 max=$tr2 > direct.su
file_rcv=shots/shots_${xsrc}_rp.su
suwind key=tracl min=1 max=601 < $file_rcv > shotz0.su
sudiff shotz0.su direct.su > refl.su
(( ishot = $ishot + 1))
sushw < refl.su key=fldr a=$ishot | \
suwind itmax=1023 >> shots/refl_rp.su
done
marchenko/demo/shots.scr 0 → 100755
+ 77
− 0
View file @ d4d85049
#!/bin/bash
#PBS -N fdelmod
#PBS -q long
#PBS -l nodes=1
#PBS -k eo
#PBS -j eo
export PATH=$HOME/bin:$HOME/src/OpenSource/bin:$PATH:
which makewave
which makemod
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
dt=0.0005
makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefw.su nt=4096 t0=0.3 scale=0 scfft=1
./model.scr
mkdir -p shots
mkdir -p jobs
dxshot=10
ishot=0
nshots=601
zsrc=0
while (( ishot < nshots ))
do
(( xsrc = -3000 + ${ishot}*${dxshot} ))
echo ishot=$ishot xsrc=$xsrc zsrc=$zsrc
cat << EOF > jobs/pbs_$ishot.job
#!/bin/bash
#
#PBS -q medium
#PBS -N mod_${xsrc}
#PBS -j eo
#PBS -m n
#PBS -l nodes=1
#PBS -V
cd /home/thorbcke/data/Kees/MultElim/ModelImageBackprop
export OMP_NUM_THREADS=1
file_rcv=shots/shots_${xsrc}.su
$HOME/bin/fdelmodc \
file_cp=syncl_cp.su ischeme=1 iorder=4 \
file_den=syncl_ro.su \
file_src=wavefw.su \
file_rcv=\$file_rcv \
src_type=7 \
src_orient=1 \
src_injectionrate=1 \
rec_type_vz=1 \
rec_type_p=1 \
rec_int_vz=2 \
rec_delay=0.3 \
dtrcv=0.004 \
verbose=2 \
tmod=4.394 \
dxrcv=10.0 \
xrcv1=-3000,-3000,-3000 xrcv2=3000,3000,3000 \
zrcv1=0,1000,1600 zrcv2=0,1000,1600 \
xsrc=$xsrc zsrc=$zsrc \
ntaper=400 \
left=4 right=4 top=4 bottom=4
EOF
qsub jobs/pbs_$ishot.job
(( ishot = $ishot + 1))
done
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment