Proof of amplitude equivalence

This section discusses the connection between the laws of traveltime transformation and the laws of the corresponding amplitude transformation. The change of the wave amplitudes in the OC process is described by the first-order partial differential transport equation (6). The general solution of this equation can be found by applying the method of characteristics. It takes the explicit integral form  

$$A_n\left(t_n\right)=A_0\left(t_0\right)\,\exp{\left(\int_{t_... ...tial \tau_n \over \partial h} \right)^{-1}\, d\tau_n\right)}\;.$$ (30)

The integral in (30) is defined on a curved time ray, and A_n(t_n) stands for the amplitude transported along this ray. In the case of a plane dipping reflector, the ray amplitude can be immediately evaluated by substituting the explicit traveltime and time ray formulas from the preceding section into (30). The amplitude expression in this case takes the simple form  

$$A_n\left(t_n\right)=A_0\left(t_0\right)\,{t_0 \over t_n}\;.$$ (31)

In order to consider the more general case of a curvilinear reflector, we need to take into account a connection between the traveltime derivatives in (30) and the geometric quantities of the reflector. As follows directly from the trigonometry of the incident and reflected rays triangle (Figure 1),

$$\begin{eqnarray} h & = & {r-s \over 2}= D\,{{\cos{\alpha}\,\sin{\gamma}\,\cos{\g... ...os^2{\alpha}-\sin^2{\gamma}}}\;, \\ y_0 & = & x+D\,\sin{\alpha}\;,\end{eqnarray}$$ (32) (33) (34)

where D is the length of the normal ray. Let $\tau_0=2\,D/v$ be the zero-offset reflection traveltime. Combining (32) and (34) with (10) we can get the following relationship:  

$$a={\tau_n\over\tau_0}={{\cos{\alpha}\,\cos{\gamma}}\over \le... ...a}}}\right)^{1/2}= {h\,\over\sqrt{h^2-\left(y-y_0\right)^2}}\;,$$ (35)

which interprets the ``DMO smile'' (28) found by Deregowski and Rocca 1981 in geometric terms. Equation (35) allows a convenient change of variables in (30). Let the reflection angle $\gamma$ be a parameter monotonously increasing along a time ray. In this case, each time ray is uniquely determined by the position of the reflection point, which in turn is defined by the values of D and $\alpha$. According to this change of variables we can differentiate (35) along a time ray to get  

$${{d\tau_n}\over\tau_n}=-{{\sin^2{\alpha}}\over {2\,\cos^2{\g... ...amma}-\sin^2{\alpha}\right)}}\, d\left(\cos^2{\gamma}\right)\;.$$ (36)

Note also that the quantity $h\,\left(\tau_n\,{\partial \tau_n \over \partial h} \right)^{-1}$ in (30) coincides exactly with the time ray invariant C₃ found in (18). Therefore its value is constant along each time ray and equals  

$$h\,\left(\tau_n\,{\partial \tau_n \over \partial h}\right)^{-1}= -{v^2 \over 4\, \sin^2{\alpha}}\;.$$ (37)

Finally, as shown in Appendix B,  

