Commit 6c4b9695 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

[WP] (i20_gremlin_example): provide seismogram preprocessing

- provide tapers for muting
- provide muting operation
- provide arrival time picks
- provide refracted wave travel times and model
- provide plots for all of them
parent 2879e29d
......@@ -47,6 +47,7 @@ include $(patsubst %.main.cfg,%.d,$(wildcard *.main.cfg))
# standard targets
# ----------------
%.pdf: %.ps; ps2pdf $<
%.psp: %.ps; gv $<; /bin/rm -fv $<
# ======================================================================
#
......@@ -106,7 +107,45 @@ expHQ.xxx: ; echo 0.4 > $@
expLQ.xxx: ; echo 0.8 > $@
simulation%.ps: simulation%.$(FILETYPE) exp%.xxx
refractx -ty $(FILETYPE) -d $@/ps -St -0.15,0.3 -Ev \
-Sm 2 -Sr 0.37 -Eg -Se $$(cat $(word 2,$^)) -Sa 15. -Sc 5. \
-Tt "simulated data with gaussian, uncorrelated noise added" $<
-Sm 2 -Sr 0.37 -Eg -Se $$(cat $(word 2,$^)) -Sa 15. -Sc 5. \
-Tt "simulated data with gaussian, uncorrelated noise added" $<
# ======================================================================
# preprocess seismogramm data
# ---------------------------
#
# tapers, picks and model are only provided once but are required for the HQ
# and the LQ version
layer%.tap.xxx: layer.tap.cfg; /bin/cp -vpd $< $@
layer%.ttc.xxx: layer.ttc.cfg; /bin/cp -vpd $< $@
layer%.tt.xxx: layer.tt.cfg; /bin/cp -vpd $< $@
layer%.VP.mod.xxx: layer.VP.mod.cfg; /bin/cp -vpd $< $@
# display arrival time picks
arrival%.ps: simulation%.$(FILETYPE) exp%.xxx layer%.tt.xxx
refractx -ty $(FILETYPE) -d $@/cps -St -0.01,0.15 -Ev \
-Lt 3 -Cc -EP T -ET T -Et T -Fa $(word 3,$^) \
-Sm 2 -Sr 1.60 -Eg -Se $$(cat $(word 2,$^)) -Sa 30. -Sc 2. \
-Tt "simulated data and first arrivals" $<
#
# display travel time polygonial and P-wave velocity model
vpmod%.ps: simulation%.$(FILETYPE) exp%.xxx layer%.ttc.xxx layer%.VP.mod.xxx
refractx -ty $(FILETYPE) -d $@/cps -St -0.01,0.15 -Ev -TM 0.06,20.\
-Lt 3 -Cc -EP T -ET T -Et T -Fp $(word 3,$^) -Fm $(word 4,$^)\
-Sm 2 -Sr 1.60 -Eg -Se $$(cat $(word 2,$^)) -Sa 30. -Sc 2. \
-Tt "simulated data and travel time curve and model for refracted P-waves" $<
#
# mute seismograms by application of taper
muted%.$(FILETYPE): simulation%.$(FILETYPE) layer%.tap.xxx
tidofi -v -o -type $(FILETYPE) -Type $(FILETYPE) \
$@ $< poft:$(word 2,$^)
# display muted seismograms
muted%.ps: muted%.$(FILETYPE) exp%.xxx layer%.tap.xxx
refractx -ty $(FILETYPE) -d $@/cps -St -0.15,0.5 \
-Lt 3 -Cc -EP T -ET T -Et T -Ft $(word 3,$^) \
-Sm 2 -Sr 0.37 -Eg -Se $$(cat $(word 2,$^)) -Sa 15. -Sc 5. \
-Tt "muted data with taper definition" $<
# ----- END OF Makefile -----
REFRACT V4.14 REFRACTion seismics - data interpretation
earth model containing 3 layers
layer depth [m] thickness [m] velocity [m/s] slowness [s/m]
1 0.000 3.254 384.726 0.0025992531
2 3.254 7.209 861.226 0.0011611353
hs: 10.463 1602.216 0.0006241354
offset dependent time domain taper
REFRACT V4.14 REFRACTion seismics - data interpretation
taper set 1 with 4 picks
offset [m] time [s]
-0.108450 -0.253573E-02
20.3292 0.277341E-01
30.9524 0.400942E-01
72.8547 0.642847E-01
taper set 2 with 6 picks
offset [m] time [s]
0.115263 0.414383E-03
6.60295 0.129090E-01
19.8571 0.339932E-01
30.3622 0.468478E-01
41.9296 0.560666E-01
72.7367 0.735445E-01
taper set 3 with 5 picks
offset [m] time [s]
-0.562935 0.190927
5.92898 0.243020
18.7948 0.311850
34.8475 0.383540
72.9908 0.528457
taper set 4 with 3 picks
offset [m] time [s]
-0.371309 0.276536
31.5425 0.448719
72.9941 0.595710
REFRACT V4.14 REFRACTion seismics - data interpretation
72
0.0010000 0.0007157
0.0020000 0.0012662
0.0030000 0.0018167
0.0040000 0.0018167
0.0050000 0.0030554
0.0060000 0.0072394
0.0070000 0.0105425
0.0080000 0.0126070
0.0090000 0.0156349
0.0100000 0.0181122
0.0110000 0.0205896
0.0120000 0.0237551
0.0130000 0.0265077
0.0140000 0.0280217
0.0150000 0.0293980
0.0160000 0.0298109
0.0170000 0.0309119
0.0180000 0.0316001
0.0190000 0.0327011
0.0200000 0.0347656
0.0210000 0.0387569
0.0220000 0.0408214
0.0230000 0.0428858
0.0240000 0.0439869
0.0250000 0.0456385
0.0260000 0.0472900
0.0270000 0.0481158
0.0280000 0.0494921
0.0290000 0.0507308
0.0300000 0.0518319
0.0310000 0.0526576
0.0320000 0.0536211
0.0330000 0.0540339
0.0350000 0.0559608
0.0350000 0.0548597
0.0360000 0.0563737
0.0370000 0.0569242
0.0380000 0.0577500
0.0390000 0.0581629
0.0400000 0.0591263
0.0420000 0.0596768
0.0430000 0.0611908
0.0430000 0.0602273
0.0440000 0.0609155
0.0450000 0.0613284
0.0460000 0.0628423
0.0480000 0.0629800
0.0480000 0.0628423
0.0490000 0.0616037
0.0500000 0.0620165
0.0510000 0.0644939
0.0520000 0.0649068
0.0530000 0.0655950
0.0540000 0.0654573
0.0550000 0.0646315
0.0560000 0.0664207
0.0570000 0.0662831
0.0580000 0.0682099
0.0590000 0.0684852
0.0600000 0.0699992
0.0610000 0.0697239
0.0630000 0.0716507
0.0630000 0.0715131
0.0640000 0.0717883
0.0650000 0.0733023
0.0660000 0.0734399
0.0670000 0.0738528
0.0680000 0.0742657
0.0690000 0.0742657
0.0700000 0.0742657
0.0710000 0.0750915
0.0720000 0.0750915
REFRACT V4.14 REFRACTion seismics - data interpretation
3
0.0105236 0.0273534
0.0286858 0.0484422
0.0737375 0.0765606
Supports Markdown
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