Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
Code to reproduce the example in section `Plane-Wave Marchenko algorithm` of the manual.
Description of files:
1) ./model.scr computes the gridded velocity/density model
2) ./shots.scr computes 601 shots using slurm arrays: ~100 s runtime for each shot
3) ./direct.scr compute 1 shot that contains direct wave only
4) ./remove_direct.scr removes the direct wave from the data created in 2 and places all shots in one file shots/refl_rp.su
Figures of the model and the middle shot are generated by
./epsModel.scr model
5) IniFocus.scr compute plane wave responses at 600 m depth with angles -5 0 and 5 degrees.
To compute the plane wave respone with the sources all startin at t=0 on a slanted line in the model run
./Ini2Focus.scr
To generate the 12 snapshot figures
./epsModel.scr snapshots
6) To generate the Marchenko results for plane wave with angle=0 (horizontal)
./marchenkoPlaneIter.scr plane0
./epsModel.scr plane0
[ 7) Marchenko with and without time-shifts in a one-D model: requires shift to build in again in marchenko code: does not generate
correct results at the moment
cd shifts
compute 1D model, f_d and marchenko with and without shift
./model.scr
./initialFocus.scr
./marchenko.scr
./eps.scr
]
2) itertions.scr computes the intermediate results of the multiple attenutation scheme and produces all output files that are used in the manuscript.
- runtime on 4 cores is
3) epsPrimaries.scr selected output from step 2) are converted to .eps pictures that are used in the Figures to explain the method.
To reproduce the postscript files of the manuscript SU postscript plotting programs are required.
-3) epsModel.scr to generate the postscript files for the numerical model
optional scripts not needed to reproduce the figures:
+) primaries.scr computes the internal multiple attenuated (middle) shot record.
- runtime on 4 cores is ~500 s.
+) primariesPlane.scr: computes the internal moval scheme for plane-waves (see Meles 2020)
+) clean: remove all produced files and start with a clean directory
--------------------------
To reproduce the Figures in the Manuscript:
--------------------------
* Figure 2: Model + Initial wavefield
==> run model.scr to generate the data .su files: this will take 4-5 minutes. The files generate are:
- hom_cp.su, hom_ro.su
- model10_cp.su, model10_ro.su
- shot5_fd_rp.su
- shot5_hom_fd_rp.su
- shot5_rp.su
- wavefw.su
==> run './epsPrimaries.scr Figure2' to generate the postscript files of Figure 2
model_cp_line.eps => Figure 2a
model_ro_line.eps => Figure 2b
shotx0_rp.eps => Figure 2c
It also produces two extra pictures of the wavelet used in the FD modelling:
wavefw_freq.eps
wavefw.eps
--------------------------
* Figure 3: First Iteration
==> run './iterations.scr Figure34910' to compute the intermediate results of the first iterations of the Marchenko Primaries algorithm.
This will take 15 seconds. The generated files are:
- M0_276000.su
- Mi_2760##.su
- k1min_2760##.su
- v1plus_2760##.su
- iter_2760##.su (not used)
- pred_rr_276.su (not used)
- DDshot_450.su (not used): selected shot record convolved with file_src=wave.su
where ## ranges from 01 to 34
To generate the postscript files for Figure 3:
==> run './epsPrimaries.scr Figure3'
This will produce the following files:
shotx0_rp.eps => Figure 2c == Figure 3a
M0_276000_flip.eps => Figure 3b
fconvN0fulltime.eps => Figure 3c
fconvN0flip.eps => Figure 3d
Mi_276001.eps => Figure 3e
--------------------------
* Figure 4 second iteration
To generate the postscript files for Figure 4:
==> run './epsPrimaries.scr Figure4'
This will produce the following files:
fconvN1fulltime.eps => Figure 4c
fconvN1flip.eps => Figure 4d
Mi_276002.eps => Figure 4e
The window time function in Figure 5 is not reproduced.
--------------------------
* Figure 6 v1plus and convergence
==> run './iterations.scr Figure6' to compute the marchenko results with ii=200
To generate the postscript files for Figure 6:
==> run './epsPrimaries.scr Figure6'
This will produce the following files:
v1plus_200001.eps => Figure 6a
v1plus_max.eps => Figure 6b
k1min_200030.eps => Figure 6b
--------------------------
* Figure 8 To compute the convergence for a strong contrast medium:
cd strongContrast
==> run ./model.scr
==> run ./iterations.scr Figure8
To generate the postscript files for Figure 8:
==> run './epsPrimaries.scr Figure8'
This will produce the following files:
v1plusStrong_max.eps => Figure 8
Don't forget to go back to the main directory with the regular contrast results
cd ../
--------------------------
* Figure 9 iterations M_i
To generate the postscript files for Figure 9:
==> run './epsPrimaries.scr Figure9'
This will produce the following files:
Mi_276002.eps => Figure 9b
Mi_276004.eps => Figure 9d
Mi_276012.eps => Figure 9f
Mi_276020.eps => Figure 9h
Mi_276001.eps => Figure 9a
Mi_276003.eps => Figure 9c
Mi_276011.eps => Figure 9e
Mi_276019.eps => Figure 9g
--------------------------
* Figure 10 iterations M_i and k_1^-
To generate the postscript files for Figure 10:
==> run './epsPrimaries.scr Figure10'
This will produce the following files:
Mi_276002flip.eps => Figure 10a
k1min_276002.eps => Figure 10b
Mi_276004flip.eps => Figure 10c
k1min_276004.eps => Figure 10d
Mi_276010flip.eps => Figure 10e
k1min_276010.eps => Figure 10f
Mi_276020flip.eps => Figure 10g
k1min_276020.eps => Figure 10h
--------------------------
* Figure 11 iterations k_1^- for different ii 246:316:10
To generate the data
==> run ./iterations.scr Figure11
this will take ~2 minutes and generate a lot of files
To generate the postscript files for Figure 11:
==> run './epsPrimaries.scr Figure11'
This will produce the following files:
k1min_246032.eps => Figure 11a
k1min_256032.eps => Figure 11b
k1min_266032.eps => Figure 11c
k1min_276032.eps => Figure 11d
k1min_286032.eps => Figure 11e
k1min_296032.eps => Figure 11f
k1min_306032.eps => Figure 11g
k1min_316032.eps => Figure 11h
--------------------------
* Figure 13 iterations M_i and k_1^- for ii-276 T-MME scheme
To generate the data
==> run ./iterations.scr Figure13
this will take ~15 seconds
****NOTE this will overwrite the results of the MME-scheme !
To generate the postscript files for Figure 13:
==> run './epsPrimaries.scr Figure13'
This will produce the following files:
Mi_276002T.eps => Figure 13a
Mi_276004T.eps => Figure 13b
Mi_276001T.eps => Figure 13c
Mi_276003T.eps => Figure 13d
k1min_276002T.eps => Figure 13e
k1min_276004T.eps => Figure 13f
k1min_276010T.eps => Figure 13g
k1min_276020T.eps => Figure 13h