Topics in Dip-Moveout Correction and Prestack Time 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 |
Contents
- 1 E.1 Reflection point dispersal
- 2 E.2 Equations for DMO correction
- 3 E.3 Log-stretch DMO correction
- 4 E.4 The DMO ellipse
- 5 E.5 Nonzero-offset traveltime equation
- 6 E.6 Prestack frequency-wavenumber migration
- 7 E.7 Velocity analysis by wavefield extrapolation
- 8 References
- 9 See also
- 10 External links
E.1 Reflection point dispersal
We start the analysis using the geometry of Figure E-1. Place the source S at the origin (0, 0), and denote the coordinates of the points of interest as N : (x_{0} = SK, z_{0} = KN), H : (x_{1} = SB, z_{1} = BH), F : (x_{2} = SA, z_{2} = AF), and R : (x_{3} = SQ, z_{3} = QR). Our objective is to compute the coordinates for the normal-incidence reflection point N : (x_{0}, z_{0}) associated with midpoint M, and the reflection point R : (x_{3}, z_{3}) associated with the nonzero-offset source receiver separation x = SG.
First, we compute the coordinates for N : (x_{0}, z_{0}). From the geometry of Figure E-1, first, note that x_{0} = SK and z_{0} = KN. Then, from the triangle MKN, we have
which, by way of MK = x_{0} − x/2, yields
( )
Similarly, from the same triangle MKN, we have
which, by way of KN = z_{0}, yields
( )
Finally, combine equations (E-1a) and (E-1b) to get a relation between x_{0} and z_{0}
( )
Next, write the equation of the line that represents the dipping reflector in terms of the directional cosines and the normal d = SC to the line from the origin S
( )
Solve for z_{0} from equation (E-2a) and substitute into equation (E-2b)
( )
Note that the directional cosines satisfy the relation cos^{2} α + cos^{2} ϕ = 1. By making use of this relation, simplify equation (E-3) to get the final expression for the coordinate x_{0}
( )
Finally, substitute x_{0} back into equation (E-2a) to get the expression for the coordinate z_{0}
( )
The distance d = SC is given by
which can be rewritten using the geometry of Figure E-1 as follows:
Substitute the definitions SM = x/2 and D = MN, where D is the distance along the normal-incident raypath from the point N on the reflector to the midpoint location M, to obtain the relation between d and D
( )
Make the further substitution of equation (E-5) into equations (E-4a) and (E-4b), and recall that cos^{2} α + cos^{2} ϕ = 1 to derive the expressions for the coordinates of N : (x_{0}, z_{0}) in terms of D
( )
and
( )
To derive the expressions for the coordinates of R : (x_{3}, z_{3}), we note that R is the intersection of the line defined by the reflecting interface, which yields the relation
( )
and the line SRF, which yields the relation
( )
where (x_{2}, z_{2}) are the coordinates of F. From the geometry of Figure E-1, first, note that x_{2} = SA, z_{2} = AF, and x_{3} = SQ and z_{3} = QR.
Solve equation (E-7b) for x_{3} and substitute into equation (E-7a) to get
( )
The coordinates for F : (x_{2}, z_{2}) are derived in Section C.3 (equations C-18c,d)
( )
and
( )
Substitute equations (E-9a) and (E-9b) into equation (E-8), and replace d via equation (E-5) to obtain, after some involved algebra, the expression for z_{3}
( )
Back substitution of equation (E-10a) into equation (E-7a) then yields the expression for x_{3}
( )
Knowing the coordinates of N : (x_{0}, z_{0}) and R : (x_{3}, z_{3}), the distance Δ = NR can now be computed
( )
Substitute equations (E-10a) and (E-10b) for the coordinates of R : (x_{3}, z_{3}), and equations (E-6a) and (E-6b) for the coordinates of N : (x_{0}, z_{0}) into equation (E-11a) and carry out the algebraic steps to get the final expression
( )
Equation (E-11b) represents the reflection point smear, otherwise known as reflection point dispersal ^{[1]} — the distance between the normal-incidence reflection point N associated with the zero-offset raypath from midpoint M, and the reflection point R associated with the nonzero-offset raypath for a source-receiver pair separated by an offset x. For the case of zero-offset x = 0, reflection point smear Δ vanishes.
Define the two-way zero-offset time t_{0} at midpoint M as the time associated with the normal-incidence raypath MN by the relation vt_{0}/2 = D. Then use this relation in equation (E-11b) to obtain the expression for reflection point dispersal Δ in terms of t_{0}
( )
This equation is applicable to the general case of three-dimensional geometry (Figure 3.1-14), where α is the angle between the normal to the dipping reflector surface and the direction of the profile line ^{[2]}.
For the 2-D geometry of the dipping reflector Figure E-1, note that
( )
where ϕ is the dip angle of the reflector. Hence, equation (E-12) is written in terms of the reflector dip angle ϕ as follows:
( )
( )
( )
Equation (E-14), by way of equation (6), is the same as equation (10) of the main text with offset x = 2h.
E.2 Equations for DMO correction
In this section, we shall derive the traveltime equation for dip-moveout (DMO) correction by using the geometry of Figure 5.1-1. Start with equation (C-22b) of Section C.3 which defines the traveltime t from source location S to the reflection point R to the receiver location G, written in prestack data coordinates
( )
where 2h is the offset, v is the medium velocity above the reflector, ϕ is the reflector dip, and t_{0} is the two-way zero-offset time at midpoint location y_{n}. Dip-moveout correction is preceded by normal-moveout (NMO) correction using the dip-independent velocity v
( )
where t_{n} is the event time after the NMO correction. To relate event time t_{n} after the NMO correction and event time t_{0}, first, write equation (E-15a) as
( )
Set the right-hand sides of equations (E-15b) and (E-16) equal and simplify the result to obtain the dip-dependent moveout equation
( )
Equations (E-15b) and (E-17) suggest that moveout correction can, in principle, be performed in two steps:
- Apply a dip-independent moveout correction using equation (E-15b) to map event time t to event time t_{n}.
- Apply a dip-dependent moveout correction using equation (E-17) to map event time t_{n} to event time t_{0}.
This two-step moveout correction is equivalent to the one-step moveout correction using equation (E-15a) to map event time t directly to event time t_{0}.
Our goal, however, is to map event time t to τ_{0} — two-way zero-offset time, not at midpoint location y_{n} associated with the source-receiver pair S − G, but at midpoint location y_{0} associated with the normal-incidence reflection point R. This mapping requires coordinate transformation of moveout-corrected prestack data P_{n}(y_{n}, t_{n}; h) to P_{0}(y_{0}, τ_{0}; h). Therefore, we need to compute the coordinates (y_{0}, τ_{0}) in terms of the coordinates (y_{n}, t_{n}).
From the geometry of Figure 5.1-1, note that
( )
where Δ = NR is the distance along the reflector between the normal-incidence reflection point N associated with the zero-offset raypath from midpoint y_{n} and the reflection point R associated with the nonzero-offset raypath for a source-receiver pair separated by an offset 2h.
Adapt equation (E-14) of Section E.1 for the prestack data coordinates x = 2h
( )
and substitute into equation (E-18) to obtain
( )
Solve equation (E-17) for t_{0} and write the result in the following form
( )
or
( )
where
( )
Now, substitute equation (E-22) into equation (E-20) to get the final expression for y_{0}
( )
Also from the geometry of Figure 5.1-1, note that
( )
Substitute equation (E-19) for Δ and carry out the required algebra to get the desired expression for τ_{0} in terms of t_{0}
( )
Now, substitute equation (E-22) into equation (E-26)
( )
Finally, simplify equation (E-27) by way of equation (E-23) to obtain the desired expression for τ_{0} in terms of t_{n}
( )
We now remind ourselves of our objective to transform the normal-moveout-corrected prestack data P_{n}(y_{n}, t_{n}; h) from y_{n} − t_{n} coordinates to y_{0} − τ_{0} coordinates so as to obtain the dip-moveout-corrected data P_{0}(y_{0}, τ_{0}; h). The transformation is done using equations (E-24) and (E-28) which relate the input data coordinates y_{n} − t_{n} to output data coordinates y_{0} − τ_{0}. Note, however, equations (E-24) and (E-28) require knowledge of the reflector dip ϕ. To circumvent this requirement, we use the relation from Section D.1
( )
which states that the reflector dip ϕ can be expressed in terms of wavenumber k_{y} and frequency ω_{0}, which are the Fourier duals of midpoint y_{0} and event time τ_{0}, respectively. The variables y_{0} and τ_{0} defined by equations (E-24) and (E-28), respectively, are associated with the zero-offset section after DMO correction. By way of equation (E-29), these equations are recast explicitly independent of reflector dip as
( )
and
( )
where A of equation (E-23), by way of equation (E-29), is of the form
( )
Since we have switched to the Fourier transform domain in our analysis, our objective now is to compute the dip-moveout-corrected data P_{0}(k_{y}, ω_{0}; h) in the transform domain. Thus, start with the 2-D Fourier transform of P_{0}(y_{0}, τ_{0}; h)
( )
A wavefield is invariant under a coordinate transformation; hence, P_{0}(y_{0}, ω_{0}; h) = P_{n}(y_{n}, t_{n}; h). Use the transform relations given by equations (E-30) and (E-31) and the invariance relation to write equation (E-33) in terms of the normal-moveout-corrected data P_{n}(y_{n}, t_{n}; h) ^{[3]}
( )
From equation (E-30) it follows that
( )
Square both sides of equation (E-31), substitute equation (E-32) for the quantity A and simplify to get
( )
Next, differentiate the result
( )
Rearrange the terms in equation (E-32) to the form
( )
and use in equation (E-36b) together with equation (E-31) for the τ_{0} in the denominator to obtain
( )
Finally, substitute equations (E-35) and (E-37) into equation (E-34)
( )
Now, we turn our attention to the phase term in equation (E-38)
Rearrange the terms, first, as follows:
then, by way of the expression for A as in equation (E-32)
( )
Return to equation (E-38) and substitute equation (E-39) for the terms in the exponential
( )
Use the Fourier transform integral
( )
to rewrite equation (E-40) in the form
( )
Equation (E-42) describes the dip-moveout correction process which transforms the normal-moveout-corrected prestack data with a specific offset 2h from y_{n} − t_{n} domain to y_{0} − τ_{0} domain. Referring to Figure 5.1-1, we remind ourselves that y_{n} − t_{n} coordinates are associated with event time t_{n} at midpoint location y_{n}, and y_{0} − τ_{0} coordinates are associated with event time τ_{0} at midpoint location y_{0} that corresponds to the normal-incidence reflection. The lateral excursion applied by dip-moveout correction is Δy_{DMO} = |y_{n} − y_{0}| and the vertical excursion is Δt_{DMO} = t_{n} − τ_{0}.
Once dip-moveout correction is applied, the data are inverse Fourier transformed
( )
A flowchart of the dip-moveout correction in the frequency-wavenumber domain described above is presented in Figure 5.1-2.
E.3 Log-stretch DMO correction
A computationally efficient DMO correction can be formulated in the logarithmic time domain ^{[4]}^{[5]}^{[6]}^{[3]} ^{[7]}. The log-stretch time variable enables linearization of the coordinate transform equation (E-28), and as a result, the DMO correction is achieved by a simple multiplication of the input data with a phase-shift operator in the Fourier transform domain.
Define the following logarithmic variables that correspond to the time variables τ_{0} and t_{n} of equation (E-28):
( )
and
( )
Hence, the inverse relationships are given by
( )
and
( )
In Section E.2, we derived the equations for DMO correction to transform the normal-moveout-corrected data from (y_{n}, t_{n}) coordinates to (y_{0}, τ_{0}) coordinates. The objective here is to derive equations for DMO correction in the log-stretch coordinates (y_{0}, T_{0}).
Square both sides of equation (E-28)
( )
and substitute equation (E-23) for A
( )
Here, we have to make the crucial assumption that A is close to unity. Refer to equation (E-23) and note that such an assumption implies one of the following: small h (offset), small ϕ (dip), large t_{n} (time after moveout) or large v (velocity). As a result, equation (E-46b) can be approximated as
( )
The slope dτ_{0}/dy_{0} measured on a zero-offset section in (y_{0}, τ_{0}) coordinates is given by
( )
Substitute equation (E-46d) into equation (E-46c)
( )
Apply the chain rule for differentiation
( )
and, by way of equation (E-45a), obtain the relationship
( )
Next, combine equations (E-29) and (E-46d) to get
( )
By analogy, we may write
( )
where Ω_{0} is the Fourier transform variable associated with the log-stretch time variable T_{0}. Finally, substitute equation (E-47d) into equation (E-47b) to get
( )
Now, return to equation (E-46e), and use equations (E-45a,E-45b) and (E-47e) to obtain the expression
( )
Take the logarithm of both sides and simplify to get the transform relation between the input log-stretch time variable T_{n} and the output log-stretch time variable T_{0}
( )
where
( )
We now turn our attention to equation (E-24) and derive the corresponding equation in the log-stretch domain. First, substitute equation (E-23) into equation (E-24) and approximate the quantity A as before
( )
Substitute equation (E-46d) into equation (E-51a)
( )
Then, make the substitutions from equations (E-45b) and (E-47e), drop the high-order term, and simplify to get
( )
Note from equation (E-48)
( )
Use this result in equation (E-51c) accompanied with the definition of equation (E-50) to obtain the final expression that relates y_{0} and y_{n} in log-stretch coordinates
( )
Equations (E-49), (E-50), and (E-53) correspond to equations (E-31), (E-32), and (E-30) in the log-stretch domain.
We remind ourselves of the objective to compute the dip-moveout-corrected data P_{0}(k_{y}, Ω_{0}; h) in the transform domain. Thus, start with the 2-D Fourier transform of P_{0}(y_{0}, T_{0}; h)
( )
A wavefield is invariant under a coordinate transformation; hence, P_{0}(y_{0}, T_{0}; h) = P_{n}(y_{n}, T_{n}; h), the latter being the normal-moveout-corrected data in the log-stretch domain. Use the transform relations given by equations (E-49) and (E-53) and the invariance relation to write equation (E-54) in terms of the normal-moveout-corrected data P_{n}(y_{n}, T_{n}; h) ^{[3]}
( )
From equation (E-53) it follows that
( )
and from equation (E-49), we have
( )
Substitute equations (E-56a,E-56b) into equation (E-55) and rearrange the terms in the exponential to obtain
( )
Assume that the Fourier transform variable Ω_{0} in the log-stretch domain is independent of T_{n}. Then, the double integral on the right-hand side represents the 2-D Fourier transform of the normal-moveout-corrected data in the log-stretch domain as
( )
Therefore, by way of equation (E-58), equation (E-57) takes the form
( )
Equation (E-59) describes the dip-moveout correction process which transforms the normal-moveout-corrected prestack data with a specific offset 2h from the log-stretch y_{n} − T_{n} domain to y_{0} − T_{0} domain. Note that the relationship of input P_{n}(k_{y}, Ω_{0}; h) to output P_{0}(k_{y}, Ω_{0}; h) given by equation (E-59) computationally is much simpler than that of equation (E-42). Once the phase-shift in the log-stretch domain given by equation (E-59) is applied to the input and the result is inverse Fourier transformed, the data are unstretched from (y_{0}, T_{0}) coordinates back to (y_{0}, τ_{0}) coordinates using the relationship expressed by equation (E-45a).
A variation of the phase-shift term in equation (E-59) is given by Notfors and Godfrey ^{[6]}. As in most log-stretch formulations of DMO correction, this reference assumes that under DMO correction the midpoint variable is invariant; hence, by way of equation (E-53), the first term in the exponential of equation (E-59) drops out. A further approximation, ln A_{e} = A_{e} − 1, and use of the definition for A_{e} given by equation (E-50) then lead to the following expression for DMO correction:
( )
Various implementations of the log-stretch technique applied to shot records are described by Biondi and Ronen ^{[8]}, Cabrera and Levy ^{[9]}, and Zhou ^{[7]}. The latter also includes a method in double-log-stretch domain, which involves stretching both in time and midpoint coordinates. See principles of dip-moveout correction for a practical implementation of the log-stretch DMO correction described here.
E.4 The DMO ellipse
The objective is to derive the traveltime trajectory associated with the dip-moveout correction operator by using the method of stationary phase. Insert from equation (E-42) the expression for P_{0}(k_{y}, τ_{0}; h) into equation (E-43)
( )
where the total phase is given by
( )
Substitute equation (E-32) for A
( )
The main contribution to the integration in equation (E-61) occurs when the phase in equation (E-62b) stays nearly constant. We therefore determine the variation of Φ with respect to variables ω_{0} and k_{y}
( )
and
( )
Then set each variation to zero. Rearranging the terms of the resulting expressions, we have
( )
and
( )
Then, square both sides of equations (E-64a,E-64b)
( )
and
( )
Finally, sum equations (E-65a,E-65b) to obtain
( )
Equation (E-66) describes an ellipse with the following properties:
- Semi-major axis in midpoint y_{0} direction: a = h.
- Semi-minor axis in time τ_{0} direction: b = t_{n}.
The ellipse of equation (E-66) in the y_{0} − τ_{0} plane describes the impulse response of a dip-moveout operator applied to a common-offset section with offset 2h.
E.5 Nonzero-offset traveltime equation
Prestack wave extrapolation is performed using the double-square operator which enables downward continuation of common-shot and common-receiver gathers (Section D.1). Stationary phase approximation to the double-square root operator yields the nonzero-offset traveltime equation (D-35) derived in Section D.2
( )
where y, h, t, z are midpoint, offset, two-way traveltime and depth coordinates and v is the medium velocity. Here, we shall show that this equation represents an ellipse (Figure E-2) in the y − z plane for a constant t, h, and v, and derive the parameters of this ellipse. The nonzero-offset two-way time is associated with the raypath from source S to reflection point R to receiver G as sketched in Figure E-2. The origin of the y − z plane coincides with midpoint M.
Square both sides of equation (E-67) to get
( )
Combine the second and third terms on the right-hand side and simplify the terms inside the square root
( )
Perform further algebraic manipulation to collect the terms in y and z
( )
Finally, normalize by the terms on the right-hand side and rearrange the terms in the denominators
( )
Equation (E-70) represents an ellipse in the y − z plane for a constant t with the following parameters (Figure E-2):
- Semi-major axis in midpoint y direction: a = vt/2.
- Semi-minor axis in depth z direction:
- Distance from center to either focus:
- Distance from one focus to a point on the ellipse to the other focus: vt.
The ellipse of equation (E-70) in the y − z plane describes the impulse response of a nonzero-offset migration operator applied to prestack data.
When equation (E-70) is specialized to the zero-offset case, h = 0, we get
( )
which describes a circle in the y − z plane for a constant t with a radius vt/2. This circle represents the impulse response of a zero-offset migration operator applied to poststack data.
When equation (E-67) is specialized to the zero-offset case, h = 0, we get
( )
which describes the well-known diffraction hyperbola in the y − t plane for a constant z.
Reduction of the double-square-root equation (E-67) to a single-square-root equation (E-71b) can also be achieved by defining an equivalent offset h_{e} such that ^{[10]}
( )
Solving for h_{e}, we obtain ^{[10]} ^{[11]}
( )
where t is the two-way nonzero-offset traveltime of equation (E-67).
Poststack time migration can be conceptualized either by way of a semicircle superposition using equation (E-72a) or a diffraction summation along the hyperbolic traveltime trajectory using equation (E-72b). Similarly, prestack time migration can be conceptualized either by way of semi-elliptical superposition using equation (E-70) or diffraction summation over the traveltime surface described by equation (E-67). The traveltime surface is known as Cheops’ pyramid Claerbout ^{[12]} and is illustrated in Figure E-3a. The result of summation of amplitudes over the pyramidal surface is placed at its apex. The question that is of practical importance is how to define the summation paths over this surface.
Four possible choices of summation paths over the pyramidal surface of equation (E-67) to perform prestack time migration are:
- Summation curves of constant offset: Consider a set of vertical cross-sections of the traveltime pyramid illustrated in Figure E-3a parallel to the midpoint axis as shown in Figure E-3b. Sum the amplitudes along each of the constant-offset table-top traveltime curves, independently, and place the result for each at the apex of the summation curve. The summation collapses the pyramidal surface onto a hyperbolic traveltime curve, which is formed by combining the apex points of the constant-offset curves. This hyperbolic traveltime curve is orthogonal to the constant-offset summation curves.
- Summation curves of constant time: Consider a set of horizontal cross-sections of the traveltime pyramid as shown in Figure E-3c ^{[13]}^{[14]}. Sum the amplitudes along each of the constant-time curves, independently, and place the result for each at the maximum offset on the summation curve. The summation, again, collapses the pyramidal surface onto the hyperbolic traveltime curve, which is formed by combining the maximum-offset points of the constant-time curves. This hyperbolic traveltime curve is orthogonal to the constant-time summation curves. The event associated with the resulting hyperbolic moveout trajectory is on the h_{e} − t plane of a common-scatter-point (CSP) gather ^{[10]}, where the equivalent offset h_{e} is given by equation (E-72b). The scatter point corresponds to the apex of the traveltime pyramid A_{0} in Figure E-3c.
- Summation curves of constant shot: Consider a set of vertical cross-sections of the traveltime pyramid as shown in Figure E-3d ^{[15]}. Sum the amplitudes along each of these constant-shot curves, independently, and place the result for each at the apex of the summation curve. The summation collapses the pyramidal surface onto the hyperbolic traveltime curve, which is formed by combining the apex points of the constant-shot curves. This hyperbolic traveltime curve is orthogonal to the constant-shot summation curves and is on the common-receiver plane that passes through the apex of the pyramid itself.
- Summation curves of constant angle: Consider a set of slanted cross-sections of the traveltime pyramid as shown in Figure E-3e ^{[16]}. These slanted traveltime curves are associated with constant angle of incidence (the slant-stack transform). Sum the amplitudes along each of the constant-angle curves, independently, and place the result for each at the apex of the summation curve. The summation collapses the pyramidal surface onto the hyperbolic traveltime curve which is formed by combining the apex points of the constant-offset curves. This hyperbolic traveltime curve is orthogonal to the constant-offset summation curves.
E.6 Prestack frequency-wavenumber migration
We start with prestack data in midpoint y and offset h coordinates, P(y, h, z = 0, t), where t is the event time in the unmigrated position, and perform 3-D Fourier transform
( )
where k_{y}, k_{h}, and ω are the Fourier duals of the variables of the midpoint-offset coordinates, y, h, and t.
Extrapolate the prestack data from the surface z = 0 to a depth z by the following extrapolation equation which we borrow from Section D.1
( )
The vertical wavenumber k_{z} is given by
( )
where, the double-square-root (DSR) operator takes the following form in midpoint-offset coordinates (equation D-22 of Section D.1)
( )
The variables Y and H are the normalized midpoint and offset wavenumbers, respectively
( )
and
( )
Assume a horizontally layered earth model associated with a vertically varying velocity function v(z). By inverse Fourier transforming equation (E-74), we have
( )
The imaging principle t = 0 then is applied to get the image volume P(y, h, z, t = 0),
( )
This is the equation for the prestack phase-shift method. Equation (E-79) involves integration over frequency ω and 2-D inverse Fourier transformation along midpoint y and offset h axes.
We now consider the special case of constant velocity v. Stolt ^{[18]} devised a prestack migration technique that involves an efficient mapping in the 3-D Fourier transform domain from temporal frequency ω to vertical wavenumber k_{z}.
First, combine equations (E-75) and (E-76)
( )
Square both sides
( )
and simplify
( )
Define
( )
and square both sides of equation (E-82). After some algebra, it follows that
( )
which can be rewritten as
( )
Now substitute equations (E-77a,E-77b) and (E-83) into equation (E-85), and simplify to obtain the final expression for the dispersion relation for prestack wave extrapolation
( )
By setting the offset wavenumber k_{h} = 0, we obtain the special case of zero-offset dispersion relation
( )
as in equation (D-85) with x replaced by y.
By keeping the wavenumbers k_{y} and k_{h} unchanged and differentiating equation (E-86), we get
( )
By setting the offset wavenumber k_{h} = 0, we obtain the special case for zero-offset
( )
as in equation (D-86) with x replaced by y.
When equations (E-86) and (E-88) are substituted into equation (E-79), we get
( )
Finally, sum over k_{h} to obtain the image at zero offset, h = 0
( )
This is the equation for constant-velocity prestack Stolt migration. It involves two operations in the f − k domain. First, the temporal frequency ω is mapped onto the vertical wavenumber k_{z} via equation (E-86). Second, the amplitudes are scaled by the quantity
( )
The zero-offset image is then obtained by summing over the wavenumber k_{h} (equation E-91), and inverse Fourier transforming in the midpoint y direction.
E.7 Velocity analysis by wavefield extrapolation
A method of migration velocity analysis based on wavefield extrapolation is described in migration velocity analysis. The main computational steps of this method ^{[19]} are outlined below. We work with seismic data in midpoint-(half) offset (y, h) coordinates. We want to obtain a volume of focused energy at zero offset in (y, v, τ) coordinates from a prestack data set in (y, h, t) coordinates. For a midpoint location y, migration velocity function then can be picked from the corresponding (v, τ) plane.
First, a 3-D Fourier transformation is applied to the upcoming wavefield P(y, h, τ = 0, t), which is recorded at the surface
( )
where t is the two-way traveltime and
( )
is the two-way vertical time equivalent of downward continuation depth z in a medium with velocity v(z). The variables (k_{y}, k_{h}, w) are the Fourier duals of (y, h, t).
The surface wavefield given by equation (E-93) then is extrapolated down to depth τ by
( )
where
( )
and Y and H are the normalized midpoint and offset wavenumbers given by equations (E-77a) and (E-77b). The −2 term puts the expression in retarded time form. (This term was not included in the previous definition of DSR given by equation E-76.) Equation (E-95) is used recursively to extrapolate the wavefield from one depth to another in steps of Δτ.
Next, we transform the extrapolated wave field P(k_{y}, k_{h}, τ, ω) into the space-time domain. In doing so, we only need to obtain the zero-offset information (h = 0). By summing the extrapolated wavefield over k_{h} in equation (E-95), we get the wave field at zero offset, P(k, h = 0, τ, ω). By doing the 2-D inverse transform over (k_{y}, ω), we obtain
( )
Here, P(y, h = 0, τ, t) is the zero-offset section at various depth levels from which we want to extract velocity information.
Suppose that velocity v_{e} were used to extrapolate the surface wavefield down to depth τ. Equation (E-95) is written with v_{e} and τ as
( )
Now suppose that the true medium velocity v were used to extrapolate the surface wavefield down to depth τ = t. By rewriting equation (E-95) with v and t, we have
( )
Match the two extrapolated wavefields in equations (E-98a) and (E-98b) to get a relationship between v_{e}, τ, v, and t
( )
Because of the complexity of DSR [equation (E-96)], equation (E-99) does not provide an explicit expression for v in terms of the other three variables τ, t, and v_{e}. However, we can get an approximate expression by expanding the square roots in equation (E-96) in the Taylor series and retaining the terms with Y^{2} and H^{2}, only. By using this approximate form and the definitions given by equations (E-77a) and (E-77b), the following approximate relationship results
( )
This expression suggests that downward continuation with the correct (medium) velocity to a wrong depth is equivalent to downward continuation to the correct depth with the wrong velocity ^{[20]}.
The derivation of equation (E-100) assumes that v_{e} is constant. When v_{e} is depth-variable, then the relationship in equation (E-100) still holds because equation (E-95) is valid for a stratified earth model. However, quantity v_{e} in that equation is replaced by the rms velocity.
Because the approximation made to equation (E-96) is best for small ratios of offset-to-reflector depth, the accuracy of the mapping procedure based on equation (E-100) degrades at very shallow depths. Refer to migration velocity analysis for the practical considerations of the migration velocity estimation technique described here.^{[21]} ^{[22]} ^{[23]} ^{[24]} ^{[25]} ^{[26]}
References
- ↑ Deregowski and Rocca (1981), Deregowski, S. M. and Rocca, F., 1981, Geometrical optics and wave theory for constant-offset sections in layered media: Geophys. Prosp., 29, 374–406.
- ↑ Levin, 1971, Levin, F.K., 1971, Apparent velocity from dipping interface reflections: Geophysics, 36, 510–516.
- ↑ ^{3.0} ^{3.1} ^{3.2} Liner, 1990, Liner, C. L., 1990, General theory and comparative anatomy of dip-moveout: Geophysics, 55, 595–607.
- ↑ Bolondi et al., 1982, Bolondi, G., Loinger, E. and Rocca, F., 1982, Offset continuation of seismic sections: Geophys. Prosp., 30, 813–828.
- ↑ Bale and Jacubowicz, 1987, Bale, R. and Jacubowitz, H., 1987, Poststack prestack migration: 57th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 714–717.
- ↑ ^{6.0} ^{6.1} Notfors and Godfrey, 1987, Marcoux, M. O., Godfrey, R. J., and Notfors, C. D., 1987, Migration for optimum velocity evaluation and stacking: Presented at the 49th Ann. Mtg. European Asn. Expl. Geophys.
- ↑ ^{7.0} ^{7.1} Zhou et al. (1996), Zhou, B, Mason, I. M., and Greenalgh, S.A., 1996, Accurate and efficient shot-gather dip-moveout processing in the log-stretch domain: Geophys. Prosp., 43, 963–978.
- ↑ Biondi and Ronen (1987), Biondi, B. and Ronen, J., 1987, Dip moveout in shot profiles: Geophysics, 52, 1473–1482.
- ↑ Cabrera and Levy (1989), Cabrera, J. and Levy, S., 1989, Shot dip-moveout with logarithmic transformations: Geophysics, 54, 1038–1041.
- ↑ ^{10.0} ^{10.1} ^{10.2} Bancroft et al., 1998, Bancroft, J. C., Geiger, H. D., and Margrave, G., 1998, The equivalent offset method of prestack time migration: Geophysics, 63, 2042–2053.
- ↑ Margrave et al., 1999, Margrave, G., Bancroft, J. C., and Geiger, H. D., 1999, Fourier prestack migration by equivalent wavenumber: Geophysics, 64, 197–207.
- ↑ Claerbout (1985), Claerbout, J.F., 1985, Imaging the earth’s interior: Blackwell Scientific Publications.
- ↑ Bancroft and Geiger (1994), Bancroft, J. C. and Geiger, H. D., 1994, Equivalent-offset CRP gathers: 64th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 672–675.
- ↑ Bancroft et al., 1997, Bancroft, J. C., Margrave, G., and Geiger, H. D., 1997, A kinematic comparison of DMO-PSI and equivalent offset migration (EOM): 67th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1575–1578.
- ↑ Berryhill, 1996, Berryhill, J. R., 1996, System and method of seismic shot-record migration: U. S. Patent 5,544,126.
- ↑ Ottolini (1982), Ottolini, R., 1982, Migration of seismic data in angle-midpoint coordinates: Ph.D. thesis, Stanford University.
- ↑ Fowler, 1997, Fowler, P., 1997, A comparative overview of prestack time migration methods: 67th Ann. Int. Soc. Explor. Geophys. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 1571–1574.
- ↑ Stolt, 1978, Stolt, R.H., 1978, Migration by Fourier transform: Geophysics, 43, 23–48.
- ↑ Yilmaz and Chambers, 1984, Yilmaz, O. and Chambers, R., 1984, Migration velocity analysis by wave field extrapolation: Geophysics, 49, 1664–1674.
- ↑ Doherty and Claerbout (1974), Doherty, S.M. and Claerbout, J.F., 1974, Velocity analysis based on the wave equation: Stanford Expl. Proj., Rep. No. 1, Stanford University.
- ↑ Deregowski, S. M., 1982, Dip-moveout and reflector-point dispersal: Geophys. Prosp., 30, 318–322.
- ↑ Deregowski, S. M., 1986, What is DMO?: First Break, 4, 7–24.
- ↑ Hale, D., 1991, A nonaliased integral method for dip moveout: Geophysics, 56, 795–805.
- ↑ Hale, D., Hill, N. R., and Stefani, J., 1992, Imaging salt with turning seismic waves: Geophysics, 57, 1453–1462.
- ↑ Notfors, C. D. and Godfrey, R. J., 1987, Dip-moveout in the frequency-wavenumber domain: Geophysics, 52, 1718–1721.
- ↑ Rocca, F., Bolondi, G., and Loinger, E., 1982, Offset continuation of seismic sections: Geophys. Prosp., 30, 813–828.
See also
- Introduction to dip-moveout correction and prestack migration
- Principles of dip-moveout correction
- Dip-moveout correction in practice
- Prestack time migration
- Migration velocity analysis
- Exercises