Wavefield extrapolation in pseudodepth domain |

The concept of *vertical time* has a long history in seismic exploration.
Vertical time
is the vertical axis for time migration (Yilmaz, 2001; Claerbout, 1985).
It is defined as the *two-way* traveltime measured by coinciding source and receiver on the surface

For conventional time image applications, this equation is used in the context of laterally invariant media. This is however, not adequate for imaging complex strudtures. We argue that vertical time, although does not correspond to the actual two-way traveltime in complex velocites, still works fine as representation of the vertical axis. Since wavefield extrapolation is usually not zero-offset, in this work we will use the

Note that velocity does not have to be identical to the velocity in which we propagate wavefields. For example, one can choose a constant , then the vertical time is simply a scaled version of depth, . In practice, we suggest using a smooth background velocity as , because it helps regularize coordinate grids. Here we retain the name ``vertical time'' for , however, it should not be confused with the vertical time used in time processing. In other words, Equation 2 represent only a change of variable from to .

Normally, wavefields are discretized on the Cartesian mesh with equally-spaced grids. For a monochromatic wave, wavelength changes with velocity and since the grid spacing is held constant, the number of samples per wavelength increases in layers with high velocities and decreases in layers with low velocities. To avoid spatial aliasing, the maximum grid spacing is limited by , where is the lowest velocity, often located in shallow layers. As a result, the deep layers with high velocities are often oversampled. The increased sampling of the layers with high velocity raises the cost of wavefield extrapolation without enhancing image resolution. Introducing vertical time partially resolves this problem. This can be seen by taking difference of equation 2 between to time levels,

(3) |

which corresponds to a fixed sampling, . This implies that the effective sampling in depth increases with velocity.

uCT
wavelength variation with depth
(Left) and vertical time
(Right) of a
Hz sine. The velocity profile is
.
Figure 1. | |
---|---|

Figure 1 shows a comparison of vertical sampling in depth and vertical time . In the Cartesian domain on the left, the sampling of the wavefield is relatively coarse in shallow layers and becomes finer with depth. In the domain on the right, the wavefield is evenly sampled in spite of the velocity variation, for the same number of samples.

Equation 2 maps a depth point to vertical time point . The inverse mapping is also straighforward, from differentiation of inverse functions, it follows

where the integration constant is zero because .

By changing the vertical axis from to , the coordinate system is effectively changed from the Cartesian frame to a new coordinate frame , where the two coordinate systems are related by

As long as the funtions and are one-to-one and free from singularities, we can interpolate any space functions between the Cartesian and domains. Figures 2(a) to 3(b) illustrate examples of such interpolations. The velocity field in the domain is obtained by interpolating the velocity in the Cartesian domain using Equation 2. In Figures 2(a) and 2(b), the mapping velocity is computed by stacking the true velocity horizontally, and thus it is laterally constant. As a result, the coordinate system is orthogonal. In Figures 3(a) and 3(b), is computed from the true velocity, i.e. , the resulting coordinate system is nonorthogonal due to the lateral variation of .

linA,linB
A linear velocity model in (a) Cartesian domain and (b) orthogonal
domain overlaid with
domain mesh. The velocity is
.
Figure 2. |
---|

linC,linD
The same velocity model as in Figure 2(a) in (a) Cartesian domain and (b) nonorthogonal
domain overlaid with
domain mesh.
Figure 3. |
---|

The new coordinate system has mixed units of time in the vertical axis and distance in the horizontal axes. Sometime it is more convenient to have distance units in all three axes of the space domain. To achieve this we can simply scale by some velocity funciton

(6) |

For example, if is set to constant, the pseudodepth is simply a linear scaling of , . The importance of this step will be evident later as we look into anisotropic media.

Wavefield extrapolation in pseudodepth domain |

2013-04-02