Click to view images of illustrations in PDF format.
Normal MoveOut and Dip MoveOut*
By
Christopher L. Liner1
Search and Discovery Article #40037 (2002)
*Adapted for online presentation from two articles by the same author in AAPG Explorer (April and May, 1998), respectively entitled “Normal MoveOut, an Effect, Process” and “Dip MoveOut Just Isn’t ‘Normal.’” Appreciation is expressed to the author and to M. Ray Thomasson, former Chairman of the AAPG Geophysical Integration Committee, and Larry Nation, AAPG Communications Director, for their support of this online version.
1Department of Geosciences, University of Tulsa ([email protected])
|
Click here for sequence of Figures 2, 5, 6.
Click here for sequence of Figures 2, 5, 6.
Click here for sequence of Figures 2, 5, 6.
Petroleum seismology is, and always has been, changing very quickly. You might have heard whispers about exotic topics like crosswell tomography, wavelet transforms, cluster analysis, texture segmentation impedance inversion, geostatistical estimation, etc. So why in this high-tech age is someone writing about something as ancient as normal moveout? The answer involves the importance of understanding fundamental concepts, the natural lead-in that normal moveout provides to the juicier topic of dip moveout, and a chance to do it without any equations. Normal moveout has two meanings – it is both:
By itself, the term “moveout” goes back to the earliest days of reflection seismology in the 1910s. In those days a seismic shot consisted of a source (dynamite) sending waves into the earth to bounce around and return to a few geophones. The data were preserved as wiggly lines (traces) on a rotating drum of paper, or as dark lines on a photographic record. The human eye is wonderfully adept at seeing patterns and relationships in very confusing data, e.g. recognizing a face across a crowded room full of strangers. Early seismic records were like that – lots of noise, not much signal. But the signals were there, and skilled interpreters could recognize them. Some of these signal events came in straight lines across the traces; others formed curves. But whatever the shape, each kind of signal showed a delay from trace to trace as we move away from the source – and thus was known as moveout. You
can talk about Normal MoveOut (NMO) all day long without mentioning A shot record is the collection of seismic traces generated when one source shoots into many receivers, as shown in Figure 1. In this example, the upper black line is the acquisition surface and the lower one is the reflector. Dots below the reflector show subsurface reflection points. Halfway between the source and a receiver is a point on the ground called the midpoint. These are shown as black dots above the acquisition surface Where there is no dip, the midpoint is directly above the reflection point. As the offset (source to receiver distance) increases, so does the travel time from source to receiver. This characteristic delay of reflection times with increasing offset is called normal moveout. In Figure 2 (page 26), reflections can be seen in real data along with other kinds of events. There are receivers on both sides of the shot in this case. The right side has been marked-up to identify different kinds of events – direct arrivals (p-wave, s-wave, air wave, surface wave), head waves and (a few) reflections. The left side is uninterpreted. The reflection events have a hyperbolic shape characteristic of normal moveout. So now we know something about normal moveout, the effect. What about normal moveout, the process? For this we will use the acronym NMO (Normal MoveOut) since for most people this term implies the process, not the effect. What is NMO? The short answer is: A seismic processing step whereby reflection events are flattened in a common midpoint gather in preparation for stacking. (If this makes sense then you can move along to the second part of this article. Otherwise, read on.) A seismic line is generated by “rolling” the shot and receivers forward a certain distance and firing again. As shown in Figure 3, this generates a second shot record which partially overlaps the first. Note that six of the seven reflection points from the blue shot were also reflection points for the red shot. This occurs by design and is called Common MidPoint (CMP) shooting. As the shots roll along, there will be many source-receiver pairs with the same CMP location, and the CMP fold is the total number of traces that live at any given CMP (Figure 4). CMP fold can vary from as few as six (low-fold land 3-D) up to several hundred (2-D marine). The reason for gathering multifold data is that we get redundant information about the reflection point down in the earth, and this redundancy can be used to reduce noise and create a more reliable image. Our goal is to eventually process all these traces as a family and add them together (CMP stack) to make one trace that lives at this CMP location. NMO is aimed at removing the hyperbolic curvature in reflection events. Basically, it is removing the effect of offset. If this is done properly, then the reflection should come in at the same time for all offsets (since we have removed any travel time delay due to offset). In short, reflection events should be flat after NMO. Figure 5 shows the data after NMO processing (and air wave removal). In this case, we see the events are pretty well flattened by NMO, but there are a couple of interesting areas. The blue box shows some disturbing behavior along a flattened reflection. This has nothing at all to do with NMO, but is related to lateral changes in the near surface layers (termed a static problem). The red box shows what NMO does to the direct arrivals. Since these were linear and not hyperbolic, NMO has not flattened them. Also, note how fat (not flat) these events look after NMO. This is because NMO actually operates by stretching the trace – and the shallower something is, the more it stretches. Since our goal is to eventually flatten all these traces and add them together to make one trace, keeping this kind of stuff would wipe out shallow reflections. It needs to go. We get rid of these events by muting – which is nothing more than replacing the offending data with zeros. We could do this by hand, but a seismic line may contain many thousands of CMP gathers. It turns out we can let NMO itself do the muting for us. The idea is to keep track of how much stretch NMO is doing to the trace. The stretch changes down the trace – biggest at the top, smallest at the bottom. So the idea of an NMO stretch mute is to set a limit on how much stretch we are going to allow. If the stretch gets bigger than our limit, then the data values are replaced by zeros in that part of the trace. Figure 6 shows the result of allowing a 25 percent stretch in the NMO process. At this point more than one astute reader is saying, “Yeah, yeah, you academic types are all alike – ignoring the real world. What about dipping beds?” Figure 7 gives a hint of things to come. The basic problem is that when the bed dips, midpoints are not vertically above reflection points. Furthermore, the reflection points become unevenly spaced along the reflector. The fact is that NMO has a constant dip assumption built in. If every bed were dipping at the same angle, say 23 degrees, it turns out we could do NMO just fine. The real problem comes where there are many dips in the subsurface. In this case, NMO acts like a dip filter – preferentially passing some dip the processor chooses while suppressing others. This was particularly notorious in places like the Gulf of Mexico, where gently dipping beds often meet steep salt domes. You could do NMO in such a way as to get a good image of one or the other, but not both. This was the situation up until about 1981 or 1982. Then something new came on the scene. It was called – you guessed it – dip moveout. In
the field, a source and receiver are located, say, 3,000 meters apart (the
offset) and a trace is recorded. This is pre-stack data. From the notes above,
we know that NMO lives to remove offset from pre-stack data. In the computer, we
adjust this trace to simulate one that would have been recorded with no effect
at a point half-way between the source and receiver (the midpoint). This new
trace is a zero-offset trace, and adding (stacking) all such traces that live at
this midpoint yields a stack trace. All the stack traces plotted side-by-side
form the stack section, which is raw material for post-stack Figure
8 illustrates the NMO idea. On the left is a field trace with some arbitrary
offset and one reflection event. NMO assumes the reflection comes from a
horizontal interface is shown in the middle (inset) figure. This is an important
and restrictive assumption. The NMO But what if the interface is not horizontal? It is tempting to think that the reflector could be anywhere and still be consistent with the observed travel time. But this is not the case. Assume that the original travel time is 1.0 second, and we know the velocity is 3,000 m/s. The total distance traveled has got to be 3,000 meters. So all valid reflector positions have one thing in common: The total distance from source to reflection point to receiver is constant, namely 3,000 meters. Thinking back to Geometry 101, this is just the definition of an ellipse with the source and receiver at each focus. Figure 9
shows such an ellipse. Remember the goal is to remove offset and thus create a
zero-offset section. NMO gives one of many possibilities; From
Figure 8, we saw that NMO is a process that takes one trace in and gives one
trace out. In
Figure
10, we have a panel of data containing only two spikes of amplitude on
one trace (left panel). The other traces are there, but empty. NMO moves the
spike up on the same trace (middle panel). Now
here is some magic. By creating the Since
Here is the bottom line: If structure and velocity variation are not too nasty in an area, we can get away with a traditional processing sequence: NMO + In this equation,
“+” means “followed by.” However, if things get really tough down there
(e.g., subsalt), this sequence breaks down and fails to give a good image. In
this case we are compelled to do one grand process called pre-stack In
fact, PreStackMigration = NMO + ? +Stack + PostStackMigration. The
“?” turned out to be |