UBC Faculty Research and Publications

Irregular sampling: from aliasing to noise Hennenfent, Gilles 2008

You don't seem to have a PDF reader installed, try download the pdf

Item Metadata


hennenfent07eage.pdf [ 470.1kB ]
JSON: 1.0107401.json
JSON-LD: 1.0107401+ld.json
RDF/XML (Pretty): 1.0107401.xml
RDF/JSON: 1.0107401+rdf.json
Turtle: 1.0107401+rdf-turtle.txt
N-Triples: 1.0107401+rdf-ntriples.txt

Full Text

Irregular sampling: from aliasing to noise Gilles Hennenfent and Felix J. Herrmann ABSTRACT Seismic data is often irregularly and/or sparsely sampled along spatial coordinates. We show that these acquisition geometries are not necessarily a source of adversity in order to accurately reconstruct adequately-sampled data. We use two examples to illustrate that it may actually be better than equivalent regularly subsampled data. This comment was already made in earlier works by other authors. We explain this behavior by two key observations. Firstly, a noise-free underdetermined problem can be seen as a noisy well-determined problem. Secondly, regularly subsampling creates strong coherent acquisition noise (aliasing) difficult to remove unlike the noise created by irregularly subsampling that is typically weaker and Gaussian-like. Adequate sampling of seismic data is traditionally understood as evenly-distributed time and space measurements of the reflected wavefield. Moreover, the sampling rate along each axis must be equal to or above twice the highest frequency/wavenumber of the contin- uous signal being discretized (Shannon/Nyquist sampling theorem). In practice, however, seismic data is often irregularly and/or sparsely sampled along spatial coordinates, which is generally considered as a nuisance since it breaks one or both previously-stated conditions of adequate sampling. It turns out that these acquisition geometries are not necessarily a source of adversity to accurately reconstruct adequately-sampled data. This new theory, developed in the information theory community, is referred to in the literature by the terms “compressed sensing” or “compressive sampling” (see e.g. Donoho, 2006, and the com- pressed sensing resources of the Rice University DSP group at www.dsp.ece.rice.edu/cs for more references). Undersampled seismic data: Fig. 1 depicts two sub-Nyquist spatial samplings (time is adequately sampled) of some adequately-sampled synthetic seismic data (Fig. 1(a)) and their effect in the f -k domain. When the data of Fig. 1(a) is regularly subsampled by a fac- tor of two (Fig. 1(b)), the well-known phenomenon of aliasing occurs. The spectrum of the adequately-sampled data (Fig. 1(d)) is corrupted by strong aliases beyond 30 Hz (Fig. 1(e)) leaving the low frequencies unaltered. In the case of irregular subsampling by a factor of two (Fig. 1(c)), the spectrum is globally altered by some weak noise (Fig. 1(f)). From an interpolation perspective, the data of Fig. 1(b) is typically perceived as easier to handle than the data of Fig. 1(c) and a variety of algorithms are available (see e.g. Gulu- nay, 2003; Abma and Kabir, 2005, for an overview). The key idea is to use the non-aliased lower frequencies to de-alias higher frequencies (Spitz, 1991; Gulunay, 2003), which is not immediately possible in the case of irregularly sampled data since the whole spectrum is degraded. Besides, as rightly pointed out by Abma and Kabir (2005), most interpolation 2algorithms assume regular spatial sampling (e.g. f -k, t-x, and f -x prediction filter algo- rithms). Inspired by Donoho et al. (2006), we propose to look at the interpolation problem as a denoising problem or, in other words, to turn a noise-free underdetermined problem into a noisy well-determined problem (Donoho’s slogan). Indeed, Figs. 1(e) and 1(f) are noisy versions of Fig. 1(d). In the case of regular subsampling, the noise is highly structured and has strong amplitudes, which makes it difficult to remove. Conversely, irregular subsam- pling generates lower amplitude noise that looks Gaussian, which is easier to remove. It thus seems that irregular subsampling may actually help to recover more accurately seismic data in a much wider frequency band than would be possible using equivalent regular sub- sampling or any subsampling that generates structured noise. This comment, already made by other authors in earlier works, is illustrated in the next section using the Curvelet Re- covery with Sparsity-Promoting Inversion (CRSI) algorithm (Herrmann and Hennenfent, 2007). (a) (b) (c) (d) (e) (f) Figure 1: Seismic data with different spatial samplings (a,b,c) and their corresponding amplitude spectrum (d,e,f). Reconstruction using CRSI: Among the many applications of multiscale (and multidi- rectional) transforms (e.g. wavelets, seislets, curvelets, and other *lets) denoising is prob- ably the most popular. The crucial point is to find a transform which can represent the noise-free signal with few atoms (this argument is called sparsity) and at the same time poorly represents the contaminating noise. In such a transform domain the signal’s energy is concentrated in few significant coefficients whereas the noise’s energy is spread across the transform domain. Hence, it is easy to discriminate between noise and signal and to re- cover the noise-free signal (see e.g. Mallat, 1998, for more details). In fact, the sparser the representation of the noise-free signal in the transform domain, the easier the separation. CRSI uses curvelets, which arguably provide the sparsest non-adaptive representation for 3seismic data (see e.g. Herrmann and Hennenfent, 2007, and other curvelet-related contri- butions to the proceedings of this meeting), to denoise (read: interpolate) seismic data with missing traces and thus separate the seismic signal from the acquisition noise. Synthetic data example: We consider two different sparse acquisition geometries for a 2D survey (one 3D shot with these acquisition patterns would lead to similar results). Adequately-sampled data is reconstructed using 3D CRSI on the prestack volume. No minimum velocity constraint and no frequency restriction are used in CRSI to better sin- gle out the effect of the acquisition geometry on the reconstruction. Fig. 2(a) shows the prestack volume where every shot was recorded by a different subset of receivers (20% of the total number of receivers), Fig. 2(c) the reconstructed data (signal-to-noise ratio (SNR) = 16.92 dB), and Fig. 2(e) the difference cube. In comparison, Fig. 2(b) shows the prestack volume where the same subset of receivers (also 20% of the total number of receivers) recorded all the shots, Fig. 2(d) the reconstructed data (SNR = 9.26 dB), and Fig. 2(f) the difference cube. This experiment clearly illustrates our comment that irregular sampling favors in fact the recovery of seismic data over sampling that exhibits structure for the same amount of data collected. Real data example: We examine the Oseberg 3D walkaway VSP survey to further illus- trate the benefit of irregular sampling. The data (Fig. 3(a)) only contains 10% of the traces we are interested in. The irregular sampling geometry comes from the fact that the source vessel sailed in circles. The reconstructed data using 3D CRSI is presented in Fig. 3(b). It turns out that, even in this vastly underdetermined problem, the acquisition noise gen- erated by the irregular sampling geometry can be removed and the data fairly accurately recovered. Conclusions: We proposed to look at the seismic data interpolation problem from a de- noising perspective. From this standpoint, we noted that regular subsampling geometries generate acquisition noise more difficult to remove than irregular subsampling geometries do for the same amount of data collected. It thus seems that irregular subsampling may ac- tually help to recover seismic data in a much wider frequency band than would be possible using equivalent regular subsampling or any subsampling that generates structured acquisi- tion noise. We believe this new insight may lead to new acquisition strategies. On land, for example, a regular sampling may lead to aliased ground-roll that needs to be interpolated to a finer grid in order to be removed. Our observations suggest one should irregularly sample on the finer grid instead. This could lead to a better ground-roll interpolation and hence ground-roll removal. Acknowledgments: This work was in part financially supported by the NSERC Discov- ery Grant 22R81254 and CRD Grant DNOISE 334810-05 of F.J. Herrmann and was carried out as part of the SINBAD project with support, secured through ITF, from the following organizations: BG Group, BP, Chevron, ExxonMobil and Shell. The authors thank Norsk Hydro for the real dataset. References Abma, R. and N. Kabir, 2005, Comparisons of interpolation methods: The Leading Edge, 24, 984–989. 4(a) (b) (c) (d) (e) (f) Figure 2: 2D prestack data volumes with different acquisition geometries (a,b), their re- construction using 3D CRSI (c,d), and the corresponding difference cubes (e,f). 5(a) (b) Figure 3: Oseberg 3D walkaway VSP survey and its reconstruction using 3D CRSI. Donoho, D., 2006, Compressed sensing: IEEE Trans. on Information Theory, 52, 1289– 1306. Donoho, D., Y. Tsaig, I. Drori, and J.-L. Starck, 2006, Sparse solution of underdetermined linear equations by stagewise orthogonal matching pursuit. (Preprint). Gulunay, N., 2003, Seismic trace interpolation in the Fourier domain: Geophysics, 68, 355–369. Herrmann, F. and G. Hennenfent, 2007, Non-parametric seismic data recovery with curvelet frames. (submitted). Mallat, S., 1998, A wavelet tour of signal processing: Academic Pr. Spitz, S., 1991, Seismic trace interpolation in the f -x domain: Geophysics, 56, 785–794.


Citation Scheme:


Usage Statistics

Country Views Downloads
United States 3 0
Germany 1 0
Canada 1 0
City Views Downloads
Redmond 1 0
Vancouver 1 0
Unknown 1 0
Sunnyvale 1 0
Ashburn 1 0

{[{ mDataHeader[type] }]} {[{ month[type] }]} {[{ tData[type] }]}


Share to:


Related Items