Amplitude preservation for offset continuation: Confirmation for Kirchhoff data

THE KIRCHHOFF MODELING APPROXIMATION

In this section, we introduce the Kirchhoff approximate integral representation of the upward propagating response to a single reflector, with separated source and receiver points. We then show how the amplitude of this integrand is related to the zero-offset amplitude at the source receiver point on the ray, making equal angles at the scattering point with the rays from the separated source and receiver. The Kirchhoff integral representation (Bleistein, 1984; Haddon and Buchen, 1981) describes the wavefield scattered from a single reflector. This representation is applicable in situations where the high-frequency assumption is valid (the wavelength is smaller than the characteristic dimensions of the model) and corresponds in accuracy to the WKBJ approximation for reflected waves, including phase shifts through buried foci. The general form of the Kirchhoff modeling integral is

$$U_S({\bf r,s, \omega}) = \int_\Sigma R ({\bf x;r,s}) {\pa... ...{\bf s, x ,\omega}) G({\bf x,r,\omega})\right] d\Sigma \; ,$$ (1)

where ${\bf s} = (s,0,0)$ and ${\bf r} = (r,0,0)$ stand for the source and the receiver location vectors at the surface of observation; ${\bf x}$ denotes a point on the reflector surface $\Sigma$; $R$ is the reflection coefficient at $\Sigma$; $n$ is the upward normal to the reflector at the point ${\bf x}$; and $U_I$ and $G$ are the incident wavefield and Green's function, respectively represented by their WKBJ approximation,

$$U_I({\bf s,x,\omega}) = F(\omega) A_s({\bf s,x}) e^{i\omega \tau_s({\bf s,x})}\;,$$ (2)

$$G({\bf x,r,\omega}) = A_r({\bf x,r}) e^{i\omega \tau_r({\bf x,r})}\;.$$ (3)

In this equation, $\tau_s({\bf s,x})$ and $A_s({\bf s,x})$ are the traveltime and the amplitude of the wave propagating from ${\bf s}$ to ${\bf x}$; $\tau_r({\bf x,r})$ and $A_r({\bf x,r})$ are the corresponding quantities for the wave propagating from ${\bf x}$ to ${\bf r}$; and $F(\omega)$ is the spectrum of the input signal, assumed to be the transform of a bandlimited impulsive source. In the time domain, the Kirchhoff modeling integral transforms to

$$u_S({\bf r,s},t) = \int_\Sigma R ({\bf x;r,s}) {\partial ... ...u_s({\bf s,x}) - \tau_r({\bf x,r})\right)\right] d{\bf x}\;,$$ (4)

with $f$ denoting the inverse temporal transform of $F$. The reflection traveltime $\tau_{sr}$ corresponds physically to the diffraction from a point diffractor located at the point ${\bf x}$ on the surface $\Sigma$, and the amplitudes $A_s$ and $A_r$ are point diffractor amplitudes, as well.

The main goal of this paper is to test the compliance of representation (4) with the offset continuation differential equation. The OC equation contains the derivatives of the wavefield with respect to the parameters of observation (${\bf s, r}$, and $t$). According to the rules of classic calculus, these derivatives can be taken under the sign of integration in formula (4). Furthermore, since we do not assume that the true-amplitude OC operator affects the reflection coefficient $R$, the offset-dependence of this coefficient is outside the scope of consideration. Therefore, the only term to be considered as a trial solution to the OC equation is the kernel of the Kirchhoff integral, which is contained in the square brackets in equations (1) and (4) and has the form

$$k({\bf s,r,x},t) = A_{sr}({\bf s,r,x}) f\left(t - \tau_{sr}({\bf s,r,x})\right)\;,$$ (5)

where

$$\tau_{sr}({\bf s,r,x}) = \tau_s({\bf s,x}) + \tau_r({\bf x,r})\;,$$ (6)

$$A_{sr}({\bf s,r,x}) = A_s({\bf s,x}) A_r({\bf x,r})\;.$$ (7)

In a 3-D medium with a constant velocity $v$, the traveltimes and amplitudes have the simple explicit expressions

$$\tau_s({\bf s,x}) = {{\rho_s({\bf s,x})} \over v}\;,\;\; \tau_r({\bf x,r}) = {{\rho_r({\bf x,r})} \over v}\;,$$ (8)

$$A_s({\bf s,x}) = {1 \over {4 \pi \rho_s({\bf s,x})}}\;,\;\; A_r({\bf x,r}) = {1 \over {4 \pi \rho_r({\bf x,r})}}\;,$$ (9)

where $\rho_s$ and $\rho_r$ are the lengths of the incident and reflected rays, respectively (Figure 1). If the reflector surface $\Sigma$ is explicitly defined by some function $z=z({\bf x})$, then

$$\rho_s({\bf s,x}) = \sqrt{({\bf x-s})^2 + z^2({\bf x})}\;,\;\; \rho_r({\bf x,r}) = \sqrt{({\bf r-x})^2 + z^2({\bf x})}\;.$$ (10)

cwpgen
Figure 1. Geometry of diffraction in a constant velocity medium: view in the reflection plane.

We then introduce a particular zero-offset amplitude, namely the amplitude along the zero offset ray that bisects the angle between the incident and reflected ray in this plane, as shown in Figure 1. We denote the square of this amplitude as $A_0$. That is,

$$A_0 = {1 \over {(4 \pi \rho_0 )^2}}\;.$$ (11)

$A_0$ is the amplitude factor that appears in the Kirchhoff integral set up for a zero-offset reflection along the ray $\rho_0$. It is, thus, the desired output factor inside the Kirchhoff integral after DMO. As follows from formulas (7) and (9), the amplitude transformation in DMO (continuation to zero offset) is characterized by the dimensionless ratio

$${A_{sr} \over A_0} = {{\rho_0^2}\over {\rho_s \rho_r}}\;,$$ (12)

where $\rho_0$ is the length of the zero-offset ray (Figure 1).

As follows from the simple trigonometry of the triangles, formed by the incident and reflected rays (the law of cosines),

$$\sqrt{\rho_s^2 + \rho_0^2 - 2 \rho_s \rho_0 \cos{\gamma}} + \sqrt{\rho_r^2 + \rho_0^2 - 2 \rho_r \rho_0 \cos{\gamma}} =$$

$$= \sqrt{\rho_s^2 + \rho_r^2 - 2 \rho_s \rho_r \cos{2\gamma}} \;,$$ (13)

where $\gamma$ is the reflection angle, as shown in the figure. After straightforward algebraic transformations of equation (13), we arrive at the explicit relationship between the ray lengths:

$${{(\rho_s + \rho_r) \rho_0} \over {2 \rho_s \rho_r}} = \cos{\gamma}\;.$$ (14)

Substituting (14) into (12) yields

$${A_{sr} \over A_0} = {\tau_0 \over \tau_{sr}} \cos{\gamma}\;,$$ (15)

where $\tau_0$ is the zero-offset two-way traveltime ( $\tau_0 = 2 \rho_0/v$).

What we have done is rewrite the finite-offset amplitude in the Kirchhoff integral in terms of a particular zero-offset amplitude. That zero-offset amplitude would arise as the geometric spreading effect if there were a reflector whose dip was such that the finite-offset pair would be specular at the scattering point. Of course, the zero-offset ray would also be specular in this case.

Amplitude preservation for offset continuation: Confirmation for Kirchhoff data