Irregular and Sparse Sampling in Exploration Seismology

The aim of exploration seismology is to obtain an image of the subsurface by probing it with seismic waves. These waves are generated at the surface, by using a suitably powerful source, at various locations. The waves propagate downwards through the subsurface, and are reflected at interfaces between geological layers. They subsequently propagate upwards to the surface, where they are detected and recorded. The seismic response depends of course on the source and the receiver location. In practice, the response for one source location is detected by hundreds of receivers, typically arranged along several receiver lines. In common use today is three-dimensional (3-D) surface seismology, where a volume image is obtained, rather than a cross-section image as in 2-D seismology. The sources and receivers can, for the sake of simplicity, be assumed to be located in a horizontal plane, at the surface. The recorded seismic wavefield then is a function of five variables, p(x s , y s , x r , y r , t), with two horizontal coordinates, x s and y s , for the source position, two horizontal positions, x r and y r , for the receiver positions and the variable time, denoted by t. In 2-D seismology the signal P x s , x r , t) is a function of two spatial coordinates and the time coordinate. For readers, not familiar with exploration seismology, we will give a brief overview of some notions, required to understand this chapter, in the next few subsections.

[1]  B. Torrésani,et al.  Wavelets: Mathematics and Applications , 1994 .

[2]  T. Brubaker,et al.  Nonlinear Parameter Estimation , 1979 .

[3]  A. Duijndam,et al.  Nonuniform fast Fourier transform , 1999 .

[4]  Jon F. Claerbout,et al.  Velocity-stack and slant-stack stochastic inversion , 1985 .

[5]  Farrokh Marvasti,et al.  A Unified Approach to Zero-Crossings and Nonuniform Sampling of Single and Multidimensional Signals and Systems , 1987 .

[6]  A. Duijndam BAYESIAN ESTIMATION IN SEISMIC INVERSION. PART II: UNCERTAINTY ANALYSIS1 , 1988 .

[7]  A. Duijndam BAYESIAN ESTIMATION IN SEISMIC INVERSION. PART I: PRINCIPLES1 , 1988 .

[8]  David Middleton,et al.  Sampling and Reconstruction of Wave-Number-Limited Functions in N-Dimensional Euclidean Spaces , 1962, Inf. Control..

[9]  Dan Hampson,et al.  Inverse Velocity Stacking For Multiple Elimination , 1986 .

[10]  William H. Press,et al.  Numerical recipes , 1990 .

[11]  C. O. H. Hindriks,et al.  Reconstruction of band‐limited signals, irregularly sampled along one spatial direction , 1999 .

[12]  Sanjit K. Mitra,et al.  The Nonuniform Discrete Fourier Transform and Its Applications in Signal Processing , 1998 .

[13]  J. O. Hornkohl,et al.  A new approach to Fourier analysis of randomly sampled data using linear regression technique , 1981, ICASSP.

[14]  A. J. W. Duijndam,et al.  Reconstruction of 3-D seismic signals irregularly sampled along two spatial coordinates , 2000 .

[15]  Mauricio D. Sacchi,et al.  Estimation of the discrete Fourier transform, a linear inversion approach , 1996 .

[16]  T. Strohmer,et al.  Efficient numerical methods in non-uniform sampling theory , 1995 .

[17]  Josef Stoer,et al.  Numerische Mathematik 1 , 1989 .