next up previous [pdf]

Next: Conclusions Up: Seismic reflection data interpolation Previous: Constant-velocity synthetic

Shot continuation

Missing or under-sampled shot records are a common example of data irregularity (Crawley, 2000). The offset continuation approach can be easily modified to work in the shot record domain. With the change of variables $s = y - h$, where $s$ is the shot location, the frequency-domain equation (1) transforms to the equation

\begin{displaymath}
h \, \left( 2\,{\partial^2 \tilde{P} \over \partial s \part...
...} -
{\partial \tilde{P} \over {\partial s}}\right) = 0 \;.
\end{displaymath} (11)

Unlike equation (1), which is second-order in the propagation variable $h$, equation (11) contains only first-order derivatives in $s$. We can formally write its solution for the initial conditions at $s=s_1$ in the form of a phase-shift operator:
\begin{displaymath}
\widehat{\widehat{P}}(s_2) = \widehat{\widehat{P}}(s_1)\,
...
...right)\,
\frac{k_h\,h-\Omega}{2\,k_h\,h-\Omega}\right]}\;,
\end{displaymath} (12)

where the wavenumber $k_h$ corresponds to the half-offset $h$. Equation (12) is in the mixed offset-wavenumber domain and, therefore, not directly applicable in practice. However, we can use it as an intermediate step in designing a finite-difference shot continuation filter. Analogously to the cases of plane-wave destruction and offset continuation, shot continuation leads us to the rational filter
\begin{displaymath}
\hat{P}_{s+1}(Z_h) =
\hat{P}_{s} (Z_h) \frac{S(Z_h)}{\bar{S}(1/Z_h)}\;,
\end{displaymath} (13)

The filter is non-stationary, because the coefficients of $S(Z_h)$ depend on the half-offset $h$. We can find them by the Taylor expansion of the phase-shift equation (12) around zero wavenumber $k_h$. For the case of the half-offset sampling equal to the shot sampling, the simplest three-point filter is constructed with three terms of the Taylor expansion. It takes the form
\begin{displaymath}
S(Z_h) = - \left(\frac{1}{12} + i\,\frac{h}{2\,\Omega}\ri...
... i\,\frac{\Omega^2 + 18\,h^2}{12\,\Omega\,h}\right)\,
Z_h\;.
\end{displaymath} (14)

Let us consider the problem of doubling the shot density. If we use two neighboring shot records to find the missing record between them, the problem reduces to the least-squares system

\begin{displaymath}
\left[\begin{array}{c}
\mathbf{S} \\
\mathbf{\bar{S}}
...
...s-1} \\
\mathbf{S}\,\mathbf{p}_{s+1}
\end{array}\right]\;,
\end{displaymath} (15)

where $\mathbf{S}$ denotes convolution with the numerator of equation (13), $\mathbf{\bar{S}}$ denotes convolution with the corresponding denominator, $\mathbf{p}_{s-1}$ and $\mathbf{p}_{s+1}$ represent the known shot gathers, and $\mathbf{p}_s$ represents the gather that we want to estimate. The least-squares solution of system (15) takes the form
\begin{displaymath}
\mathbf{p}_s = \left(
\mathbf{S}^T\,\mathbf{S} +
\mathb...
...+
\mathbf{\bar{S}}^T\,\mathbf{S}\,\mathbf{p}_{s+1}\right)\;.
\end{displaymath} (16)

If we choose the three-point filter (14) to construct the operators $\mathbf{S}$ and $\mathbf{\bar{S}}$, then the inverted matrix in equation (16) will have five non-zero diagonals. It can be efficiently inverted with a direct banded matrix solver using the $LDL^T$ decomposition (Golub and Van Loan, 1996). Since the matrix does not depend on the shot location, we can perform the decomposition once for every frequency so that only a triangular matrix inversion will be needed for interpolating each new shot. This leads to an extremely efficient algorithm for interpolating intermediate shot records.

Sometimes, two neighboring shot gathers do not fully constrain the intermediate shot. In order to add an additional constraint, I include a regularization term in equation (16), as follows:

\begin{displaymath}
\mathbf{p}_s = \left(
\mathbf{S}^T\,\mathbf{S} +
\mathb...
...+
\mathbf{\bar{S}}^T\,\mathbf{S}\,\mathbf{p}_{s+1}\right)\;,
\end{displaymath} (17)

where $\mathbf{A}$ represents convolution with a three-point prediction-error filter (PEF), and $\epsilon$ is a scaling coefficient. The appropriate PEF can be estimated from $\mathbf{p}_{s-1}$ and $\mathbf{p}_{s+1}$ using Burg's algorithm (Burg, 1975; Claerbout, 1976; Burg, 1972). A three-point filter does not break the five-diagonal structure of the inverted matrix. The PEF regularization attempts to preserve offset dip spectrum in the under-constrained parts of the estimated shot gather.

Figure 9 shows the result of a shot interpolation experiment using the constant-velocity synthetic from Figure 4. In this experiment, I removed one of the shot gathers from the original NMO-corrected data and interpolated it back using equation (17). Subtracting the true shot gather from the reconstructed one shows a very insignificant error, which is further reduced by using the PEF regularization (right plots in Figure 9). The two neighboring shot gathers used in this experiment are shown in the top plots of Figure 8. For comparison, the bottom plots in Figure 8 show the simple average of the two shot gathers and its corresponding prediction error. As expected, the error is significantly larger than the error of shot continuation. An interpolation scheme based on local dips in the shot direction would probably achieve a better result, but it is significantly more expensive than the shot continuation scheme introduced above.

shot3
shot3
Figure 8.
Top: Two synthetic shot gathers used for the shot interpolation experiment. An NMO correction has been applied. Bottom: simple average of the two shot gathers (left) and its prediction error (right).
[pdf] [png] [scons]

shotin
shotin
Figure 9.
Synthetic shot interpolation results. Left: interpolated shot gathers. Right: prediction errors (the differences between interpolated and true shot gathers), plotted on the same scale.
[pdf] [png] [scons]

A similar experiment with real data from a North Sea marine dataset is reported in Figure 11. I removed and reconstructed a shot gather from the two neighboring gathers shown in Figure 10. The lower parts of the gathers are complicated by salt dome reflections and diffractions with conflicting dips. The simple average of the two input shot gathers (bottom plots in Figure 11) works reasonably well for nearly flat reflection events but fails to predict the position of the back-scattered diffractions events. The shot continuation method works well for both types of events (top plots in Figure 11). There is some small and random residual error, possibly caused by local amplitude variations.

elfshot3
elfshot3
Figure 10.
Two real marine shot gathers used for the shot interpolation experiment. An NMO correction has been applied.
[pdf] [png] [scons]

elfshotin
elfshotin
Figure 11.
Real-data shot interpolation results. Top: interpolated shot gather (left) and its prediction error (right). Bottom: simple average of the two input shot gathers (left) and its prediction error (right).
[pdf] [png] [scons]

Analogously to the case of offset continuation, it is possible to extend the shot continuation method to three dimensions. A simple modification of the proposed technique would also allow us to use more than two shot gathers in the input or to extrapolate missing shot gathers at the end of survey lines.


next up previous [pdf]

Next: Conclusions Up: Seismic reflection data interpolation Previous: Constant-velocity synthetic

2014-03-27