Download - Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann
![Page 1: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/1.jpg)
1
Real-time or full-precision CRS imaging using a cloud computing portal: multi-
offset GPR and SH-wave examples
Z. Heilmann*, H.P. Müller**, G. Satta* and G.P. Deidda****CRS4, Energy and Environment Sector
**ABE-geo, Burgdorf, Germany***University of Cagliari, Department of Civil and Environmental
Engineering, and Architecture
![Page 2: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/2.jpg)
2
Outline
Part 1: Cloud computing portalo Concept & Realizationo GPR example (Real-time imaging)
Part 2: CRS stack imagingo from CMP to CRS stacko 3x1 vs. 1x3 parameter searcho SH-wave data (Full-precision imaging)
Conclusions
![Page 3: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/3.jpg)
3
PART 1:EIAGRID PORTAL
+ GPR-DATA EXAMPLE
![Page 4: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/4.jpg)
Basic idea
4
EIAGRID Portal
![Page 5: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/5.jpg)
5
EIAGRID Portal
Basic idea
1. Web-browser-based interface accessible from the field
2. Real-time processing using distant HPC resources3. Remote collaboration and acquisition controlling
Features
![Page 6: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/6.jpg)
6
EIAGRID Portal
Basic idea
1. Web-browser-based interface accessible from the
field2. On-the-fly processing using a distant HPC resources
3. Remote collaboration and acquisition controlling
Features
1. Project, data and user management 2. Simplistic toolbox for data visualization and
manipulation3. Data-driven imaging method suited for parallel
computing
Components
![Page 7: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/7.jpg)
7
![Page 8: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/8.jpg)
8
Multi-offset GPR data:
![Page 9: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/9.jpg)
9
Multi-offset GPR data:● Aim: monitoring of water content and water conductivity● Target depth: 0 - 5 m● Profile length: 55 m
Instrumentation:RAMAC/GPR CU II with MC4 + 4 unshielded 200 MHz antennas
Geometry:Number of sources: 546Source spacing: 0.1mNumber of receivers: 28Receiver spacing: 0.2 mMaximum offset: 6 m
CMP gather at 10 m
![Page 10: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/10.jpg)
Data visualization tools
![Page 11: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/11.jpg)
Data visualization tools
![Page 12: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/12.jpg)
Data visualization tools
![Page 13: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/13.jpg)
![Page 14: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/14.jpg)
Data visualization tools
cm/µs
cm/µs
cm/µs
MHz
![Page 15: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/15.jpg)
Data visualization tools
MHz
cm/µs
cm/µs
cm/µs
cm/µs
cm
cm
cm
cm
µs
![Page 16: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/16.jpg)
![Page 17: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/17.jpg)
![Page 18: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/18.jpg)
CRS stacking result obtained after 4 minutes using 50 CPU
![Page 19: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/19.jpg)
Time domain imaging
Published in: Perroud, H., and Tygel, M., 2005, Velocity estimation by the common-reflection-surface (CRS) method: Using ground-penetrating radar: Geophysics, 70, 1343–1352.
Results GPR data
![Page 20: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/20.jpg)
20
GPR data
Comparison with published results:
Figure (a) CRS stack section and (b) time migrated section overlain with the migration velocity model obtained from stacking parameters after smoothing and regularization..
Figure taken from Perroud and Tygel 2005. (a) Final stack and (b) RMS velocity sections obtained by the classical NMO method; the black lines in (a) represent the limits of the CRS midpoint aperture for the first-order (inner curves) and second-order (outer curves) parameter searches.
![Page 21: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/21.jpg)
21
PART 2:CRS STACK - REVISITED
+ SH-DATA EXAMPLE
![Page 22: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/22.jpg)
From CMP to CRS stacking:
Figure taken from Perroud and Tygel 2005. NMO velocity analysis for the CMP at position x = 10 m.
22
Conventional CMP-by-CMP velocity analysis:
𝑡𝐶𝑀𝑃2 (h )=𝑡0
2+4h2
𝑣𝑁𝑀𝑂2
![Page 23: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/23.jpg)
CRS stacking operator:Hyperbolic traveltime in CRS gather:
… in CMP gather:
Fig.: Mann et al. 2007
![Page 24: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/24.jpg)
24
Stacking parameter search:
Pragmatic search: 3 x 1 parameter line search in specific
gathers (Mann et al. 1999)+ Cloud = Real-time imaging
One step search:1 x 3 parameter surface search in prestack data (Garabito et al. 2001)
+ Cloud = High-precision imaging
Figs: Mann et al. 2007
![Page 25: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/25.jpg)
Algorithm: For every sample we maximize Coh(, ,) ● S
tacked ZO-S
ection
𝜶
𝑹𝑵𝑰𝑷
𝑹𝑵
Object function: Semblance (T&K 94)
![Page 26: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/26.jpg)
26
Shallow SH-data:
![Page 27: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/27.jpg)
27
Shallow SH-data:
![Page 28: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/28.jpg)
28
Urban SH-wave data:● Aim: aquiclude survey for site remediation● Target depth: 0 - 15 m● 2D line: length 92 m
Source: mini-vibrator, dS=0.5 mReceiver: landstreamer with 47
geophones, dR=0.5 m
![Page 29: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/29.jpg)
29
Some CMP gathers:
CMP 1171 CMP 1267CMP 1069
Offset [dm] Offset [dm] Offset [dm]
Tim
e [m
s]
Tim
e [m
s]
Tim
e [m
s]
![Page 30: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/30.jpg)
30
Stacking results:
CMP stack after CVS analysis (a) versus 3x1 parameter CRS stack (b)
a)
b)
+ higher S/N
- unresolved events
![Page 31: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/31.jpg)
31
Stacking results:
CMP stack after CVS analysis (a), 3x1 parameter search CRS stack (b)
a)
b)
![Page 32: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/32.jpg)
32
Shear wave data
Stacking results:
CMP stack after CVS analysis (a), 1x3 parameter search CRS stack (b)
a)
b)
+ higher S/N
+ improved resolution
![Page 33: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/33.jpg)
Computational cost:
3x1 Parameter
2-Parameter
3-Parameter
0 10 20 30 40 50 60 70 80 90
submission-time [min] on 30 CPU runtime [min] on 30 CPUs
![Page 34: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/34.jpg)
34
1x3 versus 3x1 parameter search
High precision CRS (a) versus real-time CRS stack section (b)
a)
b)
![Page 35: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/35.jpg)
35
Shear wave data
1x3 versus 3x1 parameter search
NMO velocities calculated using and : High precision CRS (a) versus real-time CRS stack section (b)
a)
b)
![Page 36: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/36.jpg)
36
Outlook: Depth imaging
(a) Velocity model from inversion of and , (b) PostSDM
a)
b)
![Page 37: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/37.jpg)
37
Outlook: Depth imaging
(a) Velocity model from inversion of and , (b) PreSDM
a)
b)
![Page 38: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/38.jpg)
38
Outlook: Depth migration
(a) PreSDM (b) PostSDM using Dix inversion of CVS velocity
a)
b)
Top of Molasse
VSP
![Page 39: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/39.jpg)
39
Conclusions:
• Real-time CRS imaging can be applied to optimize crutial acquisition parameters directly in the field: easier use of reflection methods in near-surface.
• Full-precision CRS imaging uses a spatial operator not only for stacking but also for velocity analysis: higher resolution and more stable velocities in case of low fold, strong noise & lateral inhomogeneity.
• A cloud computing portal provides optimum computing power in a location independent way: reduced hardware requirement facilitate the use of data driven methods in near-surface imaging.
![Page 40: Near Surface Geoscience Conference 2014, Athens - Real-time or full‐precision CRS - Zeno Heilmann](https://reader035.vdocument.in/reader035/viewer/2022070520/58f097c11a28ab491a8b45cf/html5/thumbnails/40.jpg)
Thank you for your attention!
40