087 GC2014 Deterministic Marine Deghosting Tutorial and Advances
087 GC2014 Deterministic Marine Deghosting Tutorial and Advances
087 GC2014 Deterministic Marine Deghosting Tutorial and Advances
Introduction
It is well-known that marine seismic data resolution is degraded by the presence of sea-surface
reflections on both source and receiver sides. The slightly delayed reflections trailing the original source
pulse are called “ghosts”, and the resulting interference can be either constructive or destructive for
different wavelengths. This interference is manifest as a series of notches in the amplitude spectrum of
the ghosting operator, with rapid variation in the phase spectrum occurring in the vicinity of these
notches. Figure 1 illustrates this notching phenomenon for a synthetic source-side ghost, and the green
curve in Figure 2 shows notches on a real data example, where the presence of both source and
receiver ghosts conspire with anelastic attenuation to drastically reduce the effective frequency content
of the embedded wavelet. Because of the natural diversity provided by variation in propagation
directions and shot/receiver depths, as well as by imperfections in the sea-surface reflections, the
notches on this poststack display are not as deep as those shown in the prestack synthetic example in
Figure 1, Still, the notching effect is very significant, and requires compensation.
There are several acquisition-based schemes aimed at suppressing this ghost effect, including
variable-depth streamers or slanted cables (Soubaras and Dowle, 2010) and dual-sensor streamers
combined with random-depth sources (Tenghamn et al.,2007; Carlson et al., 2007). The present paper
deals with processing-based, rather than acquisition-based, solutions. A processing solution is
desirable for two reasons: first it does not require added acquisition effort, and second it is applicable to
existing legacy data which have not been acquired using any of the above specialized acquisition
schemes.
Figure 2: Average amplitude spectrum of a stacked data volume from West of Shetlands, UK. Green and red curves show the spectrum with
and without application of the deghosting strategy discussed in this paper, respectively.
There are two possible approaches to deghosting in processing. The first of these is a statistical
approach based on the classical minimum phase deconvolution algorithm of Robinson (1957), and the
second of these is a deterministic approach which entails computation of the analytic inverse of the
ghost operator. While statistical deconvolution has remained the workhorse in wavelet processing of
land data, it has been eschewed in modern marine processing flows in favour of deterministic
techniques. The present paper discusses a new deterministic method for deghosting.
Theory and/or Method
Deghosting basics
We begin by revisiting the canonical model for the normal incidence ghost. The ghost effect may be
described by the convolution of the uncontaminated trace with a ghosting operator. A one-sided ghost
operator may be written in the time domain as
Note that the functional form of equation (2) gives rise to the notches and phase pathologies shown in
Figure 1; in particular, notch depth increases as r approaches unity.
Because the ghost operator and its inverse are minimum phase, it is possible to deghost using
standard minimum phase deconvolution as noted earlier. In fact, such a scheme might seem
compelling since it carries the practical advantage of not needing to know r nor td. However, it also
imposes a well-known assumption on the distribution of the reflectivity (e.g., white reflectivity), and such
an assumption is viewed as being unnecessarily restrictive in marine processing where the wavelet
distortion mechanisms, being somewhat simpler than in the land case, admit deterministic
compensation. The assumption imposed on the reflectivity (i.e., in statistical deconvolution) runs an
especially large risk of violation in the case of estimating the deghosting operator, whose long time
domain length (and associated rapid variation across frequencies in its amplitude spectrum) makes it
particularly hard to unravel from the underlying geology. Consequently, the marine processing
community has shifted away from statistical deconvolution and deterministic approaches are now
favoured for deghosting.
The first of these addresses the fact that we are dealing with oblique, rather than normal, incidence
data in typical marine surveys. Massomzadeh et al. (2013) show that the time-delay associated with a
source (or receiver) ghost for a plane wave propagating at an angle from the vertical is given by
2𝑑 cos 𝜃
𝑡𝑑 = ,
𝑣
where d is water depth and v is water velocity. Such angle-dependence results in non-stationarity of the
ghost operator as it appears on a trace in the x-t domain, especially at larger offsets. This non-
stationarity in turn violates the convolutional model that forms the cornerstone of the derivation of our
deterministic deghosting operator, and would lead to unacceptable wavelet distortions in practice.
Fortunately, we can overcome this non-stationarity issue by performing the deghosting in the slant-
stack (-p), rather than (x,t), domain, By applying a shot-domain -p transform, we focus all energy
associated with a particular ray parameter p (i.e., with a particular emergence angle) onto a single
The second important modification is to impart frequency and ray-parameter dependence to the surface
reflection coefficient r. Our justification for such dependence is based on the fact that the sea-surface is
not a perfect mirror (Williams and Pollatos, 2012) and invokes the physically-rooted argument that
higher frequency and larger ray-parameter components of the wavefield ought to experience a less
perfect reflection at the sea surface than their lower frequency/ray-parameter counterparts After
recasting equation 3 in terms of ray parameter p, we may thus write our modified deterministic one-
sided deghosting operator as
1
𝐷(𝑤, 𝑝) = 𝑖𝑤2𝑑
, (4)
√1−𝑝2 𝑣 2
1 − 𝑟(𝑤, 𝑝)𝑒 𝑣
where it is understood that both source and receiver sides are included in the actual implementation.
Although the preceding analysis assumes a flat streamer cable. Masoomzadeh et al. (2013) describe
certain modifications (not discussed here) which allow the methodology to be extended to the case of a
linearly slanted cable.
Stochastic aspects
Selecting optimal values of shot/receiver depths and frequency-dependent reflection coefficients for use
in equation (4) may be difficult, especially in the presence of significant sea surface undulations brought
on by bad weather. Therefore, we typically perform a stochastic search for the most appropriate set of
deghosting parameters, wherein the quality of each trial parameter set is judged by the autocorrelations
of the provisionally deghosted traces.
Further improvements may be achieved by sorting the stochastically-deghosted traces into common
ray-parameter ensembles. For each ensemble a separate single “global” operator is computed using a
statistical deconvolution approach in which operator design is accomplished by averaging across all
traces within the common-p ensemble. For the reader who is well-versed in surface-consistent
deconvolution of land data, these p-dependent global operators are somewhat reminiscent of the “line”
component in that latter algorithm, This final deconvolution step can help address any remmant
deghosting imperfections which are common to the entire common-p ensemble, while safeguarding
against the unwitting removal of geology through the use of global, rather than trace-by-trace,
operators.
Examples
Conclusions
Resolution of marine seismic data may be significantly degraded by the presence of source and
receiver side ghosts. We have presented a technique for deghosting in the -p domain which is based
on a deterministic framework, but which contains stochastic elements for parameter optimization. The
method, which can be adapted to handle linearly slanted cables under certain modifications, has been
applied successfully on numerous data sets and is considered an integral part of our marine wavelet
processing toolkit.
Acknowledgements
We thank TGS for granting us permission to show the real data example. We also thank Anthony
Hardwick and Simon Baldock for their contributions in reviewing the abstract.
References
Carlson, D., W. Söllner, H. Tabti, E. Brox, and M. Widmaier, 2007, Increased resolution of seismic data from a dual-sensor streamer cable:
77th Annual International Meeting, SEG, Expanded Abstracts, 994-998, http://dx.doi.org/10.1190/1.2792572.
Masoomzadeh, H., N. Woodburn, and A. Hardwick, 2013,Broadband processing of linear streamer data: 83 rd Annual International
Meeting, SEG, Expanded Abstracts, 4635-4638. http://dx.doi.org/10.1190/segam2013-0872.1
Robinson, E., 1957, Predictive Decomposition of Seismic Traces: Geophysics, 22, 767-778., etc, etc.
Soubaras, R., and R. Dowle, 2010,Variable depth streamer—A broadband marine solution: First Break, 28, 89-96.
Tenghamn, R., S. Vaage, and C. Borresen, 2007, A dual-sensor, towed marine streamer: Its viable implementation and initial results: 77th
Annual International Meeting, SEG, Expanded Abstracts, 989-993. http://dx.doi.org/10.1190/1.2792571.
Williams, R.G., and J. Pollatos, 2012, Signal to noise—The key to increased marine seismic bandwidth: First Break, 30, 101-105.