10 Seismic Stratigraphy
10 Seismic Stratigraphy
10 Seismic Stratigraphy
Seismic Stratigraphy
Reflection seismology is compartmentalized into acquisition, processing and interpretation. Seismic stratigraphy deals with interpretation. It is the study of seismic data for the purpose of extracting stratigraphic information. Seismic stratigraphy is often divided into several sub-areas: Analysis of seismic sequence Separating out time-depositional units based on detecting unconformities or changes in seismic patterns; Analysis of seismic facies Determining depositional environment from seismic reflection characteristics; Analysis of reflection character Examining the lateral variation of individual reflection events, or series of events, to locate where stratigraphic changes occur and identify their nature; the primary tool for this is modeling by both synthetic seismograms and seismic logs.
1
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
The fold of stacking refers to the number of traces in the CDP gather and may conventionally be 6. 12, 24, 48.. CDP spacing = 1/2 x group (channel) interval maximum CDP fold = Velocity increases with depth: Time section is compressed in terms of depth at greater times. Hence depth converted profiles have better representations of geometry. Vertical Exaggeration: Usually significant, and variable in depth.
Boggs (2001), p.498
4
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Nature of a reflection
1 1 = + v vf vm
v: velocity in the saturated rock : porosity vf: fluid velocity vm: velocity of the matrix
Velocity-density relationship
= 0.23v 0.25
Gardners equation
Impedance contrasts: Reflections caused by impedance contrasts. The impedance (Z) of a rock unit is the product of its velocity (v) and density (), that is Z= v. At normal incidence, Reflection Coefficient: RC = (Z2-Z1)/(Z1+Z2)
5
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Amplitude
Reflection amplitude has to do with seismic wave height and is a function of the energy of seismic waves. On a seismic record, amplitude is measured as the distance from the mid-position of a wave to the extreme position. Amplitude is directly proportional to RC. It is also affected by the spacing between reflecting surfaces. Where bed spacing is optimum, lower energy responses are phased together constructively (constructive interference) to intensity or amplify the reflected energy and thus increase amplitude.
6
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Polarity
Positive RC produces a positive reflection, by definition and negative RC produces a negative reflection. Determine polarity from known impedance boundary, for example the water bottom (positive).
Where to Pick? Actual onset of reflection corresponds to impedance contrast or geological boundary in Minimum Phase Data. Data can be processed to Zero Phase such that peak amplitude of a symmetrical wavelet lies over impedance contrast. In any case you pick on the peak because that is what is easy and in the case of minimum phase data, make any necessary adjustment for the distance between the reflector and the geologic boundary. Most, if not all, seismic sections are displayed in minimum phase data. (Definitions: Minimum phase: a characteristic of waveforms which have their energy concentrated early in the waveform; Zero phase: a characteristic waveforms which are symmetrical.)
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Minimumphase wavelet
Zero-phase wavelet
Frequency: The frequency spectrum of the acoustic signal generated varies according to the energy sources.
8
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
The higher the frequency of the waveform, the greater or better will be the vertical resolution.
Doyle and Bennett (1998), p.285
Because velocity increases with depth and frequency decreases with depth resolution goes down. Typically tops and bottoms of beds resolved by TWT of /2. Once the bed is about /30 thick or less, reflections from the top and base effectively cancel and there is no detectable seismic response. 10
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Thin-bed effect or tuning Reflectors that are spaced more closely than onequarter thick (/4 ) of the wavelength have responses that begin to add constructively to produce a reflection with high amplitude.
/30 /4 /2
Badley (1985), p.18
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
11
Interference: the superposition of waveforms. Many (perhaps most) reflections are the interference composites resulting from several interfaces. So there is no one-to-one correspondence between seismic events and interfaces in the Earth.
Examples of interference on a zero-phase normal-polarity wavelet for a range of bed thickness and spacings.
12
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
B. Lateral resolution The seismic energy travels as wave fronts and the region on the reflector where the seismic energy is reflected constructively is known as the Fresnel zone. Lateral resolution is determined by the radius of the Fresnel Zone, which itself depends on the wavelength of the acoustic pulse and the depth of the reflector. Fresnel zone are typically hundreds of meters.
For example, for a plane that reflects interface at a depth of 1,000 m and an average velocity of 2,000 m/sec, the first Fresnel zone has a radius of 130 m for a 60-Hz component and 183 m for a 30-Hz component. For a reflection at 4,000 m with an average velocity of 3,500 m/sec, the first Fresnel zone has a radius of 375 m for a 50-Hz component and 594 m for a 20-Hz component.
Sheriff (1977) in AAPG Mem.26, p.11
13
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
hole
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
14
Vail Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
15
(3)
Another example showing the chronostratigraphic significance of seismic reflections (Sherrif and Geldart, 1995, p. 403)
16
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
next page
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
17
18
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
An idealized sequence
Vail (1987)
Seismic sequence analysis involves identification of major reflection packages that can be delineated by recognizing surfaces of discontinuity. Discontinuities may thus be recognized by interpreting systematic patterns of reflection terminations along the discontinuity surfaces.
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
19
toplap
Erosional truncation is the termination of strata against an overlying erosional surface. onlap
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
http://ic.ucsc.edu/~casey/ eart168/Lec.SeisStrat.htm
20
Apparent truncation is the termination of relatively low-angle seismic reflections beneath a dipping seismic surface, where that surface represents marine condensation.
Lapout is the lateral termination of a reflection (generally a bedding plane) at its depositional limit. Baselap is the lapout of reflections against an underlying seismic surface (which marks the base of the seismic package). Baselap can consist of onlap or downlap.
21
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Onlap is recognized on seismic data by the termination of low-angle reflections against a steeper seismic surface. Two types of onlap are recognized: marine onlap and coastal onlap. Downlap is baselap in which an initially inclined stratum terminates downdip against an initially horizontal or inclined surface. The surface of downlap represents a marine condensed unit in most cases. Santa Cruz terrace deposits downlapping onto unconformity.
http://ic.ucsc.edu/~casey/eart168/Lec.SeisStrat.htm
22
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Toplap is the termination of inclined reflections (clinoforms) against an overlying lower angle surface, where this is believed to represent the proximal depositional limit. Other term: Offlap: A conformable sequence of inclined strata, deposited during a marine regression, in which each stratum is succeeded laterally by progressively younger units (a clinoform). Clinoforms merging into toplap. Peru, a temperate water carbonate of Miocene age.
23
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
http://ic.ucsc.edu/~casey/eart168/Lec.SeisStrat.htm
24
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
25
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
26
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Examples of toplap
1. Determine the vertical and horizontal scale of the section. 2. Migrated section or not, marine or land data? 3. Identify multiples (e.g., water-bottom multiples, peg-leg multiples etc.) and mark them in light blue by convention. 4. Identify and mark reflection terminations or unconformities (e.g. onlap, downlap, truncation etc.) with arrowheads (by convention use a red pencil). 5. Identify seismic surfaces on the basis of reflection terminations (A seismic surface is a line on a seismic section where reflections terminate in a consistent manner.). In initial stages, mark seismic surfaces in yellow colour by convention. Then assign specific colour to individual seismic surface based on its type (e.g. sequence boundary, transgressive surface, maximum flooding surface etc.) or age in later stages of stratigraphic interpretation. 6. Identify sequence boundaries. Sequence boundaries are commonly marked by truncation or onlap, whereas maximum-flooding surfaces are commonly marked by downlap. 7. Carry out a similar exercise on other intersected seismic lines and tie the seismic surfaces and interpretation (i.e., ensure that the interpretation is consistent where lines cross) around the data set. 8. Mapping sequence units on the basis of thickness, geometry, orientation, or other features to see how each sequence relates to neighboring sequences. 9. Identify seismic facies for each sequence. 10. Interpretation of depositional environments and lithofacies.
28
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Practical: On the following figure, overlay a piece of tracing paper and do the following: 1. Pick reflection terminations; 2. Draw lines of seismic surfaces. The seismic data is from the Outer Moray Firth, central North Sea, showing the seismic stratigraphy of the post-Palaeocene section. The surface at around 0.7 second is dated as close to the top of the Miocene, and the high-relief surfaces in the shallower section are interpreted as glacial lowstand surfaces in the PliocenePleistocene.
29
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
sec
30
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Geological Interpretation (Sangree Widmier, 1979) Impedance (velocity-density) contrasts, Layer spacing (cause constructive and destructive interference), Fluid content Bed spacing, Fluid content Beddding or layer continuity, depositional processes
31
Seismic response for a sand with a gradational base, which results in lower amplitude. The 9-m thickness is about 1/8 wavelength.
32
2. Parallelism of reflection cycles to gross bedding, and therefore, to physical surfaces that separate older from younger sediments : Reflection configurations. Reflection configuration refers to the gross stratification patterns identified on seismic records.
Feature of Reflectors Reflection Configuration (pattern) External form and areal association of seismic facies units
Geological Interpretation (Sangree Widmier, 1979) Stratification patterns, Depositional processes, Erosion and paleotopography
33
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Principle reflection patterns 1. Parallel and subparallel: generated by strata that were probably deposited at uniform rates on a uniformly subsiding shelf or in a stable basin setting.
34
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
2. Divergent: Divergent configurations are characterized by a wedge-shaped unit in which lateral thickening of the entire unit is caused by thickening of individual reflection subunits within the main unit. Divergent configurations are interpreted to signify lateral variations in rates of deposition or progressive tilting of the sedimentary surface during deposition.
35
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
3. Prograding: Generated by strata that were deposited by lateral outbuilding or progradation to form gently sloping depositional surfaces called clinoforms. Prograding reflection configurations may include patterns of sigmoid (superposed Sshaped reflectors) and oblique, complex sigmoid-oblique, shingled, hummocky.
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
36
Sigmoid
Tangential oblique
Parallel oblique
Complex sigmoid-oblique
Mitchum et al. (1977), p.126
37
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
38
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
4. Chaotic: This pattern is interpreted to represent a disordered arrangement of reflection surfaces owing to penecontemporan eous, softsediment deformation, or possibly to deposition of strata in a variable, high-energy environment.
Mitchum et al. (1977), p.129
39
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
5. Reflection-free: This pattern may represent homogeneous, nonstratified units such as igneous masses or thick salt deposits, or highly contorted or very steeply dipping strata.
Modifying terms
40
Sediment waves
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
41
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
42
Interpretation of lithofacies and depositional environments Once the objective aspects of delineating seismic sequences and facies have been completed, the final objective is to interpret the facies in terms of lithofacies, depositional environments, and paleobathymetry. The most useful seismic parameters in seismic faces analysis are the following: 1. The geometry of reflections (reflection amplitude, continuity, frequency) and reflection terminations (onlap, downlap, erosional truncation, toplap). 1. Reflection configuration (parallel, divergent, sigmoid, or oblique) 2. Three dimensional form.
43
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
Seismic facies classification A,B,C technique for two-dimensional seismic facies analysis (Ramasayer, 1979) These codes can be marked on a map:
There is no unequivocal link between seismic facies and depositional systems, with the probable exception of the link between clinoforms and slope systems. Continuous flat-lying reflections may, for example, reflect deep-marine shales, coastal-plain topsets, alluvial plain, or lacustrine facies.
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
C-Dwn/Ob
44
45
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
46
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
47
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin
48
Sequence Stratigraphy Department of Earth Sciences National Central Univ., Taiwan Prepared by Dr. Andrew T. Lin