![]() |
![]() |
![]() |
![]() | Downward continuation | ![]() |
![]() |
This pseudo code Fourier transforms a wavefield
observed at the earth's surface ,
and then it marches that wavefield down into the earth (
)
filling up a three-dimensional function,
.
Then it selects
, the time of the exploding reflectors
by summing over all frequencies
.
(Mathematically,
this is like
finding the signal at
by summing over all
).
Turning from pseudocode to real code,
an important practical reality
is that computer memories are not big enough
for the three-dimensional function
.
But it is easy to intertwine the downward continuation
with the summation over
so a three-dimensional function need not be kept in memory.
Conjugate migration (modeling) proceeds in much the same way.
Beginning from an upcoming wave that is zero at great depth,
the wave is marched upward in steps
by multiplication with
.
As each level in the earth is passed,
exploding reflectors from that level are added into the upcoming wave.
Pseudo code for modeling the upcoming wave
is
The positive sign in the complex exponential is a combination of two negatives,
the up
coming wave and the upward extrapolation.
In principle,
the three loops on ,
, and
are interchangeable,
however, since this tutorial program uses
a velocity
that is a constant function of depth,
I speeded it by a large factor by putting the
-loop on the inside
and pulling the complex exponential out of the inner loop.
![]() |
![]() |
![]() |
![]() | Downward continuation | ![]() |
![]() |