Reflection

Version (04.02.
2005)
Reflexionsseismik 1 Seite 3
Literature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1 What is reflection seismics? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2 Historical developments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3 Use of reflection seismics. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2. Seismic waves. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.1 Principles of wave propagation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Interface: reflection, refraction, conversion. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.3 Reflection- and transmission coefficients . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.4 Geometry of wave propagation: rays . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3. Seismic Velocities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
4. Amplitudes and Attenuation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
5. Seismic data acquisition systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.2 Seismic sources . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.3 Seismic receivers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.4 Registration unit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
6. Acquisition setup . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
6.1 Single-channel measurements (Seismic profiling) . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
6.2 Multi-channel measurements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
6.3 CDP, CMP und Zero-Offset, Common offset . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
7. Seismogram . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
7.1 The seismic trace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
7.2 Events in a seismic trace. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
7.3 Resolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
8. Basic scheme for processing of reflection data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
8.1 Aim of data processing: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
8.2 Basic framework fordata processing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
9. Loading the data, Demultiplexing, Data-Format . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
9.1 Principle of demultiplexing (Sorting of the data) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
9.2 Dataformats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
10. Amplitude correction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
10.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
10.2 Trace equalisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
10.3 Automatic Gain Control - AGC . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
10.4 Spherical divergence correction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
10.5 Programmable gain curves (time and Distance-dependent) . . . . . . . . . . . . . . . . . . . 47
11. Filter (Frequencyfilter) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
11.1 Fourier transformation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
11.2 Spectrumanalysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
11.3 High cut, low cut, Band pass filter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
12. Static corrections - Part 1 (Refraction) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
12.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
12.2 Methods for static correction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
Version (04.02.2005)
Seite 4 Reflexionsseismik 1
13. Deconvolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
13.1 Convolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
13.2 Cross-correlation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
13.3 Deconvolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
14. Velocity analysis and NMO-Correction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
14.1 Normal-Moveout-Correction (NMO) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
14.2 Methods for velocity analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
14.3 Problem of stretching the data caused by NMO-correction . . . . . . . . . . . . . . . . . . . 61
14.4 Factors influencing velocity estimates . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
15. Stacking . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
15.1 Muting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
15.2 Methods of stacking . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
16. Residual statics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
16.1 Principle of residual statics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
17. Special Filters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
17.1 Fk-Filter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
17.2 τ-p-Filter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
18. Migration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
18.1 Geometrical distortion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
18.2 Methods for migration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
18.3 Special migration and extensions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
19. Postprocessing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
19.1 Time-dependent Frequencyfilter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
19.2 Deconvolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
19.3 Coherencyfilter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
19.4 Amplitude correction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
19.5 Archiving of Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
20. Interpretation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
20.1 Mapping of geological structures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
20.2 Seismic sequence-analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
20.3 Seismic facies analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
20.4 Interpretation of 3D-data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
20.5 Other Aspects of interpretation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
21. Other related Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
21.1 VSP - vertical seismic profiling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
21.2 Cross-hole seismic . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
21.3 Sidescan sonar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
21.4 Georadar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
21.5 New developments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
22. Selection of used references (Figures etc.) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
Version (04.02.2005)
Reflection seismics 1 Page 5
Literature
Seismics:
Hatton, L., Worthington, M.H. and Makin, J. (1986). Seismic data processing - Theory and
practice. Blackwell Scientific publications, Oxford, UK, 177 pp.
Sheriff, E.G. and Geldart, L.P. (1995). Exploration Seismolgy, (2nd ed.). Cambridge
University Press, Cambridge, 592 pp.
Yilmaz, Ö. (1987). Seismic data processing. SEG Tulsa, OK, 826 pp.
Yilmaz, Ö. (2001). Seismic data Analysis: Processing, Inversion, and Interpretation of

Seismic data. SEG Tulsa, OK, 2027 pp.
Interpretation of seismic data:

Emery, D. and Myers (eds.) (1996). Sequence Stratigraphy. Blackwell Science, Oxford, UK,
297 pp.
McQuilin, R. et al. (1986). An introduction to seismic interpretation. Graham and Trotman

Ltd., London, UK.
Applied Geophysics in general (with one part about Seismics):

Dobrin, M.B. and Savit, C.H. (1988). Introduction to geophysical prospecting (4.ed).
McGraw-Hill Book Comany, Ney York, 867 pp.
Kearey, P. and Brooks, M. (2002). An introduction to geophysical prospecting. Blackwell

Scientific Publications, Oxford, 262 pp.
Reynolds, J.M. (1998). An introduction to applied and environmental geophysics. John

Wiley and sons, Chichester, UK, 796 pp.
Telford, W.M., Geldart, L.P. and Sheriff, R.E (1990). Applied Geophysics (2. ed.).
Cambridge University Press, Cambridge, 770 pp.
Version (04.02.2005)
1. Introduction
The aim of this lecture is to discuss the basic principles of Reflection Seismics and to explain
the basic fundamentals.
Knowledge of the lectures “Allgemeine Geophysik I und II” is indispensable.

Knowledge of the lecture “Digitale Verarbeitung seismischer Signale” is an advantage.
The topic reflection seismics is that large that only for some parts an extensive mathematical
discussion is carried out. Therefore, this lecture only discuss qualitatively several topics which
are important for the reflection seismics. An extensive treatment of specific processingsteps and
the mathematical background will be treated by “Praktikum zur Datenverarbeitung in der
Reflexionsseismik”.
1.1 What is Reflection seismics?
In seismics, the geology is examined using seismic waves. The aim is to recognize geological
structures and, if possible, to determine the material properties of the subsurface.
Source Receiver Raw data

Q t t t t t t
geological subsoil seismic Section
Version (04.02.2005)
Page 8 Reflection seismic 1
1.2 Historical Developments
Table 1: Historical development of reflection seismics(Sheriff and Geldart, 1995)

1914 Mintrop’s mechanical seismograph 1954 Continuous velocity logging
1917 Fessenden patent on seismic method 1955 Moveable magnetic heads
1921 Seismic reflection work by Geological 1956 Central data processing
Engineering Co.
1923 Refraction exploration by Seismos in Mexico 1961 Analog deconvolution and velocity filtering
and Texas
1925 Fan-shooting method 1963 Digital data recording
Electrical refraction seismograph
Radio used for communications and/or time-
break
1926 Reflection correlation method 1965 Air-gun seismic source
1927 First well velocity survey 1967 Depth controllers on marine streamer
1929 Reflection dip shooting 1968 Binary gain
1931 Reversed refraction profiling 1969 Velocity analysis
Use of uphole phone
Truck-mounted drill
1932 Automatic gain control 1971 Instaneous floating-point amplifier
Interchangeable filters
1933 Use of multiple geophones per group 1972 Surface-consistent statics
Bright spot as hydrocarbon indicator
1936 Rieber sonograph 1974 Digitization in the field
first reproducible recording
1939 Use of closed loops to check misties 1975 Seismic stratigraphy
1942 Record sections 1976 Three-dimensional surveying
Mixing Image-ray migration (depth migration)
1944 Large-scale marine surveying 1984 Amplitude variation with offset
Use of large patterns Determining porosity from amplitude
DMO (dip-moveout) processing
1947 Marine shooting with Shoran 1985 Interpretation workstations
1950 Common-midpoint method 1986 Toiving multiple streamers
1951 Medium-range radio navigation 1988 S-wave exploration
Autopicking ot 3-D volumes
1952 Analog magnetic recording 1989 Dip and azimuth displays
1953 Vibroseis recording 1990 Acoustic positioning of streamers
Weight-dropping GPS satellite positioning
Version (04.02.2005)
1.3 Use of Reflection seismic

Reflection seismic tools are used in many different ways. Depending on the aim and the use
different depth of penetration and resolution are obtained.
Applications that do not belong to the area of the classical seismic (e.g. Georadar, echolot) use
reflection seismic principles (traveltime curves, filters etc.).
Table 2: Application of reflection seismic tools
Aim Penetration Depth

Oil/gas Exploration of layers 100 m - 5 km
Groundwater
Engineering geo-
Polution 10-500 m
physics
Archaeological investigation
Earth Crust seismic Composition of Earth/Geodynamic 1-60 km
Measurements in Echo sounder, ~0m

water high resolution Seismics <100 m
Georadar GPR Shallow investigation of the Earth 0.5-10 m
Version (04.02.2005)
2. Seismic waves
The basic of Seismics are the elastic waves, which include both P- and S-waves
Acoustic waves include only P-waves (in water)
(1) Body waves:

