Chapter 3: Wave-mode separation in the continuous wavelet
domain using combined translational and rotational data
Vertical component
P waves
Current state of the art
2
Data components:
Hydrophone
Vertical
Radial
Transverse
“Coherent noise” /
“Source generated noise”
S-waves
Wave modes:
Surface waves
“ground roll”
“Vz noise”
P-waves
Wave-mode selectivity using translational + rotational data
3
Data components:
Hydrophone
Vertical
Radial
Transverse
Yaw
Roll
Pitch
…
Wave modes:
P-waves
S-waves
Rayleigh
Love
….
Agenda
4
Agenda
5
Agenda
6
Agenda
7
Coherent Noise
8
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Ground roll / Rayleigh / Surface waves
9
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Ground roll / Rayleigh / Surface waves
10
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Data courtesy of Shell
Ground roll / Rayleigh / Surface waves
11
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Data courtesy of Shell
Ground roll suppression: Mute
12
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Data courtesy of Shell
Ground roll suppression with F-K
13
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Receivers spaced at 2m intervals
Courtesy of Gustavo Alves
Ground roll suppression with F-K
14
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Receivers spaced at 8m intervals
Courtesy of Gustavo Alves
Ground roll suppression with polarization analysis
15
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Ground roll suppression with polarization analysis
16
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Not applicable to single component data
Ground roll suppression with polarization analysis
17
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Can operate on spatially aliased data (3D surveys)
Insensitive to statics
Ground roll suppression with polarization analysis
18
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Need to have some prior model of data or noise signal attributes
Ground roll suppression with polarization analysis
19
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Ground roll suppression with polarization analysis
20
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Diallo et al., 2006
Frequency-dependent polarization template matching
21
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Agenda
22
Singular Value Decomposition (SVD)
23
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Singular Value Decomposition (SVD)
24
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Scaled singular vectors
Singular Value Decomposition (SVD)
25
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Polarization vectors
Singular Value Decomposition (SVD)
26
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Data1 | | Comp 1 | Comp 2 |
| Sample 1 | 1 | 0 |
| Sample 2 | -0.2 | 0 |
Data2 | | | |
| Sample 1 | 0 | 0.2 |
| Sample 2 | 0 | 0.6 |
Data3 | | | |
| Sample 1 | 1 | 0.2 |
| Sample 2 | -0.2 | 0.6 |
Polarization vectors
Singular Value Decomposition (SVD)
27
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Data1 | | Comp 1 | Comp 2 |
| Sample 1 | 1 | 0 |
| Sample 2 | -0.2 | 0 |
Data2 | | | |
| Sample 1 | 0 | 0.2 |
| Sample 2 | 0 | 0.6 |
Data3 | | | |
| Sample 1 | 1 | 0.2 |
| Sample 2 | -0.2 | 0.6 |
Polarization vectors
Singular Value Decomposition (SVD)
28
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Data1 | | Comp 1 | Comp 2 |
| Sample 1 | 1 | 0 |
| Sample 2 | -0.2 | 0 |
Data2 | | | |
| Sample 1 | 0 | 0.2 |
| Sample 2 | 0 | 0.6 |
Data3 | | | |
| Sample 1 | 1 | 0.2 |
| Sample 2 | -0.2 | 0.6 |
Polarization vectors
Continuous Wavelet Transform
29
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
data
mother wavelet (Morlet)
daughter wavelet
dilation
time delay
Continuous Wavelet Transform
30
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
data
mother wavelet (Morlet)
daughter wavelet
dilation
time delay
Continuous Wavelet Transform
31
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
64
16
32
8
2
4
1
Agenda
32
Chevron’s six-component Kettleman survey
33
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Chevron’s six-component Kettleman survey
34
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
35
Vertical
Radial
Transverse
Yaw
Roll
Pitch
Accelerated weight-drop source
Shot interval 6.25 m
Data courtesy of Chevron
36
Vertical
Radial
Transverse
Yaw
Roll
Pitch
Vibroseis source
Shot interval 25 m
Data courtesy of Chevron
37
Vertical
Radial
Transverse
Yaw
Roll
Pitch
Dynamite source at depth = 25 m
Shot interval 25 m
Data courtesy of Chevron
38
Vertical
Radial
Transverse
Yaw
Roll
Pitch
Dynamite source at depth = 50 m
Shot interval 25 m
Data courtesy of Chevron
Rotations from geophone differencing on a free surface
39
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Pitch
r-sensor
Pitch
geo-diff
Dynamite source
d
1
2
Data courtesy of Chevron
Rotations from geophone differencing on a free surface
40
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Pitch
r-sensor
Pitch
geo-diff
Dynamite source
3 components:
Vertical
Radial
Pitch
Data courtesy of Chevron
Vibroseis source after NMO
41
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Data courtesy of Chevron
Dynamite source at depth=25m after NMO
42
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Data courtesy of Chevron
Dynamite source at depth=50m after NMO
43
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Data courtesy of Chevron
Objectives
44
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Data courtesy of Chevron
Attributes that distinguish between wave modes
45
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Data courtesy of Chevron
Agenda
46
Dynamite source at 50 m depth before removal of shear-induced wave
47
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Dynamite source at 50 m depth after removal of shear-induced wave
48
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
3C data in the continuous wavelet domain
49
Vertical
Radial
Pitch
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
CWT
3C data in the continuous wavelet domain
50
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
component
frequency
time
CWT polarization
51
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
P-wave CWT
S-wave CWT
Singular value decomposition
Scaled singular vectors =
polarizations
1st
3rd
2nd
CWT polarization
52
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
S-wave CWT
Scaled singular vectors =
polarizations
1st
Vz
Vx
Ry
CWT polarization
53
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
S-wave polarization
P-wave polarization
S-wave CWT
P-wave CWT
Polarization template
54
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
S-wave polarization template
Vz
Vx
Ry
1
0
weight
Polarization template
55
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Some other arrival’s polarization
Vz
Vx
Ry
1
0
weight
Weighting function
56
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Continuous wavelet domain polarization filtering flow
57
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
CWT
SVD
CWT+SVD
CWT-1
polarization template
Weighting function
Dynamite source at 50 m depth before removal of shear-induced wave
58
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Dynamite source at 50 m depth after removal of shear-induced wave
59
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Polarization filter vs frequency filter
60
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Input data
Polarization filter
High-pass filter
Dynamite source at 25 m depth before removal of shear-induced wave
61
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Dynamite source at 25 m depth after removal of shear-induced wave
62
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Vibroseis source before removal of Rayleigh waves
63
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Vibroseis source after removal of Rayleigh waves
64
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Vibroseis source before removal of Rayleigh waves
65
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Vibroseis source after removal of Rayleigh waves
66
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Summary
67
Summary
68
Chapter 4: Automatic wave-mode identification with machine learning
Machine-learning classification at station 337
Manual wave-mode identification / Support Vector Machines / Automatic wave-mode identification
69
Training: Class 1=
Slow ground roll
Training: Class 1=
Fast ground roll
Training: Class 1=
Slow+fast ground roll
Class 1
Class 0
Class 1
Class 0
Class 1
Class 0
Data courtesy of Chevron
Acknowledgements
Spare slides
71
Comparison of weighting function with different components
72
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
C1
C2
C3
Cone
Comparison of weighting function with different components
73
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Six components in Kettleman dataset:
Comparison of weighting function with different components
74
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Compare filtering using:
Vertical
Radial
Transverse
+ Pitch
Vertical
Radial
Transverse
vs
Comparison of weighting function with different components
75
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
C1
C2
C3
Comparison of weighting function with different components
76
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
C3
C1
C2
Comparison of weighting function with different components
77
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
C1
C2
C3
Slice
Comparison of weighting function with different components
78
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
C1
C2
C3
Wedge
Comparison of weighting function with different components
79
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
C1
C2
C3
Cone
Comparison of weighting function with different components
80
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
C1
C2
Cone
C4
Comparison of weighting function with different components
81
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Compare filtering using:
Vertical
Radial
Pitch
Vertical
Radial
Transverse
vs
Transverse vs Pitch component, dynamite 25 m
82
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Transverse
Transverse vs Pitch component, dynamite 25 m
83
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Before removal of shear-induced wave, dynamite 25 m
84
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Transverse
Pitch
After removal of shear-induced wave with vertical, radial, transverse
85
Transverse
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
After removal of shear-induced wave with vertical, radial, pitch
86
Pitch
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Transverse
Before removal of slow ground roll, vibroseis
87
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Transverse
Pitch
After removal of slow ground roll with vertical, radial, transverse
88
Transverse
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
After removal of slow ground roll with vertical, radial, pitch
89
Pitch
Transverse
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Before removal of shear-induced wave, dynamite 50 m
90
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Transverse
Pitch
After removal of shear-induced wave with vertical, radial, transverse
91
Transverse
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
After removal of shear-induced wave with vertical, radial, pitch
92
Transverse
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Vertical
Radial
Pitch
Ground roll suppression with Stack-Array
93
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
(Morse and Hildebrandt, 1989)
Ground roll suppression with Stack-Array
94
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
(Morse and Hildebrandt, 1989)
Ground roll suppression with Stack-Array
95
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
(Morse and Hildebrandt, 1989)
Vz noise in ocean bottom node data
96
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Upgoing and downgoing wavefields
97
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
(M. Wong, 2014)
Upgoing
Downgoing
Upgoing and downgoing wavefields
98
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
(M. Wong, 2014)
Upgoing
Downgoing
PZ summation:
Combine Hydrophone and vertical geophone to separate upgoing data from downgoing data
Vz noise in ocean bottom node data
99
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Hydrophone
Vertical with Vz noise
(Zhou et al., 2011)
Vz noise in ocean bottom node data
100
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Hydrophone
Vertical with Vz noise
(Zhou et al., 2011)
P wave
S wave
Vz noise suppression with match filtering
101
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Hydrophone
Vertical with Vz noise
(Zhou et al., 2011)
Vz noise suppression with match filtering
102
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Hydrophone
Vertical with Vz noise
Vertical after Vz noise suppression
(Zhou et al., 2011)
Vz noise suppression with match filtering
103
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
Hydrophone
(Zhou et al., 2011)
Ground roll suppression with adaptive subtraction
104
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
(Edme et al., 2014)
Ground roll suppression with adaptive subtraction
105
Introduction / Polarization, CWT, SVD / Kettleman 6C survey / Polarization-based wave-mode separation
(Edme et al., 2014)
106
“FD Pressure” Idea
-Put the FD Stencil in the Earth-
(Courtesy of Fred Herkenhoff)
Thanks for listening
1C
3C
6C
email: ohad@sep.stanford.edu