0% found this document useful (0 votes)
310 views

Handout Processing 1

This document provides an overview of seismic data processing techniques. It discusses improving the signal-to-noise ratio through techniques like stacking and filtering. The document then describes the typical processing sequence, including steps like geometry assignment, editing, muting, static corrections, velocity analysis, normalization, stacking, and migration to obtain a high-resolution image of subsurface reflectors and velocities. Finally, it discusses various static correction methods used to compensate for travel time variations from factors like topography and shallow low-velocity layers.

Uploaded by

Mordekhai
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
310 views

Handout Processing 1

This document provides an overview of seismic data processing techniques. It discusses improving the signal-to-noise ratio through techniques like stacking and filtering. The document then describes the typical processing sequence, including steps like geometry assignment, editing, muting, static corrections, velocity analysis, normalization, stacking, and migration to obtain a high-resolution image of subsurface reflectors and velocities. Finally, it discusses various static correction methods used to compensate for travel time variations from factors like topography and shallow low-velocity layers.

Uploaded by

Mordekhai
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 39

1

SEISMIC
DATA
PROCESSING
(PART I)
2
Improve the signal to noise ratio: e.g. by measuring of
several channels and stacking of the data (white noise
is suppressed).
Obtain a higher resolution by adapting the waveform
of the signals.
Isolate the wanted signals (isolate reflections from
multiples and surface waves).
Obtain information about the subsurface (velocities,
reflectivity etc.).
Obtain a realistic image by geometrical correction
(Conversion from travel time into depth and correction
from dips and diffractions).
Processing of Reflection Data
3
Seismic Processing
4
Usually geared to a particular type of application
Mostly CMP reflection processing;
Land or marine, 2D or 3D.
Commercial:
ProMAX (Landmark);
Omega (Western Geophysical, marine);
Focus (Paradigm);
Amoco and almost every other company
have their own
Vista (Seismic Image Soft)
Universities:
Stanford Exploration Project;
Seismic UNIX (Colorado School of Mines);
FreeUSP (Amoco);
SIOSEIS (Scripts, marine);
Seismic Processing Systems
5
Printing/storage
Seismic
Processing
Systems
6
CMP Processing Sequence
1. Demultiplex, Vibroseis Correlation, Gain Recovery.
Conversion from file formats produced by field data loggers into
processing:
SEG-Y, SEG-2.
ProMax, Focus, Omega, SU, Vista, etc., internal formats.
2. Field Geometry.
Assignment of source-receiver coordinates, offsets, etc. in
the trace headers.
3. Edit.
Removal of bad traces (noisy channels, poorly planted
geophones, channels contaminated by power line noise,
etc.).
7
CMP Processing Sequence (Continued)
4. First Arrival Picking.
May be semi-automatic or manual.
Required for generation of refraction static, models and for
designing the mutes.
5. Elevation Static.
Based on geometry information, compensates the travel-time
variations caused by variations in source/receiver elevations.
Transforms the records as if recorded at a common
horizontal datum surface.
6. Refraction Static.
Builds a model for the shallow, low-velocity subsurface.
Compensates travel-time variations caused by the shallow
velocities.
8
CMP Processing Sequence (Continued)
7. Top, Bottom, and Surgical Mute.
Eliminates (sets amplitude=0) the time intervals where strong
non-reflection energy is present: First arrivals, ground roll,
airwave.
8. Refraction Static.
Compensates geometrical spreading.
Based on a simple heuristic relation.
9. Trace Balance.
Equalizes the variations in amplitudes caused by differences
in coupling.
In true-amplitude processing, replaced with surface-
consistent deconvolution.
9
10. Deconvolution or Wavelet Processing.
Compresses the wavelet in time, attenuates reverberations.
Converts the wavelet to zero-phase for viewing.
11. Gather, CMP Sort.
12. Moveout (Radon, -p, f-k) Filtering.
Attenuates multiples, ground roll.
13. Velocity analysis.
For each of the CMP gathers, determines the optimal
stacking velocity.
CMP Processing Sequence (Continued)
10
14. Dip Moveout (DMO) correction.
Transforms the records so that the subsequent NMO + stack
work well even in the presence of dipping reflectors.
15. Normal Moveout (NMO) correction.
Removes the effects of source-receiver separation from
reflection records.
Transforms the records as if recorded at normal incidence.
16. Residual statics.
Removes the remaining small travel time variations caused
by inaccurate statics or velocity model.
CMP Processing Sequence (Continued)
Steps 13-16 above are usually iterated 3-5 times to produce accurate
velocity and residual statics models.
Success of velocity analysis depends on the quality of DMO/NMO and
residual statics, and vice versa.
11
CMP Processing Sequence (continued)
17. CMP Stack.
Produces a zero-offset section.
Utilizes CMP redundancy to increase the Signal/Noise ratio.
Can employ various normalization ideas, e.g., diversity
stack.
18. Migration.
Transforms the zero-offset time section into a depth image.
Establishes correct extents and dips of the reflectors.
19. Frequency Filtering and Display.
Attenuates noise.
Provides best display for interpretation.
12
Unspecified blocks contain mandatory or optional
information.
Initial Process
Schematic arrangement of a seismic field tape:
13
The principle of multiplexing is already discussed in the section
measurement system. It is used when the capacities of the AD converter are
not sufficient to digitize and save all channels at the same time. This is
common for older measurement systems or for measurements with a large
time window and a lot of channels per shot. The separate values of all
channels are sorted by samples and not by channels.
It is difficult to
process the data
in this form. It is
more convenient
and illustrative
when the data is
sorted by traces.
Demultiplexing (Sorting of The Data)
14
To deal with seismic data and their processing it is helpful to know
something about the data formats which are used to store these data.
Almost every program and every large oil company has developed their
own format. In the course of time several standards have been
developed, to make an exchange possible between the different
programs. The most important standards are:
SEG-Y
SEG-D
SEG-2
The other Formats SEG-D and SEG-2 are often used for the storage of
raw data. They are suitable for multiplex data and make it possible to
save traces with different lengths. e.g. the SEG2 format saves each shot
separately.
Demultiplexing (Sorting of The Data)
15
SEG stays for Society of Exploration Geophysicists. This is the most
important society of geophysicists and seismologists in the oil
industry.
The different formats are very similar. The most important standard
today is the SEGY format. Every file consist of several parts. The
standard describes which information is put where in the file :
Structure of The SEG-Y Format
16
Assignment of source-receiver coordinates, offsets, etc.
in the trace headers:
Determine Source and receiver position for
measured data.
Calculate CMP position.
Field Geometry
17
Often, signal traces must be examined visually to detect
and correct or reject those errorneously recorded or
exceptionally noisy. For this purpose a compressed plot
is made off all records. Example of the former are blank
or dead traces, reversed polarity, signal clipping,
incorrect trace sequence, and errorneous recording time.
Trace Edit
18
Editing
19
Editing
20
Muting avoids noise wave contamination of reflections on
summed (stacked trace).
Trace Mute
21
Static corrections are applied to seismic data to compensate for the effects of
variations in elevation, weathering thickness, weathering velocity, or
reference to a datum. The objective is to determine the reflection arrival time
which would have been observed if all measurements had been made on a
(usually) flat plane with no weathering or low velocity layer present. These
corrections are based on up-hole data, refraction first-breaks, or event
smoothing.
Static Correction
22
Static Correction
Static Correction : The whole trace is corrected with the same
time shift.
Objectives 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.
Dynamic Correction : Different time windows in the trace
are corrected differently. This results
in stretching and compression of the
events (e.g. NMO Correction).
23
Reference Diagram for Land Static
Corrections
24
Datuming
25
Topographic Correction (elevation statics).
Up-hole-correction using shots in borehole.
Refraction statics : Corrections for weathered layer
Delay-Time.
Generalized reciprocal method (GRM).
Diminishing residual matrices (DRM).
Methods for Static Correction
26
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.
Topographic Correction
27
Subdivision of time shift for source and receiver.
All traces with equal source are corrected for the time
shift of the specific source.
All traces with equal receiver are corrected for the time
shift of the specific receiver.
The statics correction is the sum of the corrections for
appropriate source and receiver.
Topographic Correction
28
Correction for the weathered (low velocity) layer (area
above water layer where pores are filled with air rather
than water).
Uphole-Static
When a shot is fired, also the travel time to the surface
is recorded and from this travel time, the velocity of the
weathered layer can be estimated.
29
Correction for the weathered layer.
Using the first breaks of a certain shot (refracted
energy) a model can be constructed for the weathered
layer (velocities and depth).
When the distance between the receiver is too large,
sometimes supplementary refraction measurements are
carried out.
Refraction Statics
30
Methods to determine the velocity and depth of the
weathered layer using refractions.
Delay-Time.
GRM (generalized reciprocal method).
DRM (diminishing residual matrices).
NO STATICS NO STATICS
WITH STATICS
WITH STATICS
32
Stack of Line without The Application of
Statics
33
Example Where ProMax Routine Made a
Static Solution FACTORS WHICH AFFECT AMPLITUDE FACTORS WHICH AFFECT AMPLITUDE
35
The amplitude of a seismic signal decays with
increasing travel time.
To obtain a realistic image, this decay must be
compensated for.
In general, it is difficult to describe this amplitude decay
analytically, so an approximation is usually made.
Amplitude Correction
36
Reflection and transmission at an interface.
Geometrical spreading.
Absorption.
Receiver response.
Measurement system
Loss of Amplitude
Loss of amplitude due to:
37
Individual large amplitudes dominate the processing.
Reflections are difficult to recognize.
Strong amplitude contrasts influence the digital filtering
(especially for large travel-times).
Problem for Data Processing
OBJECTIVES OF GAIN OBJECTIVES OF GAIN
STRUCTURAL PROCESSING
(AGC or TRACEWISE BALANCE)
Best continuty
Events visable at all times
STRATIGRAPHIC PROCESSING
(RELATIVE AMPLTUDES)
Bright spots visible
Amplitudes proportional to reflection
coefficients
GAIN 6
AMPLITUDE VARIATIONS AMPLITUDE VARIATIONS
SOLUTIONS: SOLUTIONS:
1. 1. Tracewise Tracewise balance balance
All times All times
All traces All traces
2. Relative amplitude 2. Relative amplitude
Systematic Systematic
Output proportional to reflection Output proportional to reflection
Coefficients Coefficients
GAIN FUNCTIONS VARY WITH: TIME, RANGE, LATERALLY GAIN FUNCTIONS VARY WITH: TIME, RANGE, LATERALLY
GAIN 10
40
Trace equalization.
AGC (Automatic gain control).
Correction for the spherical divergence.
Programmable gain functions
Methods to Preserve Amplitude Information
GAIN FUNCTION APPLICATION GAIN FUNCTION APPLICATION
GAIN 11
GAIN FUNCTION APPLICATION GAIN FUNCTION APPLICATION
GAIN 13
EXAMPLE OF LATERAL VARIATION EXAMPLE OF LATERAL VARIATION
GAIN 14
44
The simplest method is the normalization of the different
traces. All absolute values of a trace are summed and
compared with a reference value. A scaling factor is
determined from the difference between the summation
and the reference value, which is used to multiply all data
with.
Trace Equalization
45
Normalization of amplitude for a certain time sample in a
certain time window (not for the whole trace).
Advantage:
All traces are more equal which is needed for further
processing.
(Stacking: summation of different traces)
Amplification of Amplitudes for larger travel times.
Disadvantage:
Shadow effect.
Can lead to amplification of noise.
No physical base for amplification.
Automatic Gain Control (AGC)
46
Automatic gain controls scales by increasing the amplitudes in
segments of the trace where the amplitudes are low.
It does this by using a moving time window of around a second and
dividing the amplitude of the center point by the root mean-squared
(rms) amplitude of the window.
Some reflectors have stronger amplitudes than others and this may
obscure smaller, but important, reflections.
Automatic Gain Control (AGC) boosts the amplitudes so that they
are all of similar size.
Automatic Gain Control (AGC)
47
Automatic Gain Control (AGC)
48
t . v ) t ( G
r
) t ( A = =
1
[ ]
[ ] [ ] ) 0 ( / ) 0 ( / ) ( ) (
) (
1
) (
2
2 tw tw tw rms
tw tw rms
t t v t v t G
t t v
t A = =
Homogeneous space :
Layered space :
Correction for Spherical Divergence
49
Advantage:
Physical base for amplitude correction.
Relative Amplitude difference remains equal.
Disadvantages:
Velocity function not know beforehand.
Noise sources can still remain dominant.
Correction for Spherical Divergence
50
Programmable Gain Control
51
Calculation of decay of amplitude and determine a
Gain function
Programmed Gain Curve
52
Gain is time-variant scaling based on function, g (t). Base on some
criteria, this function is defined at the time samples (shown by solid
circles) that are usually at the center of specified time gates along the
traces as indicated by 1,2,3 and 4. Gain application simply involves
multiplying g (t) by the input trace amplitudes.
Applying Amplitude Corrections
53
Scale factors are indicated by the circled numbers at the
times of application
Four Different PGC Functions
54
Restored amplitudes at late times by correcting for geometrical
spreading (unfortunately ambient noise also has been strengthened
(Yilmaz, 1987).
Corrected Field Records From a Land Survey
55
Before AGC
56
After AGC
57
All traces are normalized using a certain amplitude:
RMS
Median value
Maximum Value
Advantage:
All traces are more equal which is needed for further
processing (Stacking: summation of different traces).
Disadvantage:
No physical base for amplification.
No equalization of losses with time.
Large value in a trace can dominate.
Trace Balancing
58
Common-shit gathers after trace balancing. Balancing is time-
invariant scaling of amplitudes to a common rms level for all traces.
59
Filtering
SPLIT SPLIT- -SPREAD FIELD RECORD SPREAD FIELD RECORD
TYPES OF NOISE TYPES OF NOISE
62
CAUSES OF POOR SIGNAL CAUSES OF POOR SIGNAL- -TO TO- -NOISE NOISE
NON NON- -OPTIMAL FIELD PROCEDURES OPTIMAL FIELD PROCEDURES
STRONG COHERENT NOISE STRONG COHERENT NOISE
SCATTERING OR ABSORPTION SCATTERING OR ABSORPTION
- - VULCANICS VULCANICS
- - FRACTURED ZONE FRACTURED ZONE
NEAR NEAR- -SURFACE PROBLEMS SURFACE PROBLEMS
- - POOR COUPLING POOR COUPLING
- - ABSORPTION ABSORPTION
IMPROPER PROCESSING (STACK) IMPROPER PROCESSING (STACK)
*
TYPES OF FILTERS
Instrument Compensation
Bandpass Filters
Deconvolution
Waveshaping
Spatial (Mixing)
K-F Filters (Velocity, Wedge, Pie Slice)
NARROWING THE PASSBAND NARROWING THE PASSBAND
CONSTANT FREQUENCY BANDWIDTH CONSTANT FREQUENCY BANDWIDTH
VARYING ENDPOINTS VARYING ENDPOINTS
STEEPENING THE SLOPES STEEPENING THE SLOPES
LENGTHENING THE OPERATOR LENGTHENING THE OPERATOR
RULE OF THUMB RULE OF THUMB
FOR FILTER SPECIFICATION... FOR FILTER SPECIFICATION...
FOURIER TRANSFORM FOURIER TRANSFORM
FILTER SCAN WITH NORMALIZING AFTER FILTER FILTER SCAN WITH NORMALIZING AFTER FILTER
No Filter 4-8-16-24 6-12-24-36 8-16-32-48 10-20-40-60 14-28-56-84 18-16-72-108
74
75
F-K Filter F-K Filter
76
F-K Filter F-K Filter
X
X
T T
+K
+K -K
-K
F F
REJECT or PASS
77
Contoh Aplikasi F-K Filter pada data VSP Contoh Aplikasi F-K Filter pada data VSP
78
Filter TEST? Filter TEST?
? ?
? ?
? ?
79
Deconvolution
80
Convolution is a mathematical operation defining the change of
shape of a waveform resulting from its passage through a filter.
The asterix denotes the convolution operator. In seismic, we
obtain a response for a certain model by convolving the seismic
signal of the source with the reflectivity function.
Convolution
81
Mathematically the convolution is defined as follows:
where k = 0 ... m+n; g
i
= (i=0 ... m) and f
j
= (j= 0 ... n).
The convolution can also be performed in the Fourier
domain:
Convolution in time domain = Multiplication (of the
Amplitude and spectrum and addition of the Phase
spectrum) in Fourier domain.
Convolution
82
Example of a Convolution
83
Example of a Convolution
84
Convolution of the reflectivity function with the signal of the source
returns the seismic trace.
Convolutional Model of a Seismogram
85
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:
i
i
i xy
y x . ) (
+
=
t
t |
where xi: (i=0 ... n); y
i
: (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.
Cross-Correlation
86
Cross-Correlation Function
87
The Auto-correlation is a Cross-correlation of a function
with itself. It is mathematically defined as:
i
i
i xy
x x . ) (
+
=
t
t |
where x
i
= (i=0 ... n);
xx
() = (-m < < +m)
and
m = max. displacement.
Auto-Correlation
88
Auto-Correlation of Two Identical Waveforms
89
Auto-Correlation: Multiples
90
A gradually decaying function indicative of short-
period reverberation.
A function with separate side lobes indicative of long-
period reverberations: multiples.
Autocorrelation Functions Contain
Reverberations
91
Auto-correlation :
Cross-correlation :
Normalization of Correlation
92
Convolution was the forward operation :
Source (t) * reflectivity (t) = signal (t)
Deconvolution is the reverse operation:
Reflectivity (t) = signal (t) *
-1
source (t)
Deconvolution
95
CONVOLUTIONAL MODEL CONVOLUTIONAL MODEL
WAVELET CONTRIBUTORS WAVELET CONTRIBUTORS
97
WAVELET CONTRIBUTORS WAVELET CONTRIBUTORS
98
IDEAL DECONVOLUTION IDEAL DECONVOLUTION
99
PULSE COMPRESSION PULSE COMPRESSION
100
SUMMARY OF PURPOSES SUMMARY OF PURPOSES
101
(

=
(

) 1 (
) 0 (
) 0 ( ) 1 (
) 1 ( ) 0 (
1 id
id o
ii ii
ii ii
a
a
|
|
| |
| |
Autocorrelation
Inverse filter
Cross-Correlation
i = input
d = desired
Optimum Wiener Filters for The
Deconvolution of Seismic Data
102
Objectives:
Compress the wavelet
into a sharp minimum- or
zero-phase shape.
Broaden (flatten) its
spectrum.
Deconvolution
103
Spiking and Predictive Deconvolution
Basic Idea
SUPPRESSION OF COHERENT NOISE
AFTER PREDICTIVE DECONVOLUTION
SUPPRESSION OF COHERENT NOISE
AFTER PREDICTIVE DECONVOLUTION
INPUT OUTPUT INPUT OUTPUT
105 106
A CMP stack (a) with no deconvolution, (b) with spiking deconvolution before
stack, (c) with signature processing (minimum-phase conversion of the
measured signature) followed by spiking deconvolution, (d) with signature
processing only (conversion of the signature to a spike).
107 108
109 110
111
Sorting from shot gathers to common-midpoint
(CMP) gathers.
In a CMP gather, the reflections all come from the
same point for flat layer.
We can also make common-receiver gathers.
A common-depth-point (CDP) is similar to a CMP, but
adjusts for the dip of the layers.
CMP Sort
112
Seismic data acquisition with multifold coverage is done in shot-receiver (s,g)
coordinates.
The figure is a schematic depiction of the recording geometry and ray path
associated with flat reflector. Seismic Data processing, on the other hand,
conventionally is done in midpoint-offset (y,h) coordinates. The required
coordinate transformation is achieved by sorting the data into CMP gathers.
Based on the field geometry information, each individual trace is assigned to
the midpoint between the shot and receiver location associated with that
trace. Those trace with the same midpoint location are grouped together,
making up a CMP gather.
Albeit incorrectly, the term common depth point (CDP) and Common
Midpoint (CMP) often are used interchangeably.
CMP SORTING
113
Figure depicts the geometry of a CMP gather and ray paths associated with
a flat reflector. Note that CDP gather is equivalent to a CMP gather only
when reflectors are horizontal and velocities do not vary horizontally.
However, when there are dipping reflectors in the subsurface, these two
gathers are not equivalent and only the term CMP gather should be used.
Selected CMP gathers obtained from sorting the deconvolved shot gathers.
Common Midpoint Sorting
114
Raypaths associated with
Common Shot Gather
Raypaths associated with
Common Mid-Point Gather
Raypaths associated with
Common Offset Gather
Raypaths associated with
Common Receiver position
Gather
Common Midpoint Sorting
115
CMP Gather
116
Shot Gathers and Midpoints
117
CMP Gathers Now
118
For most recording geometries, the fold of coverage nf
for CMP stacking is given by :
s
n
n
g g
f
A
A
=
2
g : receiver group
s : shot interval
n
g
: The number of recording channels
Fold of Coverage for CMP Stacking
NMO Correction of a CMP gather NMO Correction of a CMP gather
CMP Stack a CMP gather CMP Stack a CMP gather
Stacking
Example of stacking to
attenuate multiple
Example of stacking to
attenuate multiple
SIGNAL-TO-NOISE
AMPLITUDE IMPROVEMENT
SIGNAL-TO-NOISE
AMPLITUDE IMPROVEMENT
where N = number of unique
range traces stacked
where N = number of unique
range traces stacked
STACK STACK
Stack section Stack section
124
In addition to providing an improved signal-to-noise ratio,
multifold coverage with nonzero-offset recording yields
velocity information about the subsurface. Velocity
analysis is performed on selected CMP gathers or groups
of gathers. The output from one type of velocity analysis
is a table of numbers as function of velocity versus two-
way zero-offset time (velocity spectrum). These numbers
represent some measure of signal coherency along the
hyperbolic trajectories governed by velocity, offset, and
travel time.
Velocity Analysis
125
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:
V
stack
: The velocity that returns the best stacking result.
V
rms
: 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.
126
There are different ways to determine the velocity:
(t
2
-x
2
)-Analysis.
Constant velocity panels (CVP).
Constant velocity stacks (CVS).
Analysis of velocity spectra.
For all methods, selected CMP gathers are used.
Determine The Velocity
127
The (t
2
-x
2
) -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 v
2
, the square root returns the proper
velocity.
Example of a t
2
-x
2
-Analysis :
(t2-x2)-Analysis
128
Examining the normal moveout equation, it is possible to analyze
NMO velocities by plotting reflections in T
2
X
2
space.
129
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.
CVP - Constant Velocity Panels
130
Options in The ProMax Velocity Analysis
Routine
131
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.
CVS - Constant Velocity Stacks
132
Stacking Velocity
133
Concept of Constant Velocity Stack as An Aid
to Stacking Velocity Estimation
134
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.
Velocity-Spectrum
135
Amplitude of stacking
Normalized amplitude of stacking
Semblance
Normalized Amplitude of stacking :
Semblance :
Amplitude of Stacking :

=
=
n
i
t i t
w s
1
,
|
| |
,
1
t i
n
i
t
t
w
s
ns

=
=

=
t i
t i
t
t
t
w
s
n
Semblance
2
,
2
1
W
i,t
value for i-th trace, time t
Methods to Calculate The Velocity-Spectrum
136
Mapping of the offset axis to the velocity axis
Velocity Spectrum
137
Demonstration of The Velocity Spectra
138
Constant Velocity Stack (CVS)
One Method to Determine Stacking Velocity is to Use A Constant
Velocity Stack (CVS) for Several CDP Gathers.
139
Same CVS Panel of Traces as Before Switching to Variable
Density Color for The Traces to Utilize Dynamic Range.
Constant Velocity Stack (CVS)
140
Same as Previous Color Panels with Velocity Range Now
Halved to Better Pick Correct Velocities.
Constant Velocity Stack (CVS)
141
Constant Velocity Moveout (CVP) Corrections
(1)
142
Velocity ft/s
Constant Velocity Moveout (CVP) Corrections
(2)
143
Influence of Missing Long Offset Traces on
Velocity Spectra
144
Influence of Missing Long Offset Traces and
Statics on Velocity Spectra
145
Semblance-Analysis
146
147 148
Error for High Velocities and Large Travel
Times
Static
Corrections applied to seismic data to compensate
for the effects of variations in:
- elevation,
- weathering thickness,
- weathering velocity
(near-surface velocity anomalies) or
- reference to a datum.
1. What are residual static? 1. What are residual static?
Residual Static
CMP GATHER CMP GATHER
2. Why are residual static needed?
2. Why are residual static needed?
CMP GATHERS BEFORE AFTER RESIDUAL STATICS
STAT 5
After residual statics
STACK before & after residual statics STACK before & after residual statics
Before residual statics
153
STACK before & after residual statics STACK before & after residual statics
154
THE END OF
SEISMIC DATA PROCESSING
PART I

You might also like