Fowler’s velocity-independent prestack migration
![]() | |
Series | Investigations in Geophysics |
---|---|
Author | Öz Yilmaz |
DOI | http://dx.doi.org/10.1190/1.9781560801580 |
ISBN | ISBN 978-1-56080-094-1 |
Store | SEG Online Store |
We now restate the underlying principle for migration velocity analysis:
- Starting with the prestack volume of data P(y, h, t) in midpoint y, offset h and two-way event time t in the unmigrated position, create a velocity cube — volume of data P(y, v, τ) in midpoint y, migration velocity v and two-way event time τ in the migrated position. Within the context of time migration, the output time variable τ is related to depth by way of vertical stretch: τ = 2z/v.

Although the velocity cube can be created by means of some of the migration velocity analysis techniques described in this section, a variation of the method by Fowler [1] is particularly efficient and elegant. First, we review Fowler’s method to create the velocity cube. Refer to the moveout equation (1) and recall that stacking velocity vstk is dip dependent:
( )
( )
where
( )
Use equation (11) to establish a relationship between the dip-dependent stacking velocities vstk and the dip-independent DMO velocities vDMO — velocities estimated from dip-moveout-corrected data:
( )
( )
This equation is the basis for Fowler mapping of vstk to vDMO. Note that the process is applied to data in the frequency-wavenumber domain. The Fowler mapping is then followed by constant-velocity Stolt mapping (migration principles and D.7) to map the DMO velocities vDMO to migration velocities vmig.
Stolt migration of the dip-moveout-corrected data volume in the Fourier transform domain P(ky, ω0; vDMO) involves, first, mapping from ω0 to ωτ for a specific ky by using the dispersion relation of equation (4-24b) recast as
( )
where the relationship ωτ = vmigkz/2 is used. The output of this mapping P(ky, ωτ; vmig) is then scaled by the quantity S given by equation (4-24c) recast as
( )
Again, the relationship ωτ = vmigkz/2 is used to obtain equation (47) from equation (4-24c).
Figure 5.4-35 describes a flowchart for creating the migration velocity volume P(y, vmig, τ) by Fowler and Stolt mapping.
- Start with data P(y, h, t) in coordinates of midpoint y, offset 2h and event time t in the unmigrated position, and create constant-velocity stack (CVS) volume P(y, vstk, tn) using a range of velocities vstk, where tn is the event time after constant-velocity normal moveout correction.
- Apply 2-D Fourier transform to obtain the CVS cube P(ky, vstk, ω) in the frequency-wavenumber domain, where ky and ω are the Fourier transform variables associated with the variables y and tn.
- Sort the CVS volume P(ky, vstk, ω) into a set of constant-velocity sections P(ky, ω; vstk).
- Perform the Fowler mapping based on equation (45) on each of the velocity sections so as to obtain the DMO velocity volume P(ky, ω0; vDMO).
- Migrate each of the constant-velocity sections of the DMO velocity volume by performing the Stolt mapping based on equations (46) and (47) so as to obtain the migration velocity volume P(ky, ωτ; vmig).
- Apply 2-D inverse Fourier transform to obtain the migration velocity volume P(y, τ, vmig).
A variation of Fowler’s sequence described above involves creating the CVS cube directly from DMO-corrected data.
- Start with data P(y, h, t) in coordinates of midpoint y, offset 2h and event time t in the unmigrated position, and apply DMO correction followed by inverse NMO correction.
- Create constant-velocity stack (CVS) volume P(y, vDMO, t0) using a range of velocities vDMO, where t0 is the event time after constant-velocity normal moveout correction (Figure 5.4-36a).
- Sort the CVS volume P(y, vDMO, t0) into a set of constant-velocity stacked sections P(y, t0; vDMO).
- Apply 2-D Fourier transform to obtain the CVS cube P(ky, vDMO, ω0) in the frequency-wavenumber domain, where ky and ω are the Fourier transform variables associated with the variables y and t0.
- Sort the CVS volume P(ky, vDMO, ω0) into a set of constant-velocity sections P(ky, ω0; vDMO).
- Migrate each of the constant-velocity sections of the DMO velocity volume by performing the Stolt mapping based on equations (46) and (47) so as to obtain the migration velocity volume P(ky, ωτ; vmig).
- Apply 2-D inverse Fourier transform to obtain the migration velocity volume P(y, τ, vmig) (Figure 5.4-36b).
Figure 5.4-36 (a) A velocity cube computed from DMO-corrected gathers; (b) the same data as in (a) after Stolt migration of each of the constant-velocity panels.
Figure 5.4-37 (a) A time slice from the migration velocity cube shown in Figure 5.4-36b, (b) the same time slice as in (a) with velocity picking. Red and dark blue curves correspond to primaries and multiples, respectively.
Figure 5.4-38 (a) Velocity strands picked from the time slices of the migration velocity cube shown in Figure 5.4-36b as demonstrated in Figure 5.4-37b; (b) the color-coded velocity field computed by interpolation of the velocity strands shown in (a).
Figure 5.4-39 (a) The prestack time-migrated section extracted from the migration velocity cube of Figure 5.4-36b coincident with the surface that corresponds to the velocity field in Figure 5.4-38b, (b) the same section as in (a) with the viewing axis perpendicular to the plane of the paper.
Figure 5.4-41 (a) A migration velocity volume, (b) the same volume viewed end-on showing the vertical variation in rms velocity at a given CMP location.
Figure 5.4-42 (a) A time slice from the migration velocity volume shown in Figure 5.4-41a with the interpretation of the rms velocity as a function of midpoint location, (b) the migration velocity volume as in Figure 5.4-41a with the picked velocity strands as in (a).
Figure 5.4-43 (a) The color-coded optimum rms velocity surface extracted from the migration velocity volume shown in Figure 5.4-41a, (b) the same view as in Figure 5.4-41b with the intersection of the rms velocity surface showing the vertical velocity variation at a given CMP location.
Figure 5.4-44 (a) The image surface associated with prestack time migration that was extracted from the migration velocity volume shown in Figure 5.4-41a by sculpting the amplitudes over the rms velocity surface shown in Figure 5.4-43a, (b) the image surface as in (a) after collapsing it onto a 2-D section associated with prestack time migration.
The migration velocity volume P(y, t, vmig) shown in Figure 5.4-36b can be visualized and interpreted to derive a migration velocity field. For spatial consistency, velocity picking should be done on time slices from the migration velocity volume as shown in Figure 5.4-37. The resulting velocity strands shown in Figure 5.4-38a are interpolated to create the migration velocity field shown in Figure 5.4-38b. This velocity field then can be used to extract from the volume the section that corresponds to prestack time migration as shown in Figure 5.4-39. An enlarged view of this section is shown in Figure 5.4-40. Note the excellent imaging of the steeply dipping fault planes which conflict with the gently dipping strata.
Alternatively, the plane of (vmig, τ) for each midpoint y can be inverse transformed to the plane of (h, τ) associated with the common-reflection-point gather derived from prestack time migration. This is then followed by conventional normal-moveout correction and stacking. The resulting section, again, represents the image from prestack time migration.
The data example shown in Figure 5.4-36 demonstrates the use of the migration velocity volume in deriving a high-fidelity image of fault planes from prestack time migration (Figure 5.4-40). The data example shown in Figure 5.4-41 demonstrates the use of the migration velocity volume in imaging steeply dipping event. The migration velocity volume was created using the procedure described above. Specifically, the DMO-corrected CMP gathers were first NMO-corrected using a range of constant velocities and a CVS volume was created. Next, each constant-velocity stacked section was migrated using the Stolt method and the constant velocity associated with that section. The resulting migration velocity volume is shown in Figure 5.4-41. Note the vertical variation in velocities on the end-on view of the volume that represents the plane of velocity versus time.
The migration velocity volume is interpreted by picking the primary velocity trend from selected time slices as shown in Figure 5.4-42a. Note the lateral variation in velocities which is captured by continuous picking along the midpoint axis. By interpolating the velocity strands resulting from the interpretation of selected time slices (Figure 5.4-42a), an rms velocity surface is generated. The picked velocity strands are shown in Figure 5.4-42b embedded within the migration velocity volume, and the rms velocity field is shown in Figure 5.4-43a as a color-coded surface extracted from within the migration velocity volume. A quality control of the rms velocity surface can be made by intersecting it with the cross-sections of the migration velocity volume at selected CMP locations as shown in Figure 5.4-43b.
The prestack time-migrated section is a by-product of the migration velocity analysis described here. Specifically, the image surface associated with the prestack time migration is obtained by sculpting the amplitudes from within the migration velocity volume over the rms velocity surface as shown in Figure 5.4-44a. The conventional 2-D display of the prestack time-migrated section is then created by simply collapsing the sculpted image surface onto a 2-D plane (Figure 5.4-44b). A close-up view of the prestack time-migrated section shows accurate imaging of the steeply dipping event (Figure 5.4-45).
References
- ↑ Fowler (1984), Fowler, P., 1984, Velocity-independent imaging of seismic reflectors: 54th Ann. Internat. Mtg., Soc. Explor. Geophys., Expanded Abstracts, 383–385.
See also
- Prestack Stolt migration
- Common-offset migration of DMO-corrected data
- Prestack Kirchhoff migration
- Velocity analysis using common-reflection-point gathers
- Focusing analysis