3D Seismic Survey Design PDF
3D Seismic Survey Design PDF
3D Seismic Survey Design PDF
Theres more to designing a seismic survey than just choosing sources and receivers and shooting away. To
get the best signal at the lowest cost, geophysicists are tapping an arsenal of technology from integration of
borehole data to survey simulation in 3D.
C. Peter Ashton
Mrsk Olie og Gas AS
Copenhagen, Denmark
Christian Dplant
Elf Aquitaine
Pau, France
nCost of marine 3D
40
35
Dollars, in thousands
Brad Bacon
Angus Mann
Nick Moldoveanu
Houston, Texas, USA
30
25
20
15
10
5
DickiIreson
Thor Sinclair
Gatwick, England
Glen Redekop
Maersk Oil Qatar AS
Doha, Qatar
For help in preparation of this article, thanks to Jack
Caldwell and Greg Leriger, Geco-Prakla, Houston, USA;
Mandy Coxon and Dominique Pajot, Geco-Prakla,
Gatwick, England; Jacques Estival, Elf Petroleum Nigeria,
Lagos, Nigeria; Dietmar Kluge, Geco-Prakla, Hannover,
Germany; Lloyd Peardon, Schlumberger Cambridge
Research, England; Lars Sonneland, Geco-Prakla,
Stavanger, Norway; and Tim Spencer, British Gas,
Reading, England.
Appreciation is expressed to Qatar General Petroleum
Corporation (QGPC) for its consent to the release of data.
QUAD-QUAD is a mark of Geco-Prakla. TWST
(Through-Tubing Well Seismic Tool) is a mark of Schlumberger.
1. For the most recent worldwide figures:
Riley DC: Special Report Geophysical Activity in
1991, The Leading Edge 12, no. 11 (November
1993): 1094-1117.
2. Personal communication: Thor Sinclair.
April 1994
0
1990
1991
1992
1993
Year
19
Temporal Aliasing
50 Hz
nTemporal and
200 Hz
16
24
spatial aliasing
caused by sampling less than
twice per cycle.
Temporal aliasing
(top) occurs when
insufficient sampling renders a 50Hz signal and a
200-Hz signal indistinguishable
(arrows represent
sample points). The
50-Hz signal is adequately sampled,
but not the 200-Hz.
32
Time, msec
Extreme Aliasing
Two-way time
Minor Aliasing
Increasing offset
Increasing offset
aliasing (bottom)
occurs when
receiver spacing is
more than half the
spatial wavelength.
With minor aliasing
(left) arrivals can be
tracked at near offsets as time
increases, but
become difficult to
follow at far offsets.
With extreme aliasing (right) arrivals
even appear to be
traveling backwards, toward near
offsets as time
increases. (Adapted
from Claerbout, reference 6.)
Two-way time
Offset
Offset
Hyperbolic
moveout
+ Stacking
velocity =
CMP gather
Corrected
CMP gather
=
Stacked
CMP trace
Reflection arrival times from different offsets are assumed to follow a hyperbola. The shape of the hyperbola is computed from
the arrivals. Traces are aligned by flattening the best-fitting
hyperbola into a straight line, then summed, or stacked. Perfect
alignment should yield maximum signal amplitude at the time
corresponding to zero offset. A wide range of evenly spaced offsets gives a better-fitting hyperbola, and so a better stack.
20
What makes a good seismic signal? Processing specialists list three vital requirementsgood signal-to-noise ratio (S/N),
high resolving power and adequate spatial
coverage of the target. These basic elements,
along with some geophysical guidelines (see
Guidelines from Geophysics, page 22 ),
form the foundation of survey design.
High S/N means the seismic trace has
high amplitudes at times that correspond to
reflections, and little or no amplitude at
other times. During acquisition, high S/N is
achieved by maximizing signal with a seismic source of sufficient power and directivity, and by minimizing noise.3 Noise can
either be generated by the sourceshotgenerated or coherent noise, sometimes
orders of magnitude stronger than deep seismic reflectionsor be random. Limitations
in the dynamic range of acquisition equipment require that shot-generated noise be
minimized with proper source and receiver
geometry. Proper geometry avoids spatial
aliasing of the signal, attenuates noise and
obtains signals that can benefit from subsequent processing. Aliasing is the ambiguity
that arises when a signal is sampled less
than twice per cycle (left ). Noise and signal
cannot be distinguished when their sampling is aliased.
A common type of coherent noise that
can be aliased comes from low-frequency
waves trapped near the surface, called surface waves. On land, these are known as
ground roll, and create major problems for
processors. They pass the receivers at a
much slower velocity than the signal, and
so need closer receiver spacing to be properly sampled. Planners always try to design
surveys so that surface waves do not contaminate the signal. But if this is not possible, the surface waves must be adequately
sampled spatially so they can be removed.
During processing, S/N is enhanced
through filters that suppress noise. Coherent
noise is reduced by removing temporal and
spatial frequencies different from those of
the desired signal, if known. Both coherent
and random noise are suppressed by stackingsumming traces from a set of sourcereceiver pairs associated with reflections at
a common midpoint, or CMP.4 The sourcereceiver spacing is called offset. To be
stacked, every CMP set needs a wide and
evenly sampled range of offsets to define the
reflection travel-time curve, known as the
normal moveout curve. Flattening that
curve, called normal moveout correction,
will make reflections from different offsets
arrive at the time of the zero-offset reflection. They are then summed to produce a
stack trace (left ). In 3D surveys, with the
April 1994
Source
Bin
Receiver
Offset Distribution
1
160
190
220
130
40
36
160
33
29
190
25
22
220
Fold
18
14
11
250
7
3
nA fold plot showing 40-fold coverage over the heart of the survey.
The edge of the survey has partial fold because several of the first
and last shots do not reach as many receivers as in the central part
of the survey.
21
a
0
Seafloor reflection
1.0
Time, sec
Seafloor multiple
2.0
Seafloor multiple
Primary reflection
3.0
Multiple
Multiple
4.0
Primary
reflection
Ghost
Near-surface
multiples
Long-path
multiple
nMultiple reflec-
or between the earths surface and the bottom of a layer of unconsolidated rock
(below, left ). Multiples can appear as later
arrivals on a seismic section, and are easy to
confuse with deep reflections (left ).6 And
because they can have the same characteristics as the desired signalsame frequency
content and similar velocitiesthey are
often difficult to suppress through filtering
and stacking. Sometimes they can be
removed through other processing techniques, called demultiple processing, but
researchers continue to look for better ways
to treat multiples.
The second characteristic of a good seismic signal is high resolution, or resolving
powerthe ability to detect reflectors and
quantify the strength of the reflection. This is
achieved by recording a high bandwidth, or
wide range of frequencies. The greater the
bandwidth, the greater the resolving power
of the seismic wave. A common objective of
seismic surveys is to distinguish the top and
bottom of the target. The target thickness
determines the minimum wavelength
required in the survey, generally considered
to be four times the thickness.7 That wavelength is used to calculate the maximum
required frequency in the bandwidth
average seismic velocity to the target
divided by minimum wavelength equals
maximum frequency. The minimum frequency is related to the depth of the target.
Lower frequencies can travel deeper. Some
seismic sources are designed to emit energy
in particular frequency bands, and receivers
normally operate over a wider band. Ideally,
sources that operate in the optimum frequency band are selected during survey
design. More often, however, surveys are
shot with whatever equipment is proposed
by the lowest bidder.
x =
V rms
4f max sin
ple 3D surveys:
y = x.
X max =
2Tv 2
f v
v
2D fold
22
xy
2Rf dx
1. Normal moveout stretch is the distortion in waveshape caused by normal moveout correction.
a
Horizontal Reflector
April 1994
Dipping Reflector
Shotpoint
Receiver
flat (top) the CMP is a common reflection point. When the reflector
dips (bottom) there is no CMP. A dipping reflector may require
changes in survey parameters, because reflections may involve
more distant sources and receivers than reflections from a flat layer.
23
Theoretical Grid
Final Grid
Checkerboard Pattern
Brick Pattern
Zigzag Pattern
Source
Receiver
Source
Receiver
USA (top, left) calls for a checkerboard of receiver lines (blue) and
source lines (red). The actual survey shot (bottom, left) came very
close to plan. Other cities present acquisition challenges. A survey
in Milan,iItaly (right) used a random arrangement of sources and
receivers. (Adapted from Bertelli et al, reference 9.)
24
Oilfield Review
nMarine acquisi-
tion geometry
showing seismic
vessels looping in
oblong circuits.
The length of
straight segments is
calculated from
fold plots, and must
include additional
lengthrun in
and run outto
allow cable to
straighten after
each turn.
tens of thousands of dollars per day to operate. Sources are clusters of air guns of different volumes and receivers are hydrophones
strung 0.5 m [1.6 ft] apart in groups of up to
48, on cables up to 6000 m [19,680 ft]
long. Sources and receivers are almost
always towed in straight lines across the target (below, right ), although other geometries
are possible. Circular surveys have been
acquired with sources and receivers towed
by vessels running in spirals or concentric
circles.10 Geco-Praklas QUAD-QUAD system tows four receiver cables and four
source arrays simultaneously, acquiring 16
lines at a time. Currents and tides can cause
the long receiver cables to deviate by calculable amountsup to 30from the towing
direction. Spacing between shotpoints is a
function of vessel speed, and can be limited
by how quickly the air guns can recover full
pressure and fire again. Access is usually
limited only by water depth, but drilling
rigs, production platforms and shipping
lanes can present navigational obstacles.
Environmental constraints also influence
marine surveys: the commercial fishing
industry is imposing limits on location of,
and seasons for, marine acquisition.11 For
example, planning in the Caspian Sea must
avoid the sturgeon breeding season or seismic surveys would wipe out caviar production for the year.
Transition zonesshallow water areas
have their own problems, and require specialized equipment and creative planning.12
Transition zones are complex, involving
shorelines, river mouths, coral reefs and
swamps. They present a sensitive environment and are influenced by ship traffic,
commercial fishing and bottom obstructions. Survey planners have to contend with
varying water depths, high environmental
noise, complex geology, wind, surf and
multiple receiver typesoften a combination of hydrophones and geophones.
One thing all surveys have in common is
that planning must be done quickly. The
9. Bertelli L, Mascarin B and Salvador L: Planning and
Field Techniques for 3D Land Acquisition in Highly
Tilled and Populated AreasTodays Results and
Future Trends, First Break 11, no. 1 (January 1993):
23-32.
10. Hird GA, Karwatowski J, Jenkerson MR and Eyres A:
3D Concentric Circle SurveyThe Art of Going in
Circles, EAEG 55th Meeting and Technical Exhibition, Stavanger, Norway, June 7-11, 1993.
11. Gausland I: Impact of Offshore Seismic on Marine
Life, EAEG 55th Meeting and Technical Exhibition,
Stavanger, Norway, June 7-11, 1993.
12. Petersen C, Brakensiek H and Papaterpos M:
Mixed-Terrain 3D Seismics in the Netherlands,
Oilfield Review 4, no. 3 (July 1992): 33-44.
April 1994
Bin
25
Specification
Define survey
objectives
Required
geophysical
parameters
Resolution
analysis
Design
No Required Yes
equal
obtainable?
Source,
template
and array
design
Planning
Preferred
survey
parameters
No
Objectives Yes Final survey
design
achieved?
Evaluation
Analysis of
existing data
Data
Type
Resolution, noise
and coverage
analysis
Parameters
to be Determined
Maximum frequencies
attainable
Prospect
description
Obtainable
geophysical
parameters
Means to
Determine Parameters
VSP processing
Reflection response
of target
VSP
Identification
of multiples origin
Source peak amplitude
Source modeling
Source volume
Resolution attainable
Process or Output
Loss modeling
Frequency
dependent losses
Peak-to-bubble ratio
Source depth
Operational,
cost and safety
constraints
Source signature
for various depths
Bandwidth at target
Apply losses to
source signatures
Target wavelet
Noise levels
Logs or
1D Models
Synthetic
CMP gathers
Synthetic shots
Migration aperture
Long-offset analysis
Normal incidence
stacks
Statics model
Signal-to-noise ratio
Establish noise
mechanisms
Refraction velocities
(near surface)
Modeled section
Build geological
2D model and
apply appropriate
target wavelet
Analysis of existing
surface seismic
Amplitude versus
time plots
Group interval
Mute, stack,
fold tests
Crossline spacing
Spatial frequency
Spatial resolution
Shotpoint interval
Noise records
Analysis of
migration
requirements
Migration aperture
Shooting direction
Record length
Migration of
synthetic
zero-offset data
Migration of
existing 2D data
Ambient noise
estimation
26
Oilfield Review
Shot Depth 28 m
April 1994
1.0
Time, sec
Shot Depth 9 m
2.0
3.0
4.0
Ground roll
Offset
300 m
Source
Receiver
27
0-10 Hz
10-20 Hz
20-30 Hz
30-40 Hz
nFiltered 2D data
showing frequency
content variation
with depth. Each
panel has been filtered to allow a different band of frequencies, called the
passband, to pass.
As the passband
rises, the maximum
depth of penetration
of seismic energy
decreases. Lower
frequencies (left)
penetrate deeper.
Higher frequencies
(right) do not propagate to deeper levels. At the target
level of 3.0 sec there
is still some 50 Hz
energy left.
1.0
Time, sec
40-50 Hz
2.0
3.0
4.0
28
Oilfield Review
5-10 Hz
10-20 Hz
20-40 Hz
40-60 Hz
60-80 Hz
80-100 Hz
Time, sec
0.5
1.0
1.5
April 1994
2400
90
3325
80
4250
70
60
5175
Level number
Measured depth, ft
6100
50
40
7025
30
7950
20
8875
10
9800
0
0
0.5
Time, sec
1.0
80
70
60
50
nVertical seismic profile (VSP) traces (left) analyzed for amplitude loss with depth (right).
Amplitudes of first arrivals recorded in a 92-level VSP are calibrated with amplitudes of a
surface reference signal to account for changes in source amplitude from level to level.
The amplitude ratio from one level to the next is plotted in decibels (dB). One dB is 20
times the log of the amplitude ratio. An amplitude ratio of 100 is equivalent to 40 dB.
Amplitudes expected from a surface seismic survey would normally be 3 dB less than
those from a VSP, and scaled by a reflection coefficient.
29
3-m throw
6-m throw
12-m throw
24-m throw
0.3
0.4
0.5
12-m
thick
0.6
24-m
thick
0.7
36-m
thick
0.8
nResolution of thin beds and small-scale faulting. Each panel shows the modeled
response of a seismic wave of 48-m [160-ft] wavelength ( ) to a different vertical fault displacing a series of thin beds of thicknesses 12 m, 24 m and 36 m. From left to right, faults
with 3-m [10-ft], 6-m [20-ft], 12-m [40-ft] and 24-m [80-ft] throws correspond to /16, /8,
/4 and /2, respectively. A fault throw of at least 12 m, corresponding to /4, can be
resolved quantitatively. At less than that, existence of a fault can be detected, but its
throw resolved only qualitatively.
-128
Amplitude
30
127
miles
km
3.2
Oilfield Review
0
1000
2000
Depth, m
Salt
3000
4000
5000
6000
0
4000
8000
12,000
16,000
20,000
24,000
28,000
Distance, m
nRaytrace modeling showing strong changes in reflection paths through salt. Traces that
would have a common midpoint in a flat-layered earth no longer bounce in the same
bin. Salt, with its ability to deform and its high seismic velocity, creates complex structure
and strong refraction, or ray bending.
TEXAS
LOUISIANA
Ship Shoal
South Addition
G U L F
0
0
miles
km
O F
M E X
I C O
100
161
Some of the advances to be made in 3D survey design have origins in other fields. VSP
design routinely models seismic raypaths
through complex subsurface structure, but
rarely does surface seismic design account
for structure. Despite considerable sophistication in 3D data processing, most 3D survey design assumes plane layer geometry in
the subsurface to calculate midpoints and
target coverage. But to estimate subsurface
coverage adequately in complicated struc15. Migration, sometimes called imaging, is a processing
step that rearranges recorded seismic energy back to
the position from which it was reflected, producing
an image of the reflector.
16. Poster C: Taking the Pulse of 3D Seismics, Middle
East Well Evaluation Review, no. 13 (1992): 6-9.
April 1994
31
0
1000
Depth, m
2000
3000
4000
5000
6000
0
4000
8000
12,000
16,000
20,000
24,000
28,000
20,000
24,000
28,000
Distance, m
0
1000
Depth, m
2000
3000
4000
5000
6000
0
4000
8000
12,000
16,000
Distance, m
nRaytrace modeling to optimize cable length. Refraction through salt may mean a longer
cable is required to image structure below. Two cable lengths, 8075 m (top) and 5425 m
(bottom) were tested using the model on the previous page. Surprisingly, in this case both
cables give similar coverage of subsalt horizons.
32
Oilfield Review