P-Waves (also Longitudinal- or compression waves)
S-Waves (also Transversal- or Shearwaves)
(2) Surface waves exist on the interface between two different media:
Rayleigh-Waves (Surface seismic wave propagated along a free surface of a semiinfinite
medium)
Love-Waves (Surface seismic channel wave SH wave
Propagation of different waves
Version (04.02.2005)
2.1 Principles of wave propagation

The propagation of seismic waves is described by the wave equation. This wave equation can
be derived from the relation between tension, elasticity, Hook’s law and Newtons law. In three
dimension the derivation of this wave equation is quite complicated. This is why a simplified
medium is assumed which exist only of one dimension. Newtons (second) law (in one
dimension) is given by
2
∂ ∂
P = –ρ 2 Uz
∂z ∂t
where P is the acoustic pressure and Uz is the displacement and ρ is the mass density. The well-
known formula F=m*a can be recognised in this equation and describes that a force acting on a
certain mass results in an acceleration of that mass.
Hooke’s law is given by
∂
U = –κ P
∂z z
where κ is the compressibility an relates the stress (force per unit area) and strain (change of
dimensions or shape). The bulk modulus k is the reciprocal of the compressibility and is given
by
1
k = ---
κ
Combining these two equations we obtain the Acoustic wave equation:
2 2
∂P 1 ∂
2
– ----2- 2 P = – w ( t )δ ( z )
∂z v ∂t
w(t) is the source signal and v is the wavespeed which is given by
1 k
v = ----------- = ---
ρκ ρ
We will not use these expressions for the wave equation, because for most topics treated in this
lecture, it is sufficient to consider the wave equations in a geometrical way, as wavefronts or as
rays. The wave equation is mainly used for modelling and inversion of seismic waves. It will be
shown that the wave velocities depend on the compressibility and the density.
Version (04.02.2005)
2.2 Interface: Reflection, Refraction, Conversion.
When a wave encounters an interface three phenomena can occur:
(1) Reflection
(2) Refraction
(3) Conversion
s
β1
p
α1 p
α1
vp1,vs1
vp2,vs2
β2
α2
s p
Interface: α=Angle of P-Waves, β=Angle of S-Waves.
Reflection
• Angle of incidence = angle of reflection (α1 = α2)
Refraction (change of direction of a seismic ray upon passing into a medium with a different
velocity.)
sin α 1 v1
-------------- = -----
sin α 2 v2
• Critical angle: α2 = 90°, the refracted ray grazes the surface of contact between two media
sin α 1 v1
---------------- = sin α 1 = -----
sin 90° v2
Conversion
• Change at the interface of p->s or s->p
• As well as refraction as reflection is possible.
Law of Snellius sin α

For all events at an interface the ratio, ----------
-
v is always the same.
This ratio is also called the Raypath parameter p. The general Form of Snell’s law is given by:
sin α 1 sin β 1 sin α 2 sin β2
-------------- = ------------- = -------------- = -------------- = p = const
v p1 v s1 v p2 v s2
Version (04.02.2005)
2.3 Reflection- and transmission-coefficients
To derive the reflection and transmission coefficients for elastic waves, the boundary conditions
at the interface are needed. These reflections coefficients depend on
• Difference in density
• Difference in velocity
• Angle of incident of the wave
and are described by the Zoeppritz-Equations.
The Reflection- and Transmissioncoefficient give the ratio between the incident amplitude A0 ,
and the reflected (AR) and transmitted (AT) amplitude, respectively. In the special case of a
incident wave perpendicular at an interface for a P-wave, a simple expressions for the reflection
and transmission coefficient is obtained.
AR v2 ρ2 – v1 ρ1 Z2 – Z1
Reflectioncoefficient: R = ------- = ----------------------------- = ------------------
A0 v2 ρ2 + v1 ρ1 Z2 + Z1
AT 2v 1 ρ 1 2Z 1
Transmissioncoefficient: T = ------- = ----------------------------- = ------------------
A0 v2 ρ2 + v1 ρ1 Z2 + Z1
The Product Z = vρ is introduced as the acoustic Impedance.
Sometimes the coefficients which describe the Energy and not the amplitudes are introduced as
Reflection- und Transmissioncoefficients:
2
( Z2 – Z1 )
Reflection coefficient: E R = -------------------------2-
( Z2 + Z1 )
4Z 1 Z 2
Transmission coefficient: E T = -------------------------2-
( Z2 + Z1 )
Obviously the total amount of energy is the same before and after the reflection and transmis-
sion, so that : E R + E T = 1
In a general case these coefficients are depending on the angle of incidence. Moreover, also
conversions between P- and S-waves occur at an interface. This results in complicated expres-
sions which will not be discussed in this report.
Version (04.02.2005)
Angle-dependent reflection- and transmission-coefficients for P- and S-waves
For the special case of v2/v1=2 and ρ2/ρ1=0.5, we have Z1=Z2. The P-wave reflection coeffi-
cient is essentially zero for small angles of incidence. As the critical angle for P-waves is
approached, the transmitted P-wave energy falls rapidly to zero and no transmitted P-wave
exists for larger incident angles.
2.4 Geometry of wave propagation: rays

The travel time can be derived by analysing the geometry of the rays that travel from source to
receiver (when the velocity and depth are known).
=> geometrical Seismic
Version (04.02.2005)
Traveltime curve
The traveltime of the different rays are plotted in an x-t diagram, a Traveltime plot. The picture
below shows such a traveltime plot.
t0
ti
Picture of a typical travel time diagram for a two-layer case.

t0 = t0-time of the Reflection; ti=Intercept-time of the Refraction; xcrit =critical distance;
xcross=”cross-over”-distance.
The important rays in the travel time diagram are:

• Direct wave
• Reflected wave
• Refracted wave
The separate elements will be discussed in the following paragraphs:
Direct wave
The direct wave travels only in the upper layer, directly from the source to the receiver. When
the upper layer has a homogeneous velocity distribution than we can assume that the wave
travels parallel to the interface. In the traveltime diagram the direct wave shows a slope of 1/v.
Reflected wave in the horizontal two-layer case

In Reflection seismic this signal is investigated and used. Therefore, we discuss this wave in
more detail.
Version (04.02.2005)
The simplest case is the horizontal two-layer configuration.
A B
o x/2 x/2 x
s
h
s
α
α
v
C
Schematic overview to calculate the traveltime expressions for a reflection.
Using pythagoras we have:
x 2
s = h + ⎛ ---⎞
2 2
⎝ 2⎠
which gives:
2 2 2
24s 4h + x
t = -------- = ---------------------
2 2
v v
The traveltime of a reflected wave goes for large x asymptotically towards the traveltime of the
direct wave (t=x/v). Rewriting last equations renders:
2 2 2
v t x
---------2 – --------2 = 1
4h 4h
This is the expression for a hyperbola.
t0
The crossing with the time-Axis (x=0) is the timezero, t0:

2 2 2 2h
t 0 v = 4h ⇒ t ( x = 0 ) = t 0 = ------
v
Version (04.02.2005)
With the t0-time the expression for the traveltime can be written as:
2
2 x 2
t = ----2- + t 0
v
Moveout
Rewriting the former equation, and using a binomial expansion of the square root
2⎛ x ⎞
2
2
t = t 0 ⎜ 1 + -----------2-⎟
⎝ t v ⎠
2
0
.
⎛ x ⎞
2
t = t0 ⎜ 1 + -----------2-⎟
⎝ t v ⎠
2
0
⎛ x ⎞
2
x
2
t ≈ t 0 ⎜ 1 + ---------------2⎟ = t 0 + ------------2-
⎝ 2t 0 v ⎠
2
2t 0 v
we obtain an approximate expression for the moveout in time between two different positions
x1 and x2.
2 2
x2 – x1
t 2 – t 1 = ---------------------
2
2t 0 v
Normal Moveout
The difference between the traveltime t for a specific receiver at a specific distance x and a
receiver at zero-offset x0=0 is given by the normal moveout which is defined as:
2
x
∆T = t x – t 0 ≈ ------------
2
-
2v t 0
This formula will later be used in the processing of seismic data.
Version (04.02.2005)
Reflections - horizontal layered medium

When more than two layers are present the rays of the reflections are more complex due to the
refraction that occurs at each interface. With the help of Snell’s law an approximate estimate of
the rays can be obtained iteratively..
v1
v2
v3
Three horizontal layers as example for a layered horizontal medium. The solid line represents
the real travelpath, where the dashed line represents the travel path of one layer with equal
travel time.
It is possible to replace the expressions with the expressions for the one layer case when for the
layers above the actual reflection point an average velocity is assumed.
A possible average velocity is given by.
n
∑ vi τi
i=1
v average = -----------------
n
-
∑ τi
i=1
where zi is the thickness of the τi is the interval travel time through the ith layer and vi
ith layer,
th
is the velocity of the i layer.
For small offsets the travel-time curve is still essentially hyperbolic so that using the hyperbolic
expression an average velocity can be obtained for the layers above the actual reflection point.
Therefore, a better approximation is obtained when instead of the arithmetic average, the arith-
metic average is taken of the square of the velocities(“rms: root-mean-square”).
n
∑ vi
2
τi
2 i=1
v rms = --------------------
n
-
∑ τi
i=1
The so-called RMS-velocity is used for all velocity analysis of reflection data, because it is
relatively easy to obtain it from a seismogram.
The velocity that is used to remove the effect of normal moveout is called stacking velocity
(This NMO-correction will be discussed later in more detail). Correcting for the moveout using
the stacking velocity and stacking the traces results in the maximum amplitude of the reflection
event. As mentioned before, the travel-time curve for reflected rays in a multi-layered ground
is not a hyperbola. However, if the maximum offset is small compared with reflector depth, the
stacking velocity closely approximates the root-mean-square velocity.
Version (04.02.2005)
To calculate the interval-velocities from the RMS-velocities, the Formula of Dix is used:
1
2 2 ---
v rmsn t n – v rmsn – 1 t n – 1 2
v i = --------------------------------------------------
tn – tn – 1
vi is the intervalvelocity of the nth layer. tn is the measured traveltime to the nth layer. When
t n ≠ t 0n (travel time for zero offset) we assume that the travel path to the individual reflections
are equal. This is in general not correct. Dix’ Formula can then result in wrong vi - values.
Reflection from a dipping reflector

The reflectors in reality are often not horizontal but dipping.
∆Tdip
-x x
90 + Θ x
h
Θ
h
Ray path and traveltime diagram for a dipping Reflektor.
With the help of the cosinus rule an expression for the traveltimes can be derived as follows:
2 2 2 2
v t = x + 4h – 4hx ⋅ cos ( 90 + Θ )
2 2 2 2
v t = x + 4h + 4hx ⋅ sin Θ
The derivation of the next formulas for the dipping reflector follows analogue the case of the
horizontal layer, with the exception that the expressions include a term with the angle Θ. We
obtain a hyperbolic expression:
2 2 2
v t ( x + 2h sin Θ )
--------------------------2- – -----------------------------------
2
= 1
( 2h cos Θ ) ( 2h cos Θ )
Where the symmetry axis does not lie on the t-axis, but is now the line x = – 2 h sin Θ .
2 2 2 2
v t = ( 2h cos Θ ) + ( x + 2h sin Θ )
2 2
2 ( 2h ) ( x + 2h sin Θ )
t = ------------
2
- 1 + -----------------------------------
2
v 4h
Version (04.02.2005)
(Where cosΘ is assumed to be 1 for small angles). After expansion as a binomial expression we
obtain: x 2 + 4hx sin Θ 1 / 2
t = t 0 1 + ---------------------------------
2
-
4h
x 2 + 4hx sin Θ
t ≈ t 0 1 + ---------------------------------
2
-
8h
2h
where t0 is given by t 0 = ------
v
The traveltime from two locations x und -x, that both have the same distance from the shot
position is not equal anymore caused by the dipping reflector. The traveltime difference
between both points is called the “dip-moveout”.
2x sin Θ
∆T dip = t x – t – x = -------------------
v
Refraction - horizontal Two-layer case
For a horizontal two-layer case the following situation is given:
1
t -----
v1
1
------
ti v
2
xc x
x Horizontal two-layer case with the rays and the

traveltime diagram above. The slope of the
h Refracted wave is 1/v2.
v1
v2
For the horizontal two-layer case the traveltime is given by:
x 2h cos Θ c
t = ----- + ----------------------
v2 v1
or with sinΘc=v1/v2, 2 2
x 2h v 2 – v 1
t = ----- + --------------------------------
v2 v1 v2
x
t = ----- + t i
where v2
and ti die Intercept-Time is. That is the time, that the Refracted wave crosses the time axis.
The “crossover distance” xc, the distance at which the travel time of the direct wave equals that
of the refracted wave is derived as follows:
xc x c 2h cos Θ c
----- = ----- + ----------------------
v1 v2 v1
and solving this for xc gives: v2 + v1
x c = 2h -----------------
v2 – v1
Note that again sinΘc=v1/v2 is used.
Version (04.02.2005)
3. Seismic Velocities
The velocities is the most important parameter in reflection seismics. Information about the
velocities is important for:
• Change of traveltime with depth

• Check of the seismic data with a model
• Correct for the geometry (Migration)
• Classification and filtering of signal and noise
• Geological and lithological Interpretation
The seismic velocity depends on the elastic parameters as follows:

4µ
k + ------
3
vp = ----------------
ρ
µ
vs = ---
ρ
with k= Bulk modulus
µ = Shear modulus
From these equations we observe that:

• vp is always larger than vs
• In Liquids, µ = 0 => no shear waves can propagate.
Poisson-Ratio
The Poisson-Ratio σ is the ratio of transverse strain to longitudinal strain and is defined as:
3k – 2µ
σ = ------------------------
2 ( 3k + µ )
and has values between 0 < σ < 0.5 (0.5 is valid for fluids).
Seismic velocities and the factors influencing it.

In Reality the factors k, µ and ρ depend in a complex way on different factors:
• Lithology (Matrix und Struktur)

• Depth
• Interstitial fluids
• Pressure
• Porosity
• Degree of compaction
• Cementation
• etc.
Most of these factors influencing the factors k, µ and ρ are empirically determined.
Version (15.02.2005)
Page 22
Version (15.02.2005)
Summary of effects of different rock properties on P- and S-wave velocities and their
ratios.(Sherriff and Geldart 1995).
Reflection seismic 1
Examples for typical velocities
material vp (km s-1)

Unconsolidated Material
Sand (dry) 0.2 - 1.0
Sand (water saturated) 1.5 - 2.0
Clay 1.0 - 2.5
Glacial till (water saturated) 1.5 - 2.5
Permafrost 3.5 - 4.0
Sedimentary rocks
Sandstone 2.0 - 6.0
Tertiary sandstone 2.0 - 2.5
Pennant sadstone (Carboiferous) 4.0 - 4.5
Cambrian quartzite 5.5 - 6.0
Limestones 2.0 - 6.0
Cretaceous chalk 2.0 - 2.5
Jurassic oolites and bioclastic limstones 3.0 - 4.0
Carbiniferous limestone 5.0 - 5.5
Dolomites 2.5-6.5
Salt 4.5 - 5.0
Anhydrite 4.5 - 6.5
Gypsum 2.0 - 3.5
Igneous / Metamorphic rocks

Granite 5.5 - 6.0
Gabbro 6.5 - 7.0
Ultramfic rocks 7.5 - 8.5
Serpentinite 5.5 - 6,5
Pore fluids
Air 0.3
Water 1.4 - 1.5
Ice 3.4
Petroleum 1.3 - 1.4
Other materials
Steel 6.1
Iron 5.8
Aluminium 6.6
Concrete 3.6
(Kearey and Brooks, 1991)
Version (15.02.2005)
Anisotropy
Seismic velocities in a solid can depend on the direction in which they travel within the solid.
(for example Granite). This is the case when the matrix in the solid is direction dependent
caused by for example gravitational forces .
Anisotropy: A seismic wave travels faster in that direction where the matrix is more dense.
Direct measurements of the seismic velocity
Velocities can be determined when measuring the time, that a seismic wave needs to travel a
certain distance. Direct measurements are carried out in a laboratory in a small probe or in situ
in a borehole
Problem:
Many times higher frequencies are used to determine the seismic velocities. It is important that
the velocity is frequency dependent, especially in a heterogeneous solid.
Version (15.02.2005)
4. Amplitudes and Attenuation
The amplitudes of a seismic signal are influenced by several factors:
Phenomena causing the degradation of a seismic wave (Reynolds, 1998).
• Divergence
• Absorption
• Scattering of Energie at an interface
(Reflection, Refraction, Conversion)
• Interference with other waves(e.g. multiple Reflections)
• Spreading of Energie
• Influence of measurement system
Spherical divergence (geometrical spreading)
When a seismic wave propagates, the energy is spread along

the surface of the wavefront. For cylindrical and spherical
waves the surface increases with increasing radius. That is
why the energy decreases. The Energy of the different waves
are proportional to: r1
• plane Wave: konstant r2

• cylindrical Wave: ~ 1/r
• spherical Wave: ~ 1/r2
The energy of the body (spherical) and interface (cylindrical)

waves are proportional to (1/r2) and (1/r), respectively. The amplitude of the body and interface
waves are proportional to (1/r) and (1/or), respectively (energy is proportional to (amplitude)2).
Version (15.02.2005)
Absorption
When elastic waves travel through the subsurface, a part of their energy is dissipated in heat.
The decrease of amplitude due to absorption appears to be exponential with distance for elastic
waves in rock which can be written as:
– αx
A = A0 e
A and A0 are the amplitudes at two locations with a distance x and α is the absorptioncoeffi-
cient.
Quality factor Q
The absorption is often described by the Quality factor Q (Q-Factor) that is given by:
2π 2π
Q = -------------------- = -------------------------------------------------------------------------------------
( ∆E ) ⁄ E Fraction of Energie, lost per Cycle
The absorption coefficient α can be considered as a first approximation proportional to the
frequency. We can write:
1 αv αλ
---- = ------- = -------
Q πf π
This is only valid for Q >> 1 (where v = Velocity and f = Frequency).
Version (15.02.2005)
6. Acquisition setup
6.1 Single channel measurements (profiling)

Only one source and receiver are used with often an equal distance between the source and
receiver. This is repeated for several positions along a line.
Typical Applications:
• simple marine System (Boomer, Sparker etc. --> Profiling)
• Georadar
6.2 Multi channel measurements

Multi channel systems use one source and several receivers, which measure at the same time.
Several spreads are possible to orient the sources and receiver:
single ended (end on) single ended (inline offset)
+l l l l l l l l l + l l l l l l l l l
split spread fan shooting
+l
l l l l l l l l l
l l l l l l l l l l l
l l l l l l + l l l l l l +
Types of reflection spreads. The symbol o and + represent source and geophone-group center
locations, respectively.
“roll along”
In land seismic, often more geophones are put in the field than the measurement system can use.
In this way the active channels can be connected when needed and there is no need to change
the total field setup.
Principle of “Roll along” - Measurement
Version (15.02.2005)
6.3 CDP, CMP and Zero-Offset, Common offset

There are different possibilities to sort the data:
• Common shot - all traces, that belong to the same shot

• Common midpoint (CMP) - all traces with the same midpoint
• Common receiver - all traces, recorded with the same geophone
• Common offset - all traces with the same offset between shot and geophone
Common mid point - CMP
Difference between CMP und CDP: For a horizontal Reflector all traces that have the same
midpoint, have also the same reflectionpoint in the subsurface. Is the layer inclined than the
traces have a different reflectionpoint.
For a horizontal layerered earth the reflectionpoint lies between source and receiver (midpoint)
Using more shots with different positions of the source and receivers several combinations of
source and receivers exists which have the same midpoint. When a horizontal layering is present
the reflection then also come from an equal point in the subsurface (Common depth point-
CDP). For an inclined layer the point of reflection for traces with the same midpoint are not
equal anymore. The nomenclature CDP is not valid anymore. However several processing
programs still use the word CDP in stead of CMP.
Zero offset
Zero offset data is characterised when the source and receiver are present on the same location.
There is no moveout. For a normal measurement this is seldom the case. When the traces are
corrected for the moveout and are stacked then a zero offset trace is obtained.
Version (15.02.2005)
Common offset
All traces with equal offset between source and receiver. This configuration is often used for
several Single channel systems. Also Georadar measurements are often carried out with a fixed
offset between source and receiver.
Fold
The fold indicates the number of traces per CDP. This is often the number of traces in a CMP.
The theoretical formula for the fold is given by:
Numberof Geophones ⋅ Distance between Geophones
Fold = --------------------------------------------------------------------------------------------------------------------------------------
2 ⋅ Distance between shots
The number of traces which are measured at a certain geophone position is called “surface fold”
Crooked line method
Sometimes access and/or structural complications make it impossible to locate lines in desired
locations. The field recording may be done in the same way as CMP surveying, except that the
line is allowed to bend. Usually, a best-fit straight line is drawn through the midpoint plot,
rectangular bins are constructed and those traces whose midpoints fall within a bin are stacked
together. The bins are often perpendicular tot the final line, but sometimes bins are oriented in
the strike direction. Because the actual midpoint locations are distributed over an area, they
contain information about dip perpendicular to the line and in effect produce a series of cross-
spreads, from which the true dip can be resolved. Lines are sometimes run crooked intentionally
to give cross-dip information.
3D-Seismic
The sources and receivers do not lie on a line, but in both horizontal directions sources and
receiver are placed. In this way not only for different offsets, but also information from different
directions (azimuths) is obtained. Here also bins are defined and all traces of the specific CMP
are gathered after processing in that bin.
Version (15.02.2005)
7. Seismogram
7.1. The seismic Trace

In Seismic reflection surveys we measure the ground motion over a short interval of time
following the triggering of a nearby seismic source. The graphical plot of the output of a single
detector is called a Seismogram or seismic trace;
From a geological subsoil to a seismic trace.
The properties of the geological subsoil (density and seismic velocity) determine the acoustic
impedance of a layer. From these impedances the reflectivity function of an interface can be
derived (see Chapter. 2). This function is convolved with the signal of the seismic wave. The
result is a seismic trace, on which also noise is add.
Remark: note that the wave travels half of the time downwards and the other half upwards. This
must be taken into account when the source signal is convolved with the reflectivity function.
(The principle of convolution is discussed more thoroughly when the data processing is
discussed).
Waveform
The most important waveforms in seismic are shown in the figure below and are the
• Minimum-Phase wavelet
• zero phase wavelet
Version (15.02.2005)
The most important waveforms in seismics:

(a) Minimum-Wavelet and
(b) Zero-Phase Wavelet.
Both waveforms are shown with normal and
reverse polarity.
Representation
The traveltime is in general shown with increasing time along the vertical increasing
downwards (larger traveltime corresponds to a larger depth).
There are several ways to represent seismic traces. The sort of representation depend on the
processing used, but also on the number of traces.
Different representations of a seismic trace.
“Wiggle”
For the simplest representation the amplitude is depicted as a curve (Wiggle).
“Variable area+wiggle”
When there are more traces then the result is disordered. The right half of the trace is drawn
black. Standard (set by SEG) is: the positive half of the wave on the right site is colored black.
Version (15.02.2005)
This is in seismics the most used representation.

“Variable area”
When a lot of traces are depicted close to each other, then most of the time only the positive half
of the traces is plotted. (e.g. for smaller version of seismic sections.). To suppress noise one
often plot only a part of the half of the waves(Variable Amplitude).
“Variable density”
For the interpretation the amplitudes are often plotted in different grayscales or colours
(“variable density”). This is standard for Georadar or seismic Interpretation. In this way the
differences in amplitudes are more clear.
7.2 Events in a seismic trace

The important elements are
• Reflections
• Refractions
• Interface waves
• Multiples
• Noise
Reflections, Refractions and interface waves are already discussed in Chapter 2.
Multiples
1) long-path multiples (occur when exceptionally large reflectioncoefficients are present):

