What can you really see in seismic?

Module 03  ·  Vertical resolution in practice

Can you see the fault?

A fault runs into the same limit as a thinning bed, with vertical offset in place of bed thickness. Push the throw down and at some point the reflectors stop showing a step, while the fault remains in the model. A quiz mode hides the answer and randomizes the case.

Why throw runs into the bed-thickness limit

A reflector cut by a fault gives two copies of the same event, one on each side, offset in time by the throw. That is the same geometry as the wedge in module 02 — two reflections a short time apart — so it fails at the same point: somewhere around a quarter of a wavelength.

Whether a fault is visible therefore depends on velocity, frequency and noise as much as on the size of the offset.

What happens below the limit

Below the limit the section gives no indication that anything is missing. The reflectors run through continuously and the fault leaves no visible step.

The reflections either side of the fault do still interfere, and that distorts the shape of the waveform. The distortion is small and hard to see by eye, and it is what coherence and similarity attributes are computed to detect.

Two questions follow: whether the fault is visible, and if it is not, whether it is absent or below the resolution limit.

Quiz mode

The fault

The layers

The data

Display

From rock layers to seismic

The two traces come from the dashed lines, one either side of the fault

1 The rock layers impedance  ·  pale = low, dark = high

2 The wavelet

3 The seismic that comes out layers convolved with the wavelet

4 Two traces

Wavelength λ = V / f80 m
λ/4 — seen as a step20 m
λ/8 — still detectable10 m
Throw as a fraction of λλ/4.0
Throw in two-way time14.3 ms
Beds in the section9
Traces across section120
Dip as drawn (true 65°)86°
Called correctly0 / 0
VerdictResolved

Try this — and what to watch for

  1. Walk the throw down. The point: it puts a number on where your eye gives up. Start at 40 m and reduce it. Where does the offset stop reading as a clean step in the reflectors? Somewhere close to λ/4. Above that you can measure the throw off the section. Below it you can only say that something happened.
  2. Keep going, and watch the trace panel. The point: this is what an attribute sees that you do not. Take the throw down to 5 m and compare the two traces. The step has gone from the section, but the traces still differ. Reflections interfering across the fault distort the waveform. Coherence and similarity attributes are built to find that residual, which is why an attribute can flag a fault that looks like nothing on the amplitude display.
  3. Trade frequency against throw. The point: the same fault can be present or absent depending on the vintage of the data. Fix the throw at 12 m and sweep frequency from 15 Hz to 70 Hz. The geology never changes, and the fault goes from invisible to obvious. Whether you can see a fault depends on your bandwidth as much as on the rocks.
  4. Make it deep. The point: it explains why deep fault maps are always simpler than shallow ones. Raise velocity to 5000 m/s and drop frequency to 15 Hz, roughly the change from 1 s to 4 s two-way time. λ goes past 300 m, and a 40 m throw that was obvious at shallow depth becomes marginal. Deep fault maps are less complete than shallow ones for reasons that have nothing to do with the interpreter.
  5. Coarsen the acquisition. The point: resolution is limited sideways as well as vertically. Push trace spacing to 50 m with a shallow fault dip. The fault plane turns into a staircase. Trace spacing limits lateral detail the same way bandwidth limits vertical detail.
  6. Push the dip past 90°. The point: normal and reverse faults look alike on seismic; the sense of slip is your interpretation, not the data's. Take it from 65° to 115° and watch the layer model. The plane now leans the other way, so the block that drops is the footwall rather than the hanging wall and the fault becomes reverse. A normal fault cuts section out at the fault; a reverse fault repeats it, and you can count the doubled beds in the model. On the seismic both look like an offset, which is the problem: sense of slip is a geological interpretation laid on top of what the data actually shows.
  7. Thin the beds. The point: most reservoir intervals are tuned, and tracking a fault through tuned reflections is far harder. Set the thinnest bed to 8 m and the thickest to 20 m, with frequency at 25 Hz. Count the beds in the layer model, then count the reflectors in the seismic below it. They no longer match. The beds are below tuning, so the reflectors merge and interfere. Tracking a fault through tuned reflections is much harder than tracking one through clean, well-separated ones, and most reservoir intervals look like the former.
  8. Run ten quiz cases on Mixed, then ten on Hard. The point: it tests the limit against your own judgement, without the answer visible. Keep score honestly. Roughly a third contain no fault at all, and throws are drawn as a fraction of the wavelength rather than in meters, so most of them sit near the limit instead of being comfortably above it. On Hard, only about one case in twenty is a clean resolved offset. Calling "no fault" on a throw below λ/16 is scored as correct, because past the detection limit that is the defensible answer. Scores are typically lower than expected, and correlated noise accounts for many of the false picks.

