Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

How Karcher "hand migrated" dipping layer Top: Monument unveiled in 1971 at Belle Isle (Oklahoma City) on 50th anniversary of first seismic reflection survey by J. C. Karcher.

1 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Middle: Two early reflection records from Belle Isle, 1921. Bottom: Karcher's interpretation of same. uses reflected energy from interfaces between subsurface layers to determine their configuration reflections recorded as two-way (down and back up) travel times, not depths fraction of incident energy reflected from interface called reflection coefficient dependent on acoustic impedance contrast across interface

Polarity of reflected wave depends on sign of reflection coefficient (unchanged polarity means compression remains compression, dilatation remains dilatation)

Hypothetical Rock Properties Rock Granite Basalt Limestone Sandstone Shale VP, km/s 5.0 5.5 6.0 4.2 2.5 , kg/m 2700 3000 2300 2500 2300
3

Vx 13,500 16,500 13,800 10,500 5,750

For these hypothetical values, limestone-granite contact will be poor reflector

Simple Zero-offset Reflection Survey

zero offset (distance between source and receiver) single layer on half-space

reflections produce a time-section approximating interface (that even geologists can understand) travel time:

2 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

2 problems: what is V? [We measure t. If we knew V, we could find d - convert time section to depth section.] energy source is expensive in time and money solution to both problems: use geophones at different offsets

Seismic Reflection Survey with Offset, Single-fold Coverage

record traces from several geophones spaced away from source (shot) display traces side-by-side, so distance between traces proportional to geophone spacing display increasing time downward (time approx. proportional depth)

Note that subsurface reflection points have half the spacing of geophones. To get complete "single-fold "coverage of the subsurface, can shoot from either end of geophone spread:

One can also use a "split-spread" arrangement, here with shot at point B, then move half the geophones forward and shoot at C:

The next two figures show recording-truck signal check for 36-channel split-spread layout:

3 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

4 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

The travel time for the primary reflection (first layer) where geophone offset = x, thickness d, velocity V

5 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

this is a hyperbola:

for multiple geophones, seismic traces look like:

direct ray is straight line, per rate equation reflections from 1st and 2nd reflectors are not flat; reflections are hyperbolas because of normal moveout (NMO): reflection time increases with x, nonlinearly normal move-out (NMO): the difference in reflection travel-times from a horizontal reflecting surface due to variations in the source-geophone distance can correct for this using NMO correction, so reflections are flat because deeper reflectors produce flatter parabolas, NMO is less:

6 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

given that :

the NMO is just t - t0, where t0 is simply 2d/V (zero-offset 2-way travel time) but we don't know d or V however,

so, plot t2 vs. x2:

measure slope to get V intercept gives d

Multiple Layers
foregoing assumed a straight-line path from source to reflector to receiver with multiple layers with different velocities, this clearly does not hold (actual path compared with straight-line assumption):

7 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

however, if depths are large compared to total geophone spread, error can be small Green Method assuming straight-line paths, one can still just use a x2-t2 plot to estimate velocities and depths however, there is a more accurate method: The Dix Equation uses special velocity called VRMS

still assumes nearly vertical incidence/straight line raypaths

given n horizontal beds, and is the one-way vertical travel-time through bed i, Dix equation states:

use x2-t2 plot to determine RMS velocities to each layer then can get interval layer velocities and thicknesses replace velocity terms by VRMS

it can be shown that Dix's equation can be solved for the individual interval velocities in fact, the following is sometimes referred to as Dix's equation:

thicknesses can then be easily determined:

Velocity Scans - another way to get stacking velocity:

8 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Signal Summing; Stacking

As seismic energy moves away from the source, there is a decrease in signal strength as the energy spreads out. This causes energy to decrease by E = E0/(2r2). Since energy is proportional to the square of the amplitude, the signal amplitude drops off like 1/r. In addition, since rocks are not truly elastic (anelastic), some energy is lost to heat with every cycle, leading to an exponential loss of energy. Higher frequencies go though more cycles to a given depth (shorter wavelength), so high frequency energy is lost with depth Taking both of these effects into consideration, we have:

reflections from significant depth have amplitudes that may be well below the noise level summing and stacking adds (coherent) signals and (random) noise improves signal-to-noise ratio (S/N, or S/(S+N)) summing n times increases signal by n summing n times increases noise by square root of n Example: reshoot a line 36 times signal increases by 36