-> Ghost reflections, where rays from a buried explosion on land (or an airgun in water) are
reflected back from the ground surface (or sea surface) to produce a reflection event, known as
a ghost reflection, that arrives a short time after the primary.
-> Water layer reverberations, where rays from a marine source are repeatedly reflected at the
sea bed and sea surface
2) Short-part multiples (“peg-leg multiple”):
Involve only a short additional path length to arrive so soon after the primary event that they
extend the overall length of the pulse. (Multiples between two interfaces of a layer)
Example for the travelpaths of multiple reflections.
Version (15.02.2005)
Noise
The S/N ratio “Signal-to-noise ratio” gives the ratio between the amplitude of a signal (e.g.
Reflections) to background noise(“random noise”) or noise sources (“coherent noise”).
One aim of the dataprocessing is to increase the S/N ratio.
7.3. Resolution
Vertical resolution
Comparison between the wavelength of a 30-Hz Signal, Big Ben and a Log of a drilling.
Vertical resolution means: How thick must a layer be, to discern the top and bottom of the
specific layer. Theoretically, a layer can be distinguished when it has a thickness of 1/4
wavelength (Rayleigh-Kriterium).
The wavelength is determined by the ratio of the velocity and the frequency of the seismic wave:
λ = v / f.
Lateral Resolution
The lateral resolution depends on the distance between the source and receiver at the surface
and the depth of the layer. Energy that is returned to a detector within half a wavelength of the
initial reflected arrival interferes constructively to build up the reflected signal, and the part of
the interface from which this energy is returned is known as the first Fresnel zone, or, simply,
Fresnel zone.
The width of the Fresnel zone represents an absolute limit on the horizontal resolution of a
reflection survey since reflections separated by a distance smaller than this cannot be individ-
ually distinguished. The width w of the Fresnel zone is related to the dominant wavelength λ of
the source and the reflector depth z by
Version (15.02.2005)
w = 2zλ
for z>> λ.
Principle of the Fresnel-Zone
Version (15.02.2005)
11. Filter (Frequency filter)
Frequency filter are the most important filters in digital signal processing. The energy of reflec-
tions are most of the time present in a certain frequency range. Specific noise sources and
background noise are commonly present in a different frequency range and a separation of noise
and reflection information is possible.
11.1 Fourier transformation

