next up previous print clean
Next: REFERENCES Up: Three-dimensional seismic data regularization Previous: Solving the Cauchy problem

The kinematics of offset continuation

 In this Appendix, I apply an alternative method to derive equation ([*]), which describes the summation path of the integral offset continuation operator. The method is based on the following considerations.

The summation path of an integral (stacking) operator coincides with the phase function of the impulse response of the inverse operator. Impulse response is by definition the operator reaction to an impulse in the input data. For the case of offset continuation, the input is a reflection common-offset gather. From the physical point of view, an impulse in this type of data corresponds to the special focusing reflector (elliptical isochrone) at the depth. Therefore, reflection from this reflector at a different constant offset corresponds to the impulse response of the OC operator. In other words, we can view offset continuation as the result of cascading prestack common-offset migration, which produces the elliptic surface, and common-offset modeling (inverse migration) for different offsets. This approach resemble that of Deregowski and Rocca 1981. It was also applied to a more general case of azimuth moveout (AMO) by Fomel and Biondi 1995b. The geometric approach implies that in order to find the summation pass of the OC operator, one should solve the kinematic problem of reflection from an elliptic reflector whose focuses are in the shot and receiver locations of the output seismic gather.

In order to solve this problem , let us consider an elliptic surface of the general form  
 \begin{displaymath}
h(x)=\sqrt{d^2-\beta\,(x-x')^2}\;,\end{displaymath} (307)
where $0 < \beta < 1$. In a constant velocity medium, the reflection ray path for a given source-receiver pair on the surface is controlled by the position of the reflection point x. Fermat's principle provides a required constraint for finding this position. According to Fermat's principle, the reflection ray path corresponds to a stationary value of the travel-time. Therefore, in the neighborhood of this path,  
 \begin{displaymath}
{\partial \tau(s,r,x) \over \partial x} = 0\;,\end{displaymath} (308)
where s and r stand for the source and receiver locations on the surface, and $\tau$ is the reflection traveltime  
 \begin{displaymath}
\tau(s,r,x) = { \sqrt{h^2(x)+(s-x)^2} \over v} + 
 { \sqrt{h^2(x)+(r-x)^2} \over v}\;.\end{displaymath} (309)

Substituting ([*]) and ([*]) into ([*]) leads to a quadratic algebraic equation on the reflection point parameter x. This equation has the explicit solution  
 \begin{displaymath}
x(s,r)= x' + {{\xi^2+H^2-h^2+\mbox{sign}(h^2-H^2)\,
\sqrt{\l...
 ...i^2-H^2-h^2\right)^2-4\,H^2\,h^2}\over
{2\,\xi\,(1-\beta)}}}\;,\end{displaymath} (310)
where h=(r-s)/2, $\xi = y-x'$, y=(s+r)/2, and $H^2=d^2\,\left({1
 \over \beta} - 1\right)$. Replacing x in equation ([*]) with its expression ([*]) solves the kinematic part of the problem, producing the explicit traveltime expression  
 \begin{displaymath}
\tau(s,r)=\left\{
 \begin{array}
{lcr}\displaystyle{
{1 \ove...
 ...er {1-\beta}}}
& \mbox{for} & h^2 < H^2
 \end{array} \right.\;,\end{displaymath} (311)
where
\begin{eqnarray}
f=\sqrt{(r-x')^2-H^2}\; & , & \;g=\sqrt{(s-x')^2-H^2}\;,
\nonum...
 ...F=\sqrt{H^2-(r-x')^2}\; & , & \;G=\sqrt{H^2-(s-x')^2}\;.
\nonumber\end{eqnarray}

The two branches of equation ([*]) correspond to the difference in the geometry of the reflected rays in two different situations. When a source-and-receiver pair is inside the focuses of the elliptic reflector, the midpoint y and the reflection point x are on the same side of the ellipse with respect to its small semi-axis. They are on different sides in the opposite case (Figure [*]).

 
ell
Figure 2
Reflections from an ellipse. The three pairs of reflected rays correspond to a common midpoint (at 0.1) and different offsets. The focuses of the ellipse are at 1 and -1.

ell
view burn build edit restore

If we apply the NMO correction, equation ([*]) is transformed to  
 \begin{displaymath}
\tau_n(s,r)=\left\{
 \begin{array}
{lcr}\displaystyle{
{1 \o...
 ...,h^2+(F+G)^2}}
& \mbox{for} & h^2 < H^2
 \end{array} \right.\;.\end{displaymath} (312)
Then, recalling the relationships between the parameters of the focusing ellipse r, x' and $\beta$ and the parameters of the output seismic gather Deregowski and Rocca (1981)

 
 \begin{displaymath}
r={ {v\,t_n} \over 2}\;,\;x'=y\;,\;
\beta={t_n^2 \over {t_n^2+{{4\,h^2} \over v^2}}}\;,\;
H=h\;,\end{displaymath} (313)
and substituting expressions ([*]) into equation ([*]) yields the expression  
 \begin{displaymath}
t_1(s_1,r_1;s,r,t_n)=\left\{
 \begin{array}
{lcr}\displaysty...
 ...^2+(F+G)^2}}
& \mbox{for} & h_1^2 < h^2
 \end{array} \right.\;,\end{displaymath} (314)
where
\begin{eqnarray}
f=\sqrt{(r_1-r)\,(r_1-s)}\;,\;g=\sqrt{(s_1-r)\,(s_1-s)}\;,
\non...
 ...\sqrt{(r-r_1)\,(r_1-s)}\;,\;G=\sqrt{(s_1-r)\,(s-s_1)}\;.
\nonumber\end{eqnarray}

It is easy to verify algebraically the mathematical equivalence of equation ([*]) and equation ([*]) in the main text. The kinematic approach described in this appendix applies equally well to different acquisition configurations of the input and output data. The source-receiver parameterization used in ([*]) is the actual definition for the summation path of the integral shot continuation operator Bagaini and Spagnolini (1993, 1996); Schwab (1993). A family of these summation curves is shown in Figure [*].

 
shc
shc
Figure 3
Summation paths of the integral shot continuation. The output source is at -0.5 km. The output receiver is at 0.5 km. The indexes of the curves correspond to the input source location.
view burn build edit restore

 

 


next up previous print clean
Next: REFERENCES Up: Three-dimensional seismic data regularization Previous: Solving the Cauchy problem
Stanford Exploration Project
12/28/2000