Next: Datuming operator Up: POSTSTACK DATUMING Previous: POSTSTACK DATUMING

## Forward modeling operator

In forward modeling, the wavefield recorded at each geophone along an irregular surface is the wavefield propagated up to the depth level where the geophones are located. It is necessary to stop the wavefield propagation after recording because the reflection coefficient at surface is almost -1 Ji and Claerbout (1992). To do so, I formulate the forward modeling operator by propagating the wavefield upward with a filter between extrapolation steps to stop the wavefield propagation where it is recorded. Then all wavefields from all depth levels where the geophones are located are summed together to produce the wavefield along the irregular surface. In order to explain the algorithm clearly and schematically, I use a simple topography model that has only eight geophone groups on an irregular surface, as illustrated in Figure 1.

 syngeometry Figure 1 Synthetic surface recording geometry. solid squares represent geophone location on an undulating surface.

Figure 2 schematically describes the forward modeling algorithm for the simple model. Wi represents upward extrapolation at the i-th depth level and F1, F2, and F3 are spatial filters for grabbing the wavefield where the geophones are located at the corresponding depth levels. The operators I-F3 and I-F2-F3 in Figure 2 stop the wavefield at the locations where it is recorded below or at the corresponding depth level and pass the wavefield at the locations where it is not yet recorded. Each small rectangle in Figure 2 represents an abstract vector that contains wavefields at the corresponding space location. The wavefield along the irregular surface is obtained by summing the wavefields that are grabbed at the various depth levels.

tpfrdschm
Figure 2
Forward modeling scheme: the schematic diagram for forward depth extrapolation when the surface is not flat. Wi represents the upward extrapolation operator at the i-th depth level. F1, F2, and F3 are spatial filters shown in the text, and I is the identity matrix.

The forward modeling scheme shown in Figure 2 can be algebraically generalized, if we divide the topography into z levels, as follows:

 (1)

 (2)

 (3)
where

In equation (1), d0 and dz are wavefields on the irregular surface and the datum level, respectively. The extrapolation operator E is followed by the spatial filter G at every depth level. We can see that the upward extrapolation operator Wi is applied to the wavefield that does not arrived at the surface because the operator Ki-1 remove the wavefield if it has arrived at any previous depth level. All wavefields that arrive at the surface are saved by the operator Fi-1 for the final output. For the simple geometry shown in Figure 1, F1, F2, and F3 are just diagonal matrices whose elements are 1 where the geophones are located and elsewhere. Thus their diagonal elements are as follows: The operator Wi in equation (2) can be any extrapolation scheme including the Kirchhoff, phase-shift, split-step, or finite-difference method. If we use the phase-shift extrapolation algorithm for Wi, we need an additional inverse Fourier transform in every extrapolation step because the operator G is in the space domain. However, all other algorithms, such as the Kirchhoff, split-step, and finite-difference methods, don't need any additional computation except the operation by G, which is the multiplication of the extrapolated wavefield by the zero/one filter.

Next: Datuming operator Up: POSTSTACK DATUMING Previous: POSTSTACK DATUMING
Stanford Exploration Project
11/17/1997