The basis of a digital frequency filter is the Fourier transformation. The principle of Fourier
says that each signal can be described by a sum of Sinus- and Cosinus functions.
The Fourier transformation transforms a time queue from time domain (amplitude as function
of time) to the frequency domain ( amplitude as function of frequency).
∞
∫ g ( t )e
– i2πft
Fourier transformation: G(f) = dt (from time to frequency)
–∞
∫ G ( f )e
i2πft
Inverse Fourier transformation: g ( t ) = df (from frequency to time)
–∞
The Function in Frequency domain G(f) represents the amplitude and the phasedifference of a
Sinus or Cosinus function with the frequency f. We distinguish the Amplitude spectrum A(f)
and the Phasespectrum φ(f):
A(f) Amplitude-Spectrum
g(t) G(f) = A(f) e iφ(f)
φ(f) Phase-Spectrum
Because in frequency domain the data consists of an amplitude and phase (or real and imagi-
nairy part), often the amplitude spectrum of the energy, the “Powerspectrum” is shown. The
Fourier transformation is numerically evaluated very efficiently by the FFT - the Fast Fourier
Transform.
11.2 Spectrumanalysis
In seismics, we often know only roughly the interesting frequency range and the frequency
range of any noise sources. To know what frequency range is interesting for further processing,
the data must be tested with different filters and a comparison of these results indicate the inter-
esting frequency range.
Version (15.02.2005)
-.
Examples of different signals and their Fourier transformation.
11.3 High cut, Low cut, Band pass and Notch filter
There are different type of filters:
• high cut, low pass
• low cut, high pass
• Bandpass
• Notch filter
Most filters are applied in the frequency domain. For example, using the Low pass filter, the
amplitude of all frequencies above a certain frequency are put to zero. Similar filtering is
applied for the High pass and Band pass filters
Notch-Filter
A notch-filter is used to suppress one specific frequency, for example 50-Hz-Noise due to
electrical powerlines.
Version (15.02.2005)
100%
0%
Low pass High pass
100%
0%
Bandpass Notchfilter
Examples of different type of filters.
Ringing
Ideal filters have very steep edges, which is in the ideal case discontinuous.
The result of such discontinuities is ringing (also known as Gibbs’ phenomena) as can be seen
in the figure below. In figure (a), the frequency domain result is present with a discontinuity at
fc. This discontinuity results in an oscillatory behaviour in the time domain as can be seen in
figure (b). When a taper is used to make the transition more smooth at fc(c), the oscillarory
behaviour is less pronounced:
Principle of ringing which occurs in time domain when sharp boundaries are present in the
frequency domain.
In practise, it is not possible to design a frequency filter which is discontinuous at the edges.
To define this tapering a slope is given as function of the frequency as follows:
• Two Frequencies in between the taper is present (Begin and End of the tapering)
• The slope of the taper is given as function of the Frequency in dB
Version (15.02.2005)
Time-dependent Frequency filter

To take into account the larger decrease in amplitude for the higher frequencies for increasing
travel time, a time-dependent frequency filter can be used. One defines different frequency
filters for different timewindows.
12. Static Correction - Part 1 (Refraction)
12.1 Introduction
Reflections in different traces are not always lying on a hyperbole for a horizontal reflector, but
sometimes they have a certain displacement due to different lengths of the raypaths.
Effect of topography on measured data. The reflections are not lying on a hyperbola. After
static corrections the reflection appears as if source and receiver had been positioned at the
datum level (Brouwer and Helbig, 1998).
Causes for the static displacements in the data:

• Topography, i.e. Source and receiver are present at different vertical positions.
• Different depths of the boreholes in which the explosives were fired
• Weathered layer with a relative slow velocity.
Version (15.02.2005)
Aim of static corrections

Adjust the seismic traces in such a way that the sources and receivers are present at one
horizontal level. To achieve this, the travel times of the separate traces are corrected.
Static Correction: The whole trace is corrected with the same time shift
Dynamic Correction: Different time windows in the trace are corrected differently. This
results in stretching and compression of the events (e.g. NMO-
Correction. This will be discussed in Section 14. )
12.2 Methods for static Correction

• Topographic Correction (elevation statics)
• “Uphole”-Correction
• Refraction statics
Topographic correction
Vertical aligning of the different elevations of sources and receivers.
Shot-Static = (Elevation of source - Elevation of reference level) / Velocity.

Receiver-Static = (Elevation of receiver - elevation of the reference level) / Velocity.
Correction time for a trace = Shot-Static + Receiver-Static
“Uphole-Static”
=> Correction for the weathered (low velocity) layer (Area above water layer where pores are
filled with air rather than water).
When a shot is fired, also the traveltime to the surface is recorded and from this traveltime, the
velocity of the weathered layer can be estimated.
Refraction statics
=>Correction for the weathered layer.
Using the first breaks of a certain shot (refracted energy) a model can be construced for the
weathered layer (velocities and depth).
When the distance between the receiver is too large, sometimes supplementary refraction
measurements are carried out
Methods to determine the velocity and depth of the weathered layer using refractions
• Delay-Time
• GRM (generalised reciprocal method)
• DRM (deminishing residual matrices)
All these methods return an averaged model. Very small displacements between the traces are
not completely corrected for. For these small displacements corrections can be carried out
which follow the stacking and the determination of a velocity model. These correcions are
called Reststatic .
Version (15.02.2005)
13. Deconvolution
13.1 Convolution
Convolution is a mathematical operation defining the change of shape of a waveform resulting
from its passage through a filter.
Convolution
y(t) = g(t) ∗ f(t)
Output Input signal Filterfunction
The asterix denotes the convolution operator. In seismics, we obtain a response for a certain
model by convolving the seismic signal of the source with the reflectivity function.
Convolution of the reflectivity function with the signal of the source returns the seismic trace.
Version (15.02.2005)
In reality the measured signal gt consists of the Convolution of several factors:

( equivalent wavelet wt)
⎫
⎪
⎪
⎪
⎬
⎪
⎪
⎪
⎭
gt = kδt ∗ st ∗ nt ∗ pt ∗ et + Noise
Impuls Near-surface Reflektivity

of source zone of source of the Earth
Source effekt Additional

modifying effects
Mathematically the convolution is defined as follows:
m
yk = ∑ gi ⋅ fk – i
i=0
where k = 0 ... m+n; gi = (i=0 ... m) and fj = (j= 0 ... n).
The convolution can also be performed in the Fourier domain:

Convolution in time domain = Multiplication (of the Amplitudenspectrum and Addition of the
Phase spektrum) in Fourier domain.
f0 f1 f2 f3
4 3 2 1
1 0 2 4x2=8
1 0 2 2x3=6
1 0 2 1x4+2x2=8
8 6 8 5 2 1
1 0 2 1x3+2x1=5
1 0 2 1x2=2
1 0 2 1x1=1
g2 g1 g0
Example of a convolution
13.2 Cross-correlation
The cross-correlation function is a measure of the similarity between two data sets. One dataset
is displaced varying amounts relative to the other and corresponding values of the two sets are
multiplied together and the products summed to give the value of cross-correlation.
The cross-correlation is defined as:
Version (15.02.2005)
φ xy ( τ ) = ∑ xi + τ ⋅ yi
i
where xi: (i=0 ... n); yi: (i= 0 ... n); φxy(τ): (-m < τ < +m) with m = max. displacements.
Similar to the convolution, the cross-correlation can also performed in the Fourier domain.:
Cross-correlation = Multiplication of Amplitudes and Subtraction of Phase spectrum.
Auto-correlation
The Auto-correlation is a Cross-correlation of a function with itself. It is mathematically defined
as:
φ xx ( τ ) = ∑ xi + τ ⋅ xi
i
where xi = (i=0 ... n);φxx(τ) = (-m < τ < +m) and m = max. displacement.
To make the cross-correlation and auto-correlation of different traces comparable they are
normalised as follows:
φ xx ( τ )
Auto-correlation φ xx ( τ ) norm = ----------------
φ xx ( 0 )
φ xy ( τ )
Cross-correlation φ xy ( τ ) norm = -------------------------------------
φ xx ( 0 )φ yy ( 0 )
13.3 Deconvolution
The aim of deconvolution is the reverse of convolution in such a way that the reflectivity
function is reconstructed. In practice one obtains not the real reflectivityfunction, but it results in
• a shortening of the Signals
• Suppression of Noise
• Suppression of Multiples.
gt = wt ∗ et et = gt ∗ w-1t
=> Inverse Filtering

However, in general the function wt is not known. Because of that it is not easy to obtain the
inverse function w-1t.
To obtain a good approximation one can use a socalled “Optimum-Filter” or Wienerfilter.
An option is to reconstruct the waveform wt using the Autocorrelation function. The auto-corre-
lation function contains all the frequency information of the original waveform, but none of the
Version (15.02.2005)
phase information. The necessary phase information comes from the minimum-phase
assumption.
Input-Function ∗ Filter = Output-Function

(known) (wanted) (known)
Another option is to determine a Filter operator, that from the signal produces the wanted signal
(e.g. a Spike, a Minimum Phase wavelet etc.). This is called the Wiener filter or least-squares
filter.
This results in a system of equations:
g0 f0 g0 f0 = y0
g1 x f1 = g1 f0 + g0 f1 = y1
… … …
Solving this system of equations yields the appropriate filter operator f.
According to the aim there are different types of deconvolution:

(1) Spiking Deconvolution: desired output function is a spike (also whitening deconvolution)
(2) Predictive Deconvolution: attempts to remove the effect of multiples
∗
Version (15.02.2005)
14. Velocity analysis and NMO-Correction
Until now we have only discussed data processing methods that improve the signal of each
separate trace. We will now sum different traces, also called stacked, to improve the signal-to-
noise ratio and to decrease the amount of data which will be processed to obtain an image of the
subsurface. Before the stacking, a certain correction is applied on the different traces by
carrying out a velocity analysis.
A good velocity model is the basis for :

• Stacking (Improvement of S/N-Ratio)
• Appropriate conversion from traveltime into depth
• Geometrical Correction (Migration)
14.1 Normal-Moveout (NMO) Correction

Principle:
NMO-
Correction
Reflectionhyperbolas horizontal Alignment Stacking

Principle of NMO-Correction. The Reflections are alligned using the correct velocity, such that
the events are horizontally. Then all the separate traces are stacked (summed).
The traveltime curve of the reflections for different offset between source and receiver is calcu-
lated using:
2
2 2 x
t = t 0 + ------------
2
-
v stack
From this formula the NMO-correction can be derived and is given by:
2
2 x
∆t = t 0 – t ( x ) with t(x) = t 0 + ------------
2
-
v stack
The Moveout ∆t is the difference in traveltime for a receiver at a distance x from the source and
the traveltime t0 for zero-offset distance.
The NMO-Correction depends on the offset and the velocity. In contrast to the static correction,
the correction along the trace can differ. The NMO-correction is also called a dynamic
correction.
Version (15.02.2005)
To obtain a flattening of the reflections, the velocity must have the correct value. When the
velocity is too low, the reflection is overcorrected; the reflection curves upwards. When the
velocity is too high, the reflection is undercorrected; the reflection curve curves downwards.
NMO-Correction of a Reflection. (a) Reflection is not corrected; (b) with proper Velocity; (c)
Velocity is too low; (d) Velocity is too high.
Remark:
Low velocities have a stronger curvature then high velocities.
14.2 Methods for Velocity analysis.

The aim of the velocity analysis is to find the velocity, that flattens a reflection hyperbola, which
returns the best result when stacking is applied. This velocity is not always the real RMS
velocity. Therefore, a distinction is made between:
• vstack: the velocity that returns the best stacking result.
• vrms: the actual RMS-velocity of a layer.
For a horizontal layer and small offsets, both velocities are similar. When the reflectors are
dipping then vstack is not equal to the actual velocity, but equal to the velocity that results in a
similar reflection hyperbola.
There are different ways to determine the velocity:

• (t2-x2)-Analysis.
• Constant velocity panels (CVP).
• Constant velocity stacks (CVS).
• Analysis of velocity spectra.
For all methods, selected CMP gathers are used.
Version (15.02.2005)
(t2-x2)-Analysis
The (t2-x2)-Analysis is based on the fact, that the Moveout-expression for the square of t and x
result in a linear event. When different values for x and t are plotted, the slope can be used to
determine v2, the square root returns the proper velocity.
Example of a t2-x2-Analysis.
CVP - “Constant velocity panels”

The NMO-correction is applied for a CMP using different constant velocities. The results of the
different velocities are compared and the velocity that results in a flattening of the hyperbolas
is the velocity for a certain reflector. .
CVS - “Constant velocity stacks”

Similar to the CVP-method the data is NMO-corrected. This is carried out for several CMP
gathers and the NMO-corrected data is stacked and displayed as a panel for each different
stacking velocity. Stacking velocities are picked directly from the constant velocity stack panel
by choosing the velocity that yields the best stack response at a selected event.
CVP and CVS both have the disadvantage that the velocity is approximated as good as the
distance between two test velocities. Both methods can be used for quality control and for
analysis of noisy data.
Velocity-Spectrum
The velocity spectrum is obtained when the stacking results for a range of velocities are plotted
in a panel for each velocity side by side on a plane of velocity versus two-way travel-time. This
can be plotted as traces or as iso-amplitudes. This method is commonly used by interactive
software to determine the velocities.
Different possible methods can be used to determine a velocity spectrum:
Version (15.02.2005)
• amplitude of stacking
• normalised amplitude of stacking
• Semblance
•
Amplitude of Stacking
n
st = ∑ w i, t
i=1
where n=number of NMO corrected traces in the CMP gather; w=amplitude value on the i-th
trace at twoway time t.
Normalised Amplitude of stacking

st
ns t = ---------------------
n
∑ w i, t
i=1
Semblance
∑ st
2
1
Semblance = --- ⋅ ----------------------
t
-
n
∑ ∑ i, t
2
w
t i
Semblance-Calculations are only used for velocity analysis, because it returns always a value
between 0 and 1.
14.3 Problem of “Stretching” of the data caused by NMO correction

NMO is a dynamic correction, that means that the values of a single trace are shifted with
different amounts. This results for larger offets in a stretching of the data and an artificial
increase of the wavelength occurs.
Version (15.02.2005)
This effect is relatively large for horizontal reflections with low velocities. To reduce the effect
of the stretching on the result of the stacking procedure, the part with severe stretching of the
data is muted from the data (“stretch-mute”).
Dynamic correction results in a stretching of

the data, which results in a artificial increase
of the wavelength.
14.4 Factors influencing velocity estimates

The accuracy of the velocity analysis is influenced by different factors:
• Depth of the Reflectors
• Moveout of the Reflection
• Spread length
• Bandwidth of the data
• S/N-Ratio
• Static Corrections
• Dip of the Reflector
• Number of traces
By a combination of CMP’s that lie close together (Super gather), the accuracy is increased
when a small number of traces per CMP are available (low coverage).
Errors due to dipping layers and unsufficient static corrections can be reduced (DMO and
Reststatics, are discussed later on).
Version (15.02.2005)
15. Stacking
Stacking is perfomed by summation of the NMO-corrected data. The result is an approximation

of a zero-offset section, where the reflections come from below the CMP position. (For a
dipping layer, the reflections do not exactly come from below the CMP.
15.1 Muting
Sometimes, the data contains still noise signals, that influence the stacking.
Gelegentlich enthalten die Daten auch nach der Bearbeitung noch Störsignale, welche die
Stapelung beeinflussen. These traces are muted before the stacking.
Typical Noise signals are:

• Refractions (first breaks)
• Surface waves
• Air wave
Three possible muting procedures can be carried out:

• top mute: A certain timewindow in the beginning of the trace is muted (first breaks)
• bottom mute: A certain timewindow in the end of the trace is muted (surface waves)
• surgical mute: A time window in the middle of the trace is muted (air wave)
15.2 Methods of Stacking

Several methods can be used to combine the different NMO-corrected traces. The most
important are:
Mean stack
All NMO-corrected traces are summed and devided by the number of traces:
n
1
s t, mean = --- ⋅
n ∑ w i, t
i=1
Weighted stack
In certain situations, unequal weighting of the traces in a gather may yield results that are better
than the CMP stack. For example when certain traces contain a lot of noise. This type of
stacking is often used to suppress multiples by weighting the large-offset data more heaviliy
than the short-offset traces, because the difference in NMO between primaries and multiples is
larger for larger offsets. A weight factor α is introduced.
n
1
s t, meanweighted = --- ⋅
n ∑ αi wi, t
i=1
Diversity stacking / Min-Max-exclude
Certain traces are muted and not included in the stacking procedure.
• When certain values differ too much from the average value they can be muted (diversity
stacking). This to reduce the influence of spikes
• Exclusion of traces with the minimum and maximum amplitudes in the stacking procedure
(min-max-exclusion or alpha-trimmed stack).
Version (15.02.2005)
16. Residual static
In an early stage the general static corrections were discussed. Especially, the correction for the
topography and the influence of the weathered layer were discussed. The aim of the static
corrections is to shift individual traces in such a way that the reflections in a common midpoint
gather lie as accurate as possible along a hyperbola.
Topographic corrections and refraction statics solve this problem only for a certain part. Most
of the times small shifts between traces remain. To correct for these small shifts the residual
static correction is applied.
16.1 Principle of Residual statics

The process of residual statics consists of shifting the separate traces in such a way that the
optimal reflections are obtained. To make sure that the traces of a single CMP are not shifted
randomly, the shift is devided in a value for the source (“source static”) and a value for the
receiver (“receiver static”). For each source and receiver a value is determined. All traces with
a certain source are corrected with the value for that source. Similarly all traces with a certain
receiver are corrected with the value for that receiver.
The resulting shift (static correction) of a trace consists of the correction value of the source and
receiver of the corresponding trace.
This processing still assumes that the static shifts are caused by the interface. Therefore, this
processing is also called surface consistent static correction.
Scheme of residual static corrections

CMP-sorted Data
1. Velocity analysis
NMO-Correction
Residual correction
inverse NMO-Correction
2. Velocity analysis
NMO-correktion with new Velocities
Stacking
Version (15.02.2005)
17. Special Filters
In the proceding chapters the most important processing steps for reflection seismics are
discussed. Several other special filters and methods can be use to enhance the signals and to
suppress noise. Two commonly used filter techniques are:
• f-k filter
• τ-p filter
Both filters use more traces at once. e.g. a whole common shot gather (CSG) or CMP (two- or
three-dimensional transformations). The filters discussed before (deconvolution, frequency
filter etc.) are applied on separate channels (one-dimensional transformation)
17.1 f-k filter

Most of the time, it is difficult to separate the reflections from noise. This separation can be
made easier when the data is not processed in space-time domain, but transformed into another
domain.
One channal (1-D-Transformation)

t-domain f-domain
Time domain Frequency domain
Fouriertransformation
More channels (2-D-Transformation)

t-x domain f-k domain
Time-space 2-D-Fourier Frequency-Wavenumber
t-x domain τ-p domain

Time-space (Intercept - slowness)
τ-p-Transformation
What is f-k?
The f-k transformation is in principle a two-dimensional Fourier transformation. Corresponding
to the transformation of the time-axis to the frequency domain, the x-axis is transformed to the
wavenumber domain. The frequency indicates the number of oscillations per second. The
Wavenumber k indicates the number of wavelengths per meter along the horizontal axis (Some
authors define k as the number of wavelenths per meter along the horizontal axis times 2π). For
waves which propagate horizontally, the transformation returns the actual wavenumber. For
waves that do not propagate horizontally, the horizontal component of the wave is transformed.
An apparent wavelength and an apparent velocity is obtained:
v
v a = -----------
sin α
λ
λ a = -----------
sin α
with α=angle of the wavefront with the interface (or the angle of the ray with the vertical).
Version (15.02.2005)
horizontal Wave
v
vapp = v / sin α
entering Wave
α
v
Wavefront
(Reflections that travel vertically, reach the geophones at the same time and therefore have an
infinite apparent velocity)
The relation between the frequency and the wavenumber is given by: f = va k, i.e. the slope of
a line in f-k domain is the apparent velocity va.
f-k spectrum
The plotting of a dataset in f-k domain is called a f-k spectrum. (Analogous to the frequency
spectrum for the one-dimensional transformation from time to frequency).
The signals are separated and plotted as function of the frequency and slope (apparent velocity)
Negative wavenumbers indicate a slope in the other direction.
Applications of the f-k filter

Using the fk-transformation it is possible to mute a certain part of the f-k spectrum (for example
the part that contains the interface waves). When the data is transformed back to the t-x domain,
the interface waves are removed from the original data.
The spectrum that is muted is often defined as a certain fan in which the slopes or velocities are
muted. (Such Filters are called “fan-filter”,”pie-filter”, “dip-filter”, Velocity filter or Moveout-
Filter)
Version (15.02.2005)
Typical Applications:
• Suppression of noise signals with specific slopes (Interface waves)
• Suppression of multiple reflections
• Elimination of Artefacts in stacked Sections (post-stack)
Problem of spatial Aliasing

Aliasing in time domain occurs when the sampling rate of the signal is not high enough. A
similar aliasing effect can occur when the apparent wavelength is smaller than twice the
distance between the geophones. The spatial Nyquist-criterium is given by:
1
k N = ----------
2∆x
where ∆x=distance between Geophones.
(kN depends on the apparent velocity and the frequency of the signal)
Spatial Aliasing can be observed in the f-k spectrum when events that are present for positive
wavenumber continue on the other side of the spectrum for negative wavenumber values.
This effect is called “wrap-around”.
The Signal in shaded area A continues for wavenumbers larger than kN. The Aliasing occurs
when the data for wavenumbes larger than kN appear in the spectrum as area B.
Note that:
Similar to the ringing effects that occur when a frequency filter is used with wrong parameters
artefacts can also appear in the data when f-k filtering is applied.
17.2 τ-p-filter
An other filter that is often used is te τ-p filter (also called tau-p filter, Radonfilter or Slant-
Stack). This filter transforms the data from (t-x) into a domain of Intercept-time τ (t0-time) and
Slowness p (p~1/v). The relation between t-x and τ-p is given by
Version (15.02.2005)
τ = t 0 – px
Each p-value indicates a certain slope in t-x-domain. The energy along a line is summed. The
point of intersection with the t-axis (x=0) gives the intercept-time τ. In this way lines in t-x
domain become a point in τ-p-domain and reflections become ellipses.
Different Elements in t-x-domain

and the corresponding elements
in τ-p-domain:
A, C: part of a reflection;
B: line;
D: Reflection.
Typical applications of the tau-p Filter

• Velocity filter
• Time-depemdent velocity filter
• Suppression of multiples
• Interpolation between traces
• Analysis of Guided Waves
Version (15.02.2005)
18. Migration
Stacking returns the first image of the subsurface. However, for a complex geometry and
dipping reflectors, this image does not resemble with the reality. For example, the stacked data
can still contain diffraction hyperbolas. The process that corrects for these effects is called
Migration (also called “imaging”).
Example of a seismic Section. (a) Stacking without Migration. (b) with Migration.
Version (15.02.2005)
18.1 Geometrical Distortion

Zero-offset traces are generated by NMO correction followed by stacking. (All following
processing steps assume Zero-Offset-Data). The Data are plotted as if the reflection point is
present directly below the CMP. In reality, the reflection for a zero-offset measurement is
incident perpendicular on the layer. When a layer is horizontal, both facts are true. However,
when a dipping reflector is present, the point of reflection is present besides the point directly
below the CMP. All possible reflection points lie on a semi-circle that has a radius that depends
on the traveltime.
Positions of all possible reflection points with equal traveltime.
Typical Structures, that cause geometrical distortion, are:

• Dipping Reflector
• Valley
• Point reflector
Dipping Reflector
For a zero-offset measurement, the reflections coming from a dipping reflector travel perpen-
dicular to the dipping interface. However, they are plotted in a stacked section as if they have
travelled perpendicular to the surface. This is why the image of a dipping reflector obtaines a
wrong dip in the stacked section.
The difference between the real dip and the dip in the stacked section is given by:
sin α real = tan α Stapelung
Version (15.02.2005)
A perpendicular reflector will be plotted with a dip of 45° abgebildet. This shows that the
maximum dip of Reflections in a Stack is 45°. Larger dips are thus due to noise signals or
other effects.
Syncline
Another effect that occurs often in stacked data is for example a syncline, a valley in stratified
rocks in which the rocks dip toward a central depression.
As shown in the figure, there are different rays coming from position B that are perpendicular
to the reflector and thus are measured by a zero-offset measurement. The different reflections
have different traveltimes, so in stead of only one reflection, three reflections are measured at
position B and thus three reflections are plotted at position B. In a stacked section, the syncline
is not directly distinguishable. In stead, we see a “bow-tie”.
Point reflector
Point reflectors appear in a stacked section as a diffraction hyperbola. This hyperbola becomes
visible, because all rays from all directions are perpendicularly to this point reflector and will
result in a reflection. The traveltime increases with increasing distance. The travel time curve
that results is a hyperbola (Diffraction hyperbola).
Diffraction hyperbolas also appear at edges.
Version (15.02.2005)
18.2 Methods for Migration

Migration is a proces that reconstructs a seismic session so that reflection events are reposi-
tioned under their correct surface location and at a corrected vertical reflection time.
The following effects is corrected for:

• point diffractions are collapsed to one point
• Location and dip of layers are adjusted
• Improvement of the resolution by focusing of the energy
Basic Corrections
The real angle of a dipping reflector can be obtained by using the equation given earlier.
However, the actual position of the dipping layer is not obtained.
A simple graphical reconstruction using arcs can be performed to migrate the data.
One draws a semicircle with the radius of the travel time through a reflection point. The real
position lies somewhere on this circle.When different semicircles are drawn for different points
on a dipping reflector, then the real position of the dipping layer can be obtained by the tangent
of these circles.
Point diffractions can be reconstructed similarly. They are present at the apex of the diffraction
hyperbola and can be obtained by the different arcs that are drawn with the radius of the travel
time from different points on the hyperbola.
Methods for Migration

There are different ways to migrate seismic data:
Version (15.02.2005)
• “wavefront charting”
This is in principle the method discussed before.
• Diffraction-Migration (Kirchhoff-Migration)
All energie is added along diffraction hyperbolas.
• Fk-migration
Correktion for slopes in the Fk-domain
• Downwards continuation
Operation that corrects for the propagation of the wave fronts.
(e.g. phase shift migration)
• wave-equation migration(FD-Migration)
Correction for the traveltime by solving the wave equation
The different methods have also different properties and differ in:
• Accuracy and type of und type of the required velocity model.
• Vertical velocity change can be taken into account
• Lateral velocity change can be taken into account
• Correction of dip
• Calculation time
Velocity model for Migration

Some methods only use one velocity for the whole dataset, whereas other methods can use a
complex velocity model with vertical and lateral velocity changes.
In general one can say:

A complex structure in the subsurface requires a more complicated velocity model.
For a starting model, often the model is used that is employed for the stacking. However, most
of the time adaptions are needed (The migration velocity for deeper layers is often between 90%
and 80% of the stacking velocity).
When the choosen velocity is too small, the diffractions are not fully collapsed into one point.
The data is undermigrated. When the choosen velocity is too high, then arc appear above the
diffraction point (“Smile”). The data is overmigrated.
Boundary effects
At the boundary of a seismic section arise often strong effects : half arcs or hyperbolas.
18.3 Special Migration and extension
3D Migration
Until now we only discussed two-dimensional data. In reality, the structures are not always
directly beneath the profile. There are often reflections from structures that lie off-line from the
seismic profile (“side swipe”). Reflections from such off-line structures can not migrated to
their correct position by using 2D migration. To correct these structures correctly, the data must
be processed and migrated in three-dimensions.
Version (15.02.2005)
Depth migration
Until now the migration is treated mainly in time domain. However, one is also interested in the
actual depth of a reflector. There are different possibilities to reach this goal:
• Change of stack in depth
• Change of migration in depth
• Depth migration
When a stacked section is converted then all geometrical distortions remain.

When a time migrated section is converted into depth then the timing of the reflections are
converted into depth using the vertical velocity model.
However for complex structures and strong velocity contrasts distortions appear, because in
reality the waves do not propagate only vertically. Different algorithms consider the travel path
through a complex structure and transform directly the stacked sections into a depth section. The
conversion from travel time into travel path is not carried out only vertically, but along the
(correct) direction of the traveling waves. Depth migration converts a seismic time section into
a depth section where the perpendicular travel time of each individual reflector is converted in
depth using the “correct” velocity model. This is the most accurate way of migration, but
requires also the most processing time.
“Dip Moveout Correction” (DMO-Correction)

(sometimes also called partial Prestack-Migration)
Until now we have assumed that the data was real zero-offset data, that is obtained by stacking
of multi-offset data. It is mentioned several times that the NMO correction for strong dips does
not result in the correct result. The traces of a CMP have a different reflection point and the apex
of the reflection hyperbola is horizontally shifted. To correct for these effects additional
processing is needed: DMO correction (“Dip Moveout”).
Version (15.02.2005)
Prestack-Migration
One step further is to use all data which is available before the actual stacking is carried out.
That means that all geometrical distortions are removed before the stacking.
=> Prestack-Migration
In principle the same correction is applied on the data as the use of DMO followed by Post-stack
migration.
DMO + Poststack-Migration = Prestack-Migration
Poststack migration migrates all stacked CMP traces, whereas prestack migration migrates all
traces of all shots. Because there are a lot more traces before the stacking compared with after,
the prestack migration requires much more computing time.
Version (15.02.2005)
19. Post processing
After succesfull stacking and migration, the data are often “post processed“ for reproduction,
e.g. interpretation purposes
Aim of post processing:

• Elimination of filter effects
• Eliminate effects that originate from stacking or migration
• Adjust Amplitudes for plotting purposes
• Improve coherency of Reflections
Frequently used operations in “Postprocessing”:

• Time dependent Frequency filter
• Deconvolution
• Coherency filter
• AGC
19.1 Time dependent frequency filter
Reason for filtering after Stacking:

• By using different filters additional noise is generated
• The frequency content of the data has changed due to different processing steps (e.g. NMO
stretching)
• Time dependent frequency filter is often not used before the stacking procedure
Time dependent frequency filter consists of filter parameters that change with travel time. In this
way, the absorption of high frequencies with increasing travel time can be adjusted for.
The filter parameters do not vary continuously, but the travel time is divided in different
windows for which the different parameters are determined
19.2 Deconvolution
Sometimes, it is tried to improve the resolution of the data or to change the waveform to a
predefined shape after stacking using deconvolution.
The deconvolution after stacking does not differ from the deconvolution before stacking.
Parameters and the working of the operations are equal.
19.3 Coherency filter

After stacking or Migration, it is possible that the data are disturbed by noise and artefacts.
Using special filters it is possible to improve the signal to noise ratio and the continuity of the
reflections.
Often used filter methods are

• fk-Filter
• fx-dekonvolution (Wiener-Levingson Filter)
• Korhunen-Loeve Filter (Eigenvector-Filter)
Version (15.02.2005)
All these Filters use several seismic traces in the calculation. (More channel-filter)
19.4 Adaption of Amplitudes

Depending on the aim of the interpretation or output, the reflections should be gained such that
also weak reflections can be distinguished
When one is interested more in structures, then all reflectons are amplified using an AGC gain
with a small window. When the reflection characteristic, the difference in amplitude, should be
preserved then an AGC gain with large window is used.
Adjustments of the amplitudes is also often needed when the data are plotted on a printer.
19.5 Archiving of Data

When all processing steps are carried out, the data is plotted or saved. This can be done by
• digital archiving and Visualisation on a computer screen
• Plotting on paper
Digital archiving
• Archiving on Magnetic tapes or CD-ROMs
• Data can be reproduced using different processing parameters or different time slices
• Interpretation on screen with an interactive program
Often the Data are archived in SEG-Y Format, which is generally used everywhere.
The data can be easily exchanged with other programs.
Printing on paper
• Safe archive
• Interpretation by hand (without expensive software-package)
• Improved view for long profiles
Components of a paper print:

Except for the data, also additional information should be mentioned on the paper print, such
that everybody knows what processing was performed on the data. It is e.g. important whether
the data is only stacked or also migrated. Such information is mentioned in the socalled
“sidelabel”.
Version (15.02.2005)
Example of a “Sidelabel” and its components.
It is common, additional to the information in the sidelabel, to show also the parameters that
change along the profiel. This information is then plotted above the data (“Header Plots”).
In the header plot different information can be given, e.g.:

• Velocity functions
• Fold
• Topography
• Point of intersection with other profiles
• Location of borehole loggings
Version (15.02.2005)
20. Interpretation
In this lecture only a general overview of the interpretation of seismic data is given.
Further information can be obtained for example from the course “Seismische Faziesanalyse”,
that is usually given in the Summer semester.
One can distinguish different aims for the processing of the data:
• Mapping of geological Structures
• Seismostratigraphy
• Seismic Facies analyse
• Modeling
In addition there are special interpretation possibilities for 2D and 3D Data.
20.1 Mapping of geological Structures

For the mapping of geological structures in seismic data, the position of main horizons and
disturbances are mapped and the form and position of faults are identified.
Aim:
• to obtain the geological profile
• to obtain depth charts of horizons and disturbances
Foregoing steps:
• “Picking” of the beginning of Reflections
• Determination of their travel time
When more seismic lines are available:

• Draw travel times in a map
• Generate isolines
For more crossing profiles, a horizon on a certain profile should be tracked to the other profiles
to obtain a closed curve. In this way, one can check whether the same phase is correctly tracked
or that one has followed a wrong reflection after a disturbance.
“Fence-Diagram”
Three-dimensional picture of crossing profiles
Version (15.02.2005)
20.2 Seismic Sequence-Analysis (seismic sequence stratigraphy)
Aim
• Identification of sequences in seismic data
• Determine the sequence of the sedimentation in time
• Analyse sealevel fluctuations.
Sequence of a coast.
The fundaments of sequence-stratigraphy were investigated in the 70” by a research group of

the oil company Exxon. Different concepts were defined which are now standard to describe
seismic data.
An elaborate describtion can be found in the articles from the Exxon-Group. Especially
Mitchum, R. M., Jr. and P. R. Vail and S. Thompson, III (1977). Seismic stratigraphy and global changes of sea
level, part 2: The depositional sequence as a basic unit for stratigraphic analysis. In: Payton (edt.), Seismic
stratigraphy - applications to hydrocarbon exploration. Am. Assoc. Pet. Geol., Memoir, 26, Tulsa, Oklahoma,
53-62.
Vail, P. R. and R. M Mitchum, Jr. and R. G. Todd and J. M. Widmier and S. {Thompson, III} and J. B. Sangree and
J. N. Bubb and W. G. Hatlelid (1977). Seismic stratigraphy and global changes of sea level. In: Payton (edt.),
Seismic stratigraphy - applications to hydrocarbon exploration. Am. Assoc. Pet. Geol., Memoir, 26, Tulsa,
Oklahoma, 149-212.
Sequences are on the upper- or/and lower part terminated by “unconformities” or concor-
dancen. The figure below shows the most important types:
Version (15.02.2005)
When different sequences can be identified in a seismic section, one can construct the time
sequence of the sediment (Chronostratigraphy). One can then draw conclusions about the
interpretations for different phases of relative rise or fall of the sealevel (Transgression und
Regression).
20.3 Seismic Facies Analysis

(Only a short overview is given)
Additional to the boundaries of a seismic sequences, one can also investigate the reflection
characteristics inside a sequence. Areas with similar refection character correspond to a seismic
facies. Not only the time sequence of the sedimentation can be obtained, but it is also possible
that conclusions can be drawn about the sedimentation in the environment.
Similar to the boundaries of a sequence, there are different concepts, to describe the character
of reflections.
Geometry of reflections in a sequence (Selection).
Additional to the description of reflections in a sequence also different shapes of seismic facies
exist:
Three-dimensional shapes of seismic facies units.
Version (15.02.2005)
20.4 Interpretation of 3D-data

For a three-dimensional interpretation we need three-dimensional data. Special systems are
needed for a usefull analysis and visualisation of 3D-Data.
Advantage of 3D-Interpretation
• Arbitrary slices through data is possible
• Horizontal time slices can be generated and interpretated
• Combination of different plotting methods
• Analysis of the actual dip and position of structures
Visualiation methods of three-dimensional Data

• Data-Cube
• Arbitrary slices
• Horizontal time slices
• “Chair”-Diagram
20.5 Other Aspects of interpretation
AVO - Amplitude variation with offset

Investigation of the angle-dependence of the reflections for different distancew between the
source and receiver. A variation in amplitude depends on the Poisson-ratio, density and the
seismic velocity and can result in a conclusion concerning the oil and gas content.
One needs data with a high quality. Additionaly, the data processing must be adapted such that
the amplitudes are preserved
Attribute-Analysis
The common amplitudes are replaced by another parameter, e.g.:
• Instantaneous phase
• Instantaneous frequency
• Envelope
Using these attributes certain aspects of the data can show up more clearly.
Typical example for such an analysis is the emphasis and investigation of so-called “bright
spots”. These bright spots have exceptional strong amplitudes, that are due to the high
Impedance contrast of oil or gas with water.
Interpretation of “deep” reflections

In the processing of seismic Data from the deep Crust, one obtains in general no continous
reflections, but patterns with different reflectivities. In such cases the strong phases are often
plotted as lines and interpreted.
Version (15.02.2005)
21. Other related methods
Finally, some methods and concepts are discussed that are not discussed thoroughly in the
lecture, but have to do with reflection seismics.
21.1 VSP - “vertical seismic profiling”

For a standard seismic measurement, the source and receiver are always present on the surface.
However, it is also possible to carry out measurements where either the receiver or the source
are present in a borehole (In general, the receiver is present in the bore hole).
One can distinguish:

• zero-offset VSP: The source is present at the surface close to the bore hole.
• offset VSP: The source has a certain distance from the bore hole.
• walk-away VSP: Measurements are carried out with different distances between the source
and the bore hole .
Typical offset VSP measurement.
In general, the VSP measurements are expensive (one needs a bore hole), but there are also
certain advantages compared with the standard seismic measurement.
• Depth of the receiver is known

-> Improved velocity depth model.
• Smaller travel times
-> Improved resolution, less attenuation
• Improved separation between primaries and multiples
• Improved deconvolution due to the measurement of the direct wave.
In a VSP seismogram one can distinguish between upgoing and downgoing waves. The direct
wave (downgoing wave) enables the determination of a velocity model. Reflections are upgoing
waves.
A separation of up- and down-going waves is possible by using a fk-filter.
Version (15.02.2005)
Synthetic Seismogram of a VSP-Measurement

(Kearey and Brooks, 1991).
The reflections from the upgoing waves can also be aligned horizontally by shifting with the
traveltime of the direct wave. This data can then be stacked and a seismic trace is obtained that
can be compared with a general surface recorded seismic trace.
(a) Synthetic VSP section with downgoing waves removed by fk-filter; (b) Correction with the
traveltime of the direct wave ; (c) Stacked seismogram produced by stacking in the shaded
corridor zone of (b) to avoid muliple events.
Version (15.02.2005)
Comparison of the zero-offset VSP stacked

section (last figure (c)) and a conventional
seismic section based on surface profiling.
21.2 Cross-hole seismic

In stead of using only one borehole, also more boreholes for sources and receivers can be used.
This method is called Cross-hole seismics.
In general, only the travel time between different source-receiver combinations is used to obtain
a velocity model of the subsurface between the boreholes (Tomography).
Version (15.02.2005)
21.3 Sidescan sonar

Acoustic measurement methods like the echosounder or Sidescan Sonar use a similar principle
as reflection seismic. Both methods are used to investigate the bottom of the sea.
A sidescan sonar insonifies the sea floor to one or both sides of the survey vesssel by beams of
high frequency sound (30-110 kHz) . Sea bed features facing towards the survey vessel, such as
rock outcrops or sedimentary bedforms, reflect energy back towards the transducers while in the
case of features facing away from the vesser, or a featureless sea floor, the acoustic energy is
reflected away form the transducers..
21.4 Georadar
The Georadar uses an electromagnetic signal in stead of elastic waves. Often, only one source
and one receiver are used with a common offset. A radargram contains reflections which are
comparable with seismic results. The processing (filtering, amplitude correction, NMO and
Migration) are similar.
21.5 New developments
4D-Seismics
(time-lapse-seismics)
Repetition of 3-D seismics at a certain location at different times to analyse the changes in the
subsurface. Is mainly used in the oilindustry.
3-C Seismic (3-components Seismic)

Simultaneously measureming of x-, y- and z-components of the seismic wave fields.
4-C Seismic (4-components Seismic)

For marine measurements, regular 3-components geophones are placed together with a hydro-
phone on the bottom of the ocean.
Converted waves
Analysis and processing of P-S or S-P converted Waves, to obtain additional information of the
subsurface.
Version (15.02.2005)
22. Selection of used References (Figures etc.)
Brouwer, J. and Helbig, K. (1998). Shallow high-resolution reflection seismics. Handbook

of Geophysical Exploration, Volume 19. Elsevier Science, Amsterdam, 391 pp.
Domenico, S.N. and Danbom, S.H., 1987. Shear-wave technology in petroleum exploration
- past, current and future. In: Danbom, S.H. and Domenico, S.N. (eds.), Shear-wave
exploration. Society of exploration geophysists, Tulsa, OK, USA.
Emery, D. and Myers, K.J. (edt.) (1996). Sequence Stratigraphy. Blackwell Science, Oxford,
UK, 297 pp.
McQuillin, R., Bacon, M. and Barclay, W. (1984). An introduction to seismic introduction -

Reflection seismic in Petroleum Exploration. Graham and Trotman, London, UK, 287 pp.
Meissner, R and R. K. Bortfeld, 1990. DEKORP-Atlas. Springer-Verlag. Berlin.
Sheriff, E.G. and Geldart, L.P. (1995). Exploration Seismolgy, (2nd ed.). Cambridge
University Press, Cambridge, 592 pp.
Yilmaz, Ö. (1987). Seismik data prosessing. SEG Tulsa, OK, 826 pp.
Kearey, P. and Brooks, M. (1991). An inrtoduction to geophysical prospecting. Blackwell

Scientific Publications, Oxford, 254 pp.
Reynolds, J.M. (1998). An introduction to applied and environmental geophysics. John

Wiley and sons, Chichester, UK, 796 pp.
Telford, W.M., Geldart, L.P. and Sheriff, R.E (1990). Applied Geophysics (2. ed.).
Cambridge University Press, Cambridge, 770 pp.
Version (15.02.2005)
Version (15.02.2005)

Reflection

Uploaded by

Copyright:

Available Formats

Reflection

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Reflection

Uploaded by

Copyright:

Available Formats

Version (04.02.

Yilmaz, Ö. (2001). Seismic data Analysis: Processing, Inversion, and Interpretation of

Interpretation of seismic data:

McQuilin, R. et al. (1986). An introduction to seismic interpretation. Graham and Trotman

Applied Geophysics in general (with one part about Seismics):

Kearey, P. and Brooks, M. (2002). An introduction to geophysical prospecting. Blackwell

Reynolds, J.M. (1998). An introduction to applied and environmental geophysics. John

Knowledge of the lectures “Allgemeine Geophysik I und II” is indispensable.

1.1 What is Reflection seismics?

Source Receiver Raw data

geological subsoil seismic Section

1.2 Historical Developments

Table 1: Historical development of reflection seismics(Sheriff and Geldart, 1995)

1.3 Use of Reflection seismic

Table 2: Application of reflection seismic tools

Aim Penetration Depth

Measurements in Echo sounder, ~0m

(1) Body waves:

Propagation of different waves

2.1 Principles of wave propagation

Combining these two equations we obtain the Acoustic wave equation:

w(t) is the source signal and v is the wavespeed which is given by

2.2 Interface: Reflection, Refraction, Conversion.

When a wave encounters an interface three phenomena can occur:

Law of Snellius sin α

2.3 Reflection- and transmission-coefficients

The Product Z = vρ is introduced as the acoustic Impedance.

Angle-dependent reflection- and transmission-coefficients for P- and S-waves

2.4 Geometry of wave propagation: rays

=> geometrical Seismic

Picture of a typical travel time diagram for a two-layer case.

The important rays in the travel time diagram are:

The separate elements will be discussed in the following paragraphs:

Reflected wave in the horizontal two-layer case

The simplest case is the horizontal two-layer configuration.

Using pythagoras we have:

The crossing with the time-Axis (x=0) is the timezero, t0:

This formula will later be used in the processing of seismic data.

Reflections - horizontal layered medium

Reflection from a dipping reflector

Ray path and traveltime diagram for a dipping Reflektor.

x Horizontal two-layer case with the rays and the

• Change of traveltime with depth

The seismic velocity depends on the elastic parameters as follows:

From these equations we observe that:

Seismic velocities and the factors influencing it.

• Lithology (Matrix und Struktur)

Examples for typical velocities

material vp (km s-1)

Igneous / Metamorphic rocks

(Kearey and Brooks, 1991)

Direct measurements of the seismic velocity

4. Amplitudes and Attenuation

The amplitudes of a seismic signal are influenced by several factors:

Phenomena causing the degradation of a seismic wave (Reynolds, 1998).

Spherical divergence (geometrical spreading)

When a seismic wave propagates, the energy is spread along

• plane Wave: konstant r2