$$\tau_n\,\left({\partial^2 \tau_n \over \partial y^2}- {\part... ...}\,\left({\sin^2{\alpha}+DK}\over {\cos^2{\gamma}+DK}\right)\;,$$ (38)

where K is the reflector curvature at the reflection point. Substituting (36), (37), and (38) into (30) transforms the integral to the form  

$$= -{1 \over 2}\,\int_{\cos^2{\gamma_0}}^{\cos^2{\gamma}} \le... ...ver {\cos^2{\gamma'}+DK}}\right)\,d\left(\cos^2{\gamma'}\right)$$ (39)

which we can evaluate analytically. The final formula for the amplitude transformation takes the form

$$\begin{eqnarray} A_n & = & A_0\,{\sqrt{\cos^2{\gamma}-\sin^2{\alpha}}\over \sqrt... ...eft( {\cos^2{\gamma_0}+DK}\over{\cos^2{\gamma}+DK}\right)^{1/2}\;.\end{eqnarray}$$ (40)

In case of a plane reflector, the curvature K is zero, and (40) coincides with (31). Equation (40) can be rewritten as  

$$A_n={{c\,\cos{\gamma}}\over{\tau_n\,\sqrt{\cos^2{\gamma}+DK}}}\;,$$ (41)

where c is constant along each time ray (it may vary with the reflection point location on the reflector but not with the offset). We should compare equation (41) with the known expression for the reflection wave amplitude of the leading ray series term in 2.5-D media:  

$$A={{C_R(\gamma) \Psi}\over G}\;,$$ (42)

where C_R stands for the angle-dependent reflection coefficient, G is the geometric spreading  

$$G=v \tau {\sqrt{\cos^2{\gamma}+DK}\over \cos{\gamma}}\;,$$ (43)

and $\Psi$ includes other possible factors (such as the source directivity) that we can either correct or neglect in the preliminary processing. It is evident that the curvature dependence of the amplitude transformation (41) coincides completely with the true geometric spreading factor (43), and that the angle dependence of the reflection coefficient is not provided by the offset continuation process. If the wavelet shape of the reflected wave on seismic sections (R_n in (4)) is described by the delta function, then, as follows from the known properties of this function,  

$$A\,\delta\left(t-\tau(y,h)\right)=\left\vert{{dt_n} \over {d... ...ight) = {t \over t_n}\,A\,\delta\left(t_n-\tau_n(y,h)\right)\;,$$ (44)

which leads to the equality  

$$A_n=A\,{t \over t_n}\;.$$ (45)

Combining (45) with (42) and (41) allows us to evaluate the amplitude after continuation from some initial offset h₀ to another offset h₁, as follows:  

$$A_1={{C_R(\gamma_0) \Psi_0}\over G_1}\;.$$ (46)

Equation (46) indicates that the OC process described by equation (1) is amplitude-preserving in the sense that corresponds to the so-called Born DMO Bleistein (1990); Liner (1991). This means that the geometric spreading factor from the initial amplitudes is transformed to the true geometric spreading on the continued section, while the reflection coefficient stays the same. This remarkable dynamic property allows AVO (amplitude versus offset) analysis to be performed by a dynamic comparison between true constant-offset sections and the sections transformed by OC from different offsets. With a simple trick, the offset coordinate is transferred to the reflection angles for the AVO analysis. As follows from (35) and (10),  

$${\tau_n^2 \over \tau\,\tau_0}=\cos{\gamma}\;.$$ (47)

If we include the ${t_n^2 \over t\,t_0}$ factor in the DMO operator (continuation to zero offset) and divide the result by the DMO section obtained without this factor, the resultant amplitude of the reflected events will be directly proportional to $\cos{\gamma}$, where the reflection angle $\gamma$ corresponds to the initial offset. Of course, this conclusion is rigorously valid for constant-velocity 2.5-D media only.

Black et al. 1993 recently suggested a definition of true-amplitude DMO different from that of Born DMO. The difference consists of two important components:

1.

True-amplitude DMO addresses preserving the peak amplitude of the image wavelet instead of preserving its spectral density. In the terms of this paper, the peak amplitude corresponds to the initial amplitude A instead of the spectral density amplitude A_n. A simple correction factor $t \over t_n$would help us take the difference between the two amplitudes into account. Multiplication by $t \over t_n$ can be easily done at the NMO stage.

2.

Seismic sections are multiplied by time to correct for the geometric spreading factor prior to DMO (or in our case, offset continuation) processing.

As follows from (43), multiplication by t is a valid geometric spreading correction for plane reflectors only. It is amplitude-preserving offset continuation based on the OC equation (1) that is able to correct for the curvature-dependent factor in the amplitude. To take into account the second aspect of Black's definition, we can consider the wave field $\hat{P}$ such that  

$$\hat{P}\left(y,h,t_n\right)=t\,P\left(y,h,t_n\right)\;.$$ (48)

Substituting (48) into the OC equation (1) transforms the latter to the form  

$$h \, \left( {\partial^2 \hat{P} \over \partial y^2} - {\part... ... t_n \, \partial h}}\, -{\partial \hat{P} \over \partial h}\; .$$ (49)

Equations (49) and (1) differ only with respect to the first-order term $\partial \hat{P} \over \partial h$.This term affects the amplitude behavior but not the traveltimes, since the eikonal-type equation (5) depends on the second-order terms only. Offset continuation operators based on (49) conform to Black's definition of true-amplitude processing.