9 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

noise increases by 6 (square root) S/N improves by 36/6 = 6 different means of improving S/N: geophone groups "Geophones are rarely used singly. Normally several (as many as 20 or more) are electrically connected to each other in a group in such a way that the outputs of the individual phones are effectively summed. The information from each group must be transmitted via cables to the recording truck. In modern land recording with 48, 96, or more group recordings, the cables are long and heavy and often add noise to the recording, especially in the presence of powerlines or water." - Dobrin and Savit, Introduction to Geophysical Prospecting, 4th ed. "geophone" is actually a group of geophones "tied together" in a geophone group signals in parallel, fed into one channel of system signal usually has small incident angle, reaches all geophones together (coherent) surface noise sweeps across geophones, tends to cancel multiple shots dynamite in hole vibroseis trucks: multiple trucks, all in sync; shake several times multiple hammer blows, shotgun blasts, etc. multi-fold coverage Example: 4 geophones (channels); move shot and geophones one geophone spacing and reshoot:

note that subsurface reflection points twice as closely spaced as geophones had we moved shot (array) 2 geophone spacings, only get single-fold coverage

try this with different numbers of geophones and shot spacings to find a simple formula to calculate fold-coverage (coverage = f(# geophones, # geophone spacings moved between shots) Rock Velocities
P, S velocities of various rocks (See Sydney Clark's GSA Memoir for much more) velocity vs. density velocity vs. density, part 2 velocity vs. rock age for sandstones and shales

Data Collection
Source

10 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Geophones

11 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Geophones (~ $100 each) have a typical natural frequency of 10 Hz Response is relative good over a range from about 2 Hz to 100 Hz

12 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Recording Digitally Analog signal from phones (continuous voltage vs. time) sampled (typically every 2 msec) by an A-to-D (A/D) converter originally integer recording example: 16-bit recording; 216 = 65,536, or -32767 to +32768 dynamic range of about 4.5 orders of magnitude poor resolution at low amplitudes floating pointing recording single precision, 4 bytes, 6 digits of resolution, dynamic range 10+/-32

Processing Steps
At one time, greatest computing power was owned by government (mostly DOD) Petroleum companies ranked second most seismic reflection processing is computer intensive, but requires intelligent "operator" input at many steps in the process AGC: automatic gain control early-arriving reflections may be orders of magnitude larger in amplitude than later ones AGC looks at average amplitude in a sliding time window and boosts (or attenuates) amplitude to a constant value over that window AGC causes loss of true amplitude information modern floating-point recording allows full amplitude information to be retained retaining "relative, true amplitude" done with linear or quadratic increase of gain with time:

13 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Filtering remove or attenuate certain frequencies to reduce noise and improve S/N notch filter common at 60 Hz filtering often done with FFT filters must be "ramped" to avoid ringing (Gibbs phemonemon)

14 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Statics Removal

refraction statics - requires reversed profile up-hole shooting vertical velocity distribution near the surface determined by shooting up the hole" (Geophysical Services, Inc.):

15 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

automatic statics

16 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Migration

17 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

18 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Synthetic Seismograms

19 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Directions Seismic Reflection is Heading

20 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Percentage of seismic activity involving various techniques. (Data from SEG annual Geophysical Activity Reports, pre-1981 data are for U.S. activity, post-1980 for worldwide activity, 3-D data from Dutt, 1992, adjusted according to judgment expressed in Goodfellow, 1991.)

Seismic Attributes
reflections are not only information available in seismic data already seen value of preserving "relative, true" amplitude preservation and display of velocity data can reveal info otherwise missed: Conventional B&W section on which carbonate bank would be missed:

21 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

Color display in which colors are keyed to interval velocity estimates (1000 ft/s increments):

Close-up of carbonate bank sequence seen above:

22 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

3-D Seismic Reflection

3-D representation ("data cube") migration out of plane of section (side-swipe) geology is, after, 3-D much more expensive! (~n2)

23 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

3-D sesimic time slices at time ranging from 1060 ms to 1260 ms:

24 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

3D visualization (caves, virtual reality, etc.) wavelet processing

Summary of Processing Steps

25 de 26

01/02/2012 15:20

Seismic Reflection

http://principles.ou.edu/seismic_explo/reflect/reflect.html

26 de 26

01/02/2012 15:20