Sesi - 07A - Seismic Processing - Basic
Sesi - 07A - Seismic Processing - Basic
by
TAUFIK ALIF MANAN
Geophysics Physics
University of Indonesia
SEISMIC
DATA
PROCESSING
An Overview
Processing of Reflection Data
• Improve the signal to noise ratio: e.g. by measuring of several
channels and stacking of the data (white noise is suppressed).
• 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);
Flow Chart of 2D Land Seismic Reprocessing
Data Input
Field Report
Field
-Observer’s report
Tape
-Geometry information
Reformat
Build geometry database
MERGE reformatted field tape
With geometry database
Preprocessing
- First Break Mute
- Editing
- Low / high cut filter
- Amplitude correction (Spherical Divergence)
- Deconvolution (Surface Consistent Predictive Decon)
- Refraction Static applied
Brute Stack
Sorting to CMP gather
using single velocity
1st Velocity Analysis
1st Residual Static Corr.
2nd Velocity Analysis
2nd Residual Static Corr.
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.).
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
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.
CMP Processing Sequence (continued)
10. Deconvolution or Wavelet Processing.
• Compresses the wavelet in time, attenuates reverberations.
• Converts the wavelet to zero-phase for viewing.
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.
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.
• Shot-Static
= (Elevation of source - Elevation of reference level) /
Velocity.
• Receiver-Static
= (Elevation of receiver - elevation of the reference
level) / Velocity.
• 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.
• Geometrical spreading.
• Absorption.
• Receiver response.
• Measurement system
Problem for Data Processing
• Trace equalization.
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)
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.
Filtering
SPLIT-SPREAD FIELD RECORD
TYPES OF NOISE
CAUSES OF POOR SIGNAL-TO-NOISE
NON-OPTIMAL FIELD PROCEDURES
SCATTERING OR ABSORPTION
- VULCANICS
- FRACTURED ZONE
NEAR-SURFACE PROBLEMS
- POOR COUPLING
- ABSORPTION
IMPROPER PROCESSING (STACK)
TYPES OF FILTERS
Instrument Compensation
Bandpass Filters
Deconvolution
Waveshaping
Spatial (Mixing)
X
-K +K
F
T
REJECT or PASS
Contoh Aplikasi F-K Filter pada data VSP
Filter TEST?
?
Deconvolution
Deconvolution
• In a CMP gather, the reflections all come from the same point for flat
layer.
* * *
Albeit incorrectly, the term common depth point (CDP) and Common
Midpoint (CMP) often are used interchangeably.
Common Midpoint Sorting
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.
CMP Gathers
Common Midpoint Sorting
ng g Δg : receiver group
nf Δs : shot interval
2s ng : The number of recording channels
Example of a t2-x2-Analysis :
Examining the normal moveout equation, it is possible to analyze
NMO velocities by plotting reflections in T2 X2 space.
CVP - Constant Velocity Panels
• 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.
Stacking Velocity
Concept of Constant Velocity Stack as An Aid
to Stacking Velocity Estimation
Velocity-Spectrum
• Amplitude of stacking
• Normalized amplitude of stacking
• Semblance
n
Amplitude of Stacking : st wi ,t
i 1
| st |
Normalized Amplitude of stacking : nst n
i 1
wi ,t |
1 s 2
t
Semblance : Semblancet t
n w 2
i ,t
t i
Velocity ft/s
Influence of Missing Long Offset Traces on
Velocity Spectra
Influence of Missing Long Offset Traces and
Statics on Velocity Spectra
Semblance-Analysis
Error for High Velocities and Large Travel
Times
Velocity analysis and NMO-Correction
• Introduction
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.
From this formula the NMO-correction can be derived and is given by:
2 x2
t t0 t ( x) with t ( x) t 2
0
The Moveout Δt is the difference in traveltime vstack for a receiver at a
distance x from the source and the traveltime t0 for zero-offset distance.
Change of Moveout with Depth
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.
2. 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 heavily than
the short-offset traces, because the difference in NMO between
primaries and multiples is larger for larger offsets. A weight factor α is
introduced.
1 n
st ,meanweighted i wi ,t
n i 1
Methods of Stacking (Continued)
3. Diversity stacking / Min-Max-exclude
This processing still assumes that the static shifts are caused by the
interface. Therefore, this processing is also called surface
consistent static correction.
Methods to Determine The Residual Statics
• Delay-Time.