If you remember nothing else

  • Throw and bed thickness hit the same limit. The λ/4 and λ/8 numbers from the wedge model apply here without modification.
  • A fault you cannot see is not necessarily absent. Sub-resolution throw still distorts the waveform, and that distortion is what seismic attributes are designed to pick up.
  • Noise sets the practical limit. Textbook resolution limits assume clean data. At 3:1 you will miss faults you would pick easily at 15:1, and you will call a few that aren't there.
  • Two vintages over the same field will map different numbers of faults. That is usually a difference in bandwidth rather than a disagreement about the geology.
Method, assumptions, and where to read more

A stack of horizontal reflectors is offset across a planar normal fault. Bed thicknesses are drawn at random from the range you set and converted to two-way time using the interval velocity, so moving the velocity slider stretches the same geology rather than generating new geology. Each trace is the reflectivity convolved with a zero-phase Ricker wavelet. The noise is band-limited (random reflectivity through the same wavelet) and correlated from trace to trace, so it resembles noise on a migrated section rather than static. Signal-to-noise is quoted as reflection-coefficient amplitude over noise RMS.

Fault dip runs from 5° to 175°. Below 90° the plane leans one way and the block that drops is the hanging wall, giving a normal fault that omits section. Above 90° the plane leans the other way, the same block is now the footwall, and the fault is reverse and repeats section. Both come out of one construction: the two blocks translate rigidly past each other, and a boundary is drawn wherever it ends up. Reverse faults with large throw at low angles produce heavily repeated sections, which is correct but can look surprising.

This is a post-migration cartoon. It leaves out fault-plane reflections, tip diffractions, drag folding, amplitude changes from juxtaposition, and the smearing that imperfect migration puts on every real fault. Those effects make faults easier to see in some settings and harder in others, but none of them removes the resolution floor, which is what this module is about.

The two panels show different quantities. The layer model is shaded by acoustic impedance, a property of the rock: pale for low, dark for high. The seismic is colored by reflection coefficient, a property of the boundary between two rocks, so color appears only where impedance changes and a thick uniform bed is blank in the middle.

Polarity follows the SEG convention. Under normal polarity, the positive standard used in the United States, an increase in acoustic impedance is a positive reflection coefficient and is displayed as a peak: black on a gray section, blue on the default color map. A decrease is a trough, shown red. Reverse polarity flips that, so an impedance increase becomes a trough. Either way the color map itself runs trough to peak and never changes; polarity acts on the traces, not on the colors.

The practical check is that a pale-over-dark boundary in the layer model — impedance increasing downward — should line up with a blue event in the section under normal polarity, and a red one under reverse. The legend swatches are sampled through the current polarity and color map, so they stay correct in grayscale and in colorblind-safe mode too.

The display is vertically exaggerated, as nearly all seismic displays are. The stat strip gives the fault dip as drawn alongside its true dip so you can compare the two. Fault dips read off a section are routinely too steep for this reason.

Next up

04 · How small a feature can I map in plan view?

Everything so far has been about resolution downward — how thin a bed, how small a throw. The next module turns the problem on its side and asks how small a feature can be before it disappears in map view. The answer involves a patch of ground the size of several football fields, and one processing step that shrinks it.

Open module 04 → All modules
← The wedge model Next: horizontal resolution →