next up previous print clean
Next: Linking Snell waves to Up: DOWNWARD CONTINUATION Previous: Continuation of a dipping

Downward continuation with Fourier transform

One of the main ideas in Fourier analysis is that an impulse function (a delta function) can be constructed by the superposition of sinusoids (or complex exponentials). In the study of time series this construction is used for the impulse response of a filter. In the study of functions of space, it is used to make a physical point source that can manufacture the downgoing waves that initialize the reflection seismic experiment. Likewise observed upcoming waves can be Fourier transformed over t and x.

Recall in chapter [*], a plane wave carrying an arbitrary waveform, specified by equation ([*]). Specializing the arbitrary function to be the real part of the function $ \exp [ - i \omega (t-t_0 ) ]$ gives  
 \begin{displaymath}
\hbox{moving cosine wave} \ \ =\ \ \cos
\left[ \ \omega \lef...
 ...a \ +\ 
{z \over v }\ \cos \, \theta \ -\ t \ \right) \ \right]\end{displaymath} (8)

Using Fourier integrals on time functions we encounter the Fourier kernel $\exp(-i\omega t)$.To use Fourier integrals on the space-axis x the spatial angular frequency must be defined. Since we will ultimately encounter many space axes (three for shot, three for geophone, also the midpoint and offset), the convention will be to use a subscript on the letter k to denote the axis being Fourier transformed. So kx is the angular spatial frequency on the x-axis and $ \exp ( i k_x x ) $ is its Fourier kernel. For each axis and Fourier kernel there is the question of the sign before the i. The sign convention used here is the one used in most physics books, namely, the one that agrees with equation (8). Reasons for the choice are given in chapter [*]. With this convention, a wave moves in the positive direction along the space axes. Thus the Fourier kernel for (x , z , t)-space will be taken to be  
 \begin{displaymath}
\hbox{Fourier kernel} \ =\ 
e^{ i \, k_x x} \ e^{ i \, k_z z...
 ...\omega t} 
\ = \ 
\exp [ i ( k_x x \ +\ k_z z \ -\ \omega t ) ]\end{displaymath} (9)

Now for the whistles, bells, and trumpets. Equating (8) to the real part of (9), physical angles and velocity are related to Fourier components. The Fourier kernel has the form of a plane wave. These relations should be memorized!  
 \begin{displaymath}
\vbox{\offinterlineskip
 
 \hrule
 \halign {&\vrule ... (10)

A point in $(\omega,k_x,k_z)$-space is a plane wave. The one-dimensional Fourier kernel extracts frequencies. The multi-dimensional Fourier kernel extracts (monochromatic) plane waves.

Equally important is what comes next. Insert the angle definitions into the familiar relation $ \sin^2 \theta + \cos^2 \theta = 1 $.This gives a most important relationship:  
 \begin{displaymath}
k_x^2 \ +\ k_z^2 \ \ =\ \ { \omega^2 \over v^2 }\end{displaymath} (11)
The importance of (11) is that it enables us to make the distinction between an arbitrary function and a chaotic function that actually is a wavefield. Imagine any function u(t,x,z). Fourier transform it to $U(\omega , k_x , k_z ) $.Look in the $(\omega,k_x,k_z)$-volume for any nonvanishing values of U. You will have a wavefield if and only if all nonvanishing U have coordinates that satisfy (11). Even better, in practice the (t,x)-dependence at z=0 is usually known, but the z-dependence is not. Then the z-dependence is found by assuming U is a wavefield, so the z-dependence is inferred from (11).

Equation (11) also achieves fame as the ``dispersion relation of the scalar wave equation,'' a topic developed more fully in IEI.


next up previous print clean
Next: Linking Snell waves to Up: DOWNWARD CONTINUATION Previous: Continuation of a dipping
Stanford Exploration Project
12/26/2000