Arbitrary source and receiver positioning in finite‐difference schemes using Kaiser windowed sinc functions

In finite‐difference methods a seismic source can be implemented using either initial wavefield values or body forces. However, body forces can only be specified at finite‐difference nodes, and, if using initial values, a source cannot be located close to a reflecting boundary or interface in the model. Hence, difficulties can exist with these schemes when the region surrounding a source is heterogeneous or when a source either is positioned between nodes or is arbitrarily close to a free surface.A completely general solution to these problems can be obtained by using Kaiser windowed sinc functions to define a small region around the true source location that contains several nodal body forces. Both monopole and dipole point sources can be defined, enabling many source types to be implemented in either acoustic or elastic media. Such a function can also be used to arbitrarily locate receivers. If the number of finite‐difference nodes per wavelength is four or more (and with a source region half‐width of o...

[1]  Z. Alterman,et al.  Source of Finite Extent, Applied Force and Couple in an Elastic Half-Space , 1970 .

[2]  Henrik Schmidt,et al.  Finite-difference and frequency-wavenumber modeling of seismic monopole sources and receivers in fluid-filled boreholes , 1994 .

[3]  Z. Alterman,et al.  Propagation of elastic waves in layered media by finite difference methods , 1968 .

[4]  Carlos A. Cunha.  Elastic modeling in discontinuous media , 1993 .

[5]  Olav Holberg,et al.  COMPUTATIONAL ASPECTS OF THE CHOICE OF OPERATOR AND SAMPLING INTERVAL FOR NUMERICAL DIFFERENTIATION IN LARGE-SCALE SIMULATION OF WAVE PHENOMENA* , 1987 .

[6]  James B. Cole.  A nearly exact second-order finite difference time-domain wave propagation algorithm on a coarse grid , 1994 .

[7]  R. M. Alford,et al.  ACCURACY OF FINITE‐DIFFERENCE MODELING OF THE ACOUSTIC WAVE EQUATION , 1974 .

[8]  Jacob Aboudi,et al.  Numerical Simulation of Seismic Sources , 1971 .

[9]  I. Štekl,et al.  Accurate viscoelastic modeling by frequency‐domain finite differences using rotated operators , 1998 .

[10]  R. Pratt,et al.  INVERSE THEORY APPLIED TO MULTI‐SOURCE CROSS‐HOLE TOMOGRAPHY.: PART 1: ACOUSTIC WAVE‐EQUATION METHOD1 , 1990 .

[11]  Rune Mittet,et al.  General Source And Receiver Positions In Coarse-grid Finite-difference Schemes , 1999 .

[12]  Johan O. A. Robertsson,et al.  A hybrid wave propagation simulation technique for ocean acoustic problems , 1996 .

[13]  M. Landrø,et al.  Implementing measured source signatures in a coarse-grid, finite-difference modeling scheme , 1993 .

[14]  C. Shin,et al.  An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator , 1996 .

[15]  Alan V. Oppenheim,et al.  Discrete-Time Signal Pro-cessing , 1989 .

[16]  G. McMechan,et al.  Line sources for seismic modeling by finite differences in inhomogeneous media , 1987 .

[17]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .