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

small updates and changes

parent 0db27382
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
fdelmodc/demo/back_injrate.scr 0 → 100755
+ 107
− 0
View file @ 7fc510d1
#!/bin/bash
export PATH=:$HOME/src/OpenSource/bin:$HOME/bin64:$PATH:
dx=1.0
dt=0.0002
makemod sizex=700 sizez=700 dx=$dx dz=$dx cp0=1500 ro0=1000 \
orig=-350,-350 file_base=water.su verbose=2 \
intt=def x=-150,150 z=-150,-150 poly=1 cp=1950 ro=3700 \
intt=def x=-150,150 z=150,150 poly=1 cp=1500 ro=1000
makewave w=fw fmin=0 flef=5 frig=60 fmax=65 dt=$dt file_out=wavefw.su nt=4096 t0=0.4 scale=0
makewave fp=20 dt=$dt file_out=wave.su nt=4096 t0=0.1 scale=1
#circular receivers
../fdelmodc \
file_cp=water_cp.su ischeme=1 iorder=4 \
file_den=water_ro.su \
file_src=wave.su \
file_rcv=inj_rate.su \
src_type=1 \
src_injectionrate=1 \
rec_type_vz=1 \
rec_type_p=1 \
rec_int_p=0 \
dtrcv=$dt \
rec_delay=0.0 \
verbose=4 \
tmod=0.8000 \
xrcv1=-350,-350,-350,350 xrcv2=350,350,-350,350 zrcv1=-350,350,-350,-350 zrcv2=-350,350,350,350 \
dxrcv=$dx,$dx,0,0 dzrcv=0,0,$dx,$dx \
xsrc=0 zsrc=0 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
# oxrcv=0 ozrcv=0 dphi=1.0 rrcv=300 \
#middle trace of snapshot number 11 at t=0
#suwind key=fldr min=11 max=11 < snapinj_rate_sp.su | suwind key=tracf min=350 max=350 | suxgraph
#suwind key=fldr min=11 max=11 < snapinj_rate_sp.su | suwind key=tracf min=350 max=350 | suxgraph x1beg=-100 x1end=100 style=normal
../fdelmodc \
file_cp=water_cp.su ischeme=1 iorder=4 \
file_den=water_ro.su \
file_src=inj_rate_rp.su \
file_rcv=back_inj_rate.su \
file_snap=snapinj_rate.su \
grid_dir=1 \
src_type=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_p=0 \
dtrcv=$dt \
rec_delay=0.0 \
verbose=4 \
tmod=1.1000 \
xrcv1=-300 xrcv2=300 \
zrcv1=0 zrcv2=0 dxrcv=1 \
tsnap1=0.000 tsnap2=0.90 dtsnap=0.01 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
exit
../fdelmodc \
file_cp=water_cp.su ischeme=1 iorder=4 \
file_den=water_ro.su \
file_src=wave.su \
file_rcv=inj.su \
src_type=1 \
src_injectionrate=0 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_p=0 \
dtrcv=0.0004 \
rec_delay=0.0 \
verbose=4 \
tmod=0.4000 \
oxrcv=0 ozrcv=0 dphi=1.0 rrcv=300 \
xsrc=0 zsrc=0 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
../fdelmodc \
file_cp=water_cp.su ischeme=1 iorder=4 \
file_den=water_ro.su \
file_src=inj_rp.su \
file_rcv=back_inj.su \
file_snap=snapinj.su \
grid_dir=1 \
src_type=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_p=0 \
dtrcv=0.0004 \
rec_delay=0.0 \
verbose=4 \
tmod=0.5000 \
xrcv1=-300 xrcv2=300 \
zrcv1=0 zrcv2=0 dxrcv=1 \
tsnap1=0.200 tsnap2=0.50 dtsnap=0.01 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
fdelmodc/demo/back_injrate_planes.scr 0 → 100755
+ 138
− 0
View file @ 7fc510d1
#!/bin/bash
export PATH=:$HOME/src/OpenSource/bin:$HOME/bin64:$PATH:
dx=0.5
dt=0.0001
fast=""
dx=2.5
dt=0.00050
fast="fast"
makemod sizex=6000 sizez=1250 dx=$dx dz=$dx cp0=1500 ro0=1000 \
orig=-3000,0 file_base=tutodelfast.su verbose=2 \
intt=def x=-3000,-2000,-1000,0,1000,2000,3000 z=240,130,250,300,350,380,320 poly=2 cp=1950 ro=4500 grad=0 \
intt=def x=-3000,-2200,-1500,0,1300,2100,3000 z=620,640,590,600,740,700,600 poly=2 cp=2000 ro=1400 grad=0 \
intt=def x=-3000,-1800,0,2200,3000 z=920,1000,900,1000,1010 poly=2 cp=2300 ro=1600 grad=0
# intt=def x=-3000,-1800,0,2200,3000 z=920,1000,900,1000,1010 poly=2 cp=1500 ro=1000 grad=0
#makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefwfast.su nt=4096 t0=0.4 scale=0 scfft=1
export OMP_NUM_THREADS=8
makewave fp=30 dt=$dt file_out=wave${fast}.su nt=4096 t0=0.1 scale=1
file_cp=tutodel${fast}_cp.su
file_ro=tutodel${fast}_ro.su
fdelmodc \
file_cp=$file_cp ischeme=1 iorder=4 \
file_den=$file_ro \
file_src=wave${fast}.su \
file_rcv=inj_rate_plane.su \
src_type=1 \
src_injectionrate=1 \
rec_type_vz=1 \
rec_type_p=1 \
rec_int_vz=2 \
dtrcv=$dt \
rec_delay=0.1 \
verbose=4 \
tmod=3.4000 \
xrcv1=-3000,-3000 xrcv2=3000,3000 zrcv1=0,1200 zrcv2=0,1200 \
dxrcv=2.5,2.5 dzrcv=0,0 \
xsrc=0 zsrc=800 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
# xrcv1=-3000,-3000,-3000,3000 xrcv2=3000,3000,-3000,3000 zrcv1=0,1200,0,0 zrcv2=0,1200,1200,1200 \
# dxrcv=2.5,2.5,0,0 dzrcv=0,0,2.5,2.5 \
# 1500/2300
scale=0.65217391304347826086
suwind key=tracl min=1 max=2401 < inj_rate_plane_rp.su > inj_rate_plane0_rp.su
suwind key=tracl min=2402 max=4802 < inj_rate_plane_rp.su | sugain scale=$scale >> inj_rate_plane0_rp.su
scale=-1.0
suwind key=tracl min=1 max=2401 < inj_rate_plane_rvz.su | sutaper tr1=100 tr2=100 ntr=2401 > inj_rate_plane0_rvz.su
suwind key=tracl min=2402 max=4802 < inj_rate_plane_rvz.su | sutaper tr1=100 tr2=100 ntr=2401 | sugain scale=$scale >> inj_rate_plane0_rvz.su
#good results with scaled Vz src_type=1 src_injectionrate=0 and Vz snapshots
# suxmovie < snapinj_rate_plane1vzrate0_sp.su loop=1 clip=30
#good results with scaled Vz src_type=1 src_injectionrate=1 and Vz snapshots
# suxmovie < snapinj_rate_plane1vzrate1_sp.su loop=1 clip=4000
# at 3.2000 seconds the focus is at t=0
file_snap=snapinj_rate_plane0_rvzrate1
fdelmodc \
file_cp=$file_cp ischeme=1 iorder=4 \
file_den=$file_ro \
file_src=inj_rate_plane0_rvz.su \
file_rcv=back_inj_rate_plane.su \
file_snap=$file_snap.su \
grid_dir=1 \
src_type=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_p=0 \
dtrcv=$dt \
rec_delay=0.0 \
verbose=4 \
tmod=4.0000 \
tsnap1=2.000 tsnap2=4.00 dtsnap=0.05 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
# tsnap1=3.200 tsnap2=3.50 dtsnap=0.004 \
# tsnap1=3.290 tsnap2=3.31 dtsnap=0.0005 \
suop2 tutodelfast_cp.su SrcRecPositions.su op=sum w1=1 w2=2000 | \
supsimage hbox=4 wbox=6 labelsize=10 \
x1beg=0 x1end=1250.0 legend=0 threecolor=1 \
d1s=1.0 d2s=1.0 \
wrgb=1.0,.5,0 grgb=0,.7,1.0 brgb=0,1.0,0 \
bps=24 bclip=2400 wclip=1500 \
n1tic=5 x2beg=-3000 f2num=-3000 d2num=1000 x2end=3000 > tutodelBackacq_cp.eps
sumax < inj_rate_plane0_rvz.su mode=abs outpar=nep
clip=`cat nep | awk '{print $1/3}'`
sugain epow=1.2 < inj_rate_plane0_rvz.su | \
supsimage hbox=4 wbox=8 labelsize=10 linewidth=0.0 \
n1tic=2 d2=2.5 f1=0.0 x1beg=0 x1end=3.004 \
f2=0 f2num=0 d2num=1000 clip=$clip > inj_rate_rvz.eps
sumax < inj_rate_plane0_rp.su mode=abs outpar=nep
clip=`cat nep | awk '{print $1/3}'`
sugain epow=1.2 < inj_rate_plane0_rp.su | \
supsimage hbox=4 wbox=8 labelsize=10 linewidth=0.0 \
n1tic=2 d2=2.5 f1=0.0 x1beg=0 x1end=3.004 \
f2=0 f2num=0 d2num=1000 clip=$clip > inj_rate_rp.eps
sumax < ${file_snap}_sp.su mode=abs outpar=nep
clip=`cat nep | awk '{print $1/10}'`
echo $clip
for fldr in 27 32 22 17 12 7
do
suwind key=fldr min=$fldr max=$fldr < ${file_snap}_sp.su | \
supsimage hbox=4 wbox=6 labelsize=10 \
x1beg=0 x1end=1200.0 clip=$clip \
n1tic=5 x2beg=-1000 f2num=-1000 d2num=500 x2end=1000 > snapinj_rate_$fldr.eps
done
suwind key=fldr min=25 max=25 < ${file_snap}_sp.su | \
suwind key=gx min=0 max=0 | \
supsgraph hbox=6 wbox=2 labelsize=10 style=seismic \
x1beg=0 d1=2.5 f1=0 > snapinj_rate_27_trace0_t0.eps
suwind key=fldr min=25 max=25 < ${file_snap}_sp.su | \
suwind key=gx itmin=321 itmax=321 | sustrip > trace.bin
suaddhead < trace.bin n1=2401 dt=25 | supsgraph hbox=2 wbox=6 labelsize=10 \
f1=-3000 d1=$dx f1num=-3000 d1num=1000 style=normal > snapinj_rate_27_z800_t0.eps
fdelmodc/demo/back_injrate_planes_1D.scr 0 → 100755
+ 128
− 0
View file @ 7fc510d1
#!/bin/bash
export PATH=:$HOME/src/OpenSource/bin:$HOME/bin64:$PATH:
dx=0.5
dt=0.0001
fast=""
dx=2.5
dt=0.00050
fast="fast"
makemod sizex=6000 sizez=1250 dx=$dx dz=$dx cp0=1500 ro0=1000 \
orig=-3000,0 file_base=tutodelfast.su verbose=2 \
intt=def x=-3000,3000 z=240,240 cp=1950 ro=4500 grad=0 \
intt=def x=-3000,3000 z=620,620 cp=1500 ro=1000 grad=0
makewave w=fw fmin=0 flef=5 frig=80 fmax=100 dt=$dt file_out=wavefwfast.su nt=4096 t0=0.4 scale=0 scfft=1
export OMP_NUM_THREADS=8
makewave fp=30 dt=$dt file_out=wave${fast}.su nt=4096 t0=0.1 scale=1
file_cp=tutodel${fast}_cp.su
file_ro=tutodel${fast}_ro.su
fdelmodc \
file_cp=$file_cp ischeme=1 iorder=4 \
file_den=$file_ro \
file_src=wave${fast}.su \
file_rcv=inj_rate_plane.su \
src_type=1 \
src_injectionrate=1 \
rec_type_vz=1 \
rec_type_p=1 \
rec_int_vz=2 \
rec_int_vx=2 \
dtrcv=$dt \
rec_delay=0.0 \
verbose=4 \
tmod=4.4000 \
xrcv1=-3000,-3000 xrcv2=3000,3000 zrcv1=0,1200 zrcv2=0,1200 \
dxrcv=2.5,2.5 dzrcv=0,0 \
xsrc=0 zsrc=800 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
# xrcv1=-3000,-3000,-3000,3000 xrcv2=3000,3000,-3000,3000 zrcv1=0,1200,0,0 zrcv2=0,1200,1200,1200 \
# dxrcv=2.5,2.5,0,0 dzrcv=0,0,2.5,2.5 \
scale=1
scale=0.40760869
suwind key=tracl min=1 max=2401 < inj_rate_plane_rp.su > inj_rate_plane0_rp.su
suwind key=tracl min=2402 max=4802 < inj_rate_plane_rp.su | sugain scale=$scale >> inj_rate_plane0_rp.su
scale=-1.0
suwind key=tracl min=1 max=2401 < inj_rate_plane_rvz.su > inj_rate_planez0_rvz.su
suwind key=tracl min=2402 max=4802 < inj_rate_plane_rvz.su | sugain scale=$scale >> inj_rate_planez0_rvz.su
#suwind key=tracl min=2402 max=4802 < inj_rate_plane_rvz.su | sugain scale=$scale > inj_rate_planez1200_rvz.su
#good results with scaled Vz src_type=1 src_injectionrate=0 and Vz snapshots
# suxmovie < snapinj_rate_plane1vzrate0_sp.su loop=1 clip=30
#good results with scaled Vz src_type=1 src_injectionrate=1 and Vz snapshots
# suxmovie < snapinj_rate_plane1vzrate1_sp.su loop=1 clip=4000
fdelmodc \
file_cp=$file_cp ischeme=1 iorder=4 \
file_den=$file_ro \
file_src=inj_rate_planez0_rvz.su \
file_rcv=back_inj_rate_planez0.su \
file_snap=snapinj_rate1_planevz0.su \
grid_dir=1 \
src_type=1 \
src_injectionrate=1 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_p=0 \
dtrcv=$dt \
rec_delay=0.0 \
verbose=4 \
tmod=5.0000 \
tsnap1=4.200 tsnap2=4.50 dtsnap=0.004 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
exit
fdelmodc \
file_cp=$file_cp ischeme=1 iorder=4 \
file_den=$file_ro \
file_src=inj_rate_planez1200_rvz.su \
file_rcv=back_inj_rate_planez1200.su \
file_snap=snapinj_rate_planevz1200.su \
grid_dir=1 \
src_type=1 \
src_injectionrate=0 \
rec_type_vz=0 \
rec_type_p=1 \
rec_int_p=0 \
dtrcv=$dt \
rec_delay=0.0 \
verbose=4 \
tmod=4.0000 \
tsnap1=3.200 tsnap2=3.50 dtsnap=0.004 \
ntaper=400 \
left=4 right=4 top=4 bottom=4
suop2 snapinj_rate_planevz0_sp.su snapinj_rate_planevz1200_sp.su op=sum w1=1 w2=1 > snapinj_rate_planevzsum_sp.su
exit
sumax < inj_rate_plane0_rvz.su mode=abs outpar=nep
clip=`cat nep | awk '{print $1/3}'`
echo $clip
sugain epow=1.2 < inj_rate_plane0_rvz.su | \
supsimage hbox=4 wbox=8 labelsize=10 linewidth=0.0 \
n1tic=2 d2=2.5 f1=0.0 x1beg=0 x1end=3.004 \
f2=0 f2num=0 d2num=1000 clip=$clip > inj_rate_rvz.eps
for fldr in 27 32 22 17 12 7
do
suwind key=fldr min=$fldr max=$fldr < snapinj_rate_plane1vzrate0_sp.su | \
supsimage hbox=4 wbox=6 labelsize=10 \
x1beg=0 x1end=1200.0 clip=30 \
n1tic=5 x2beg=-3000 f2num=-3000 d2num=1000 x2end=3000 > snapinj_rate_$fldr.eps
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