Where is the zero‐velocity layer?

Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 266-269 ◽  
Author(s):  
Samuel H. Gray

The zero‐velocity layer was introduced in Higginbotham et al. (1985) to increase the maximum dip imaging capability of finite‐difference depth migration. Beasley and Lynn (1992) adapted the idea to improve the imaging, again using finite‐difference depth migration, of seismic data acquired in areas of irregular topography. Beasley and Lynn's application improves upon the conventional method of processing, which is to time shift the data from the acquisition surface to a horizontal datum, and then migrate using the near‐surface velocity above the surface and the best estimate of seismic velocity below the surface. The conventional procedure typically produces artifacts in the shallow part of the section that are characteristic of overmigration. To reduce these artifacts, velocities are often reduced for the migration step. The use of the zero‐velocity layer overcomes the need to adjust the migration velocities. Here, a component of the migration velocity is set to zero in the layer between the datum and the surface. The function of the zero‐velocity layer in migration is to remove the elevation‐static correction applied in shifting the data to the flat datum. Only after the data have migrated through the zero‐velocity layer to the irregular recording surface does the migration begin to act in its customary sense, moving energy from trace to trace.

Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 25-38 ◽  
Author(s):  
Xianhuai Zhu ◽  
Burke G. Angstman ◽  
David P. Sixta

Through the use of iterative turning‐ray tomography followed by wave‐equation datuming (or tomo‐datuming) and prestack depth migration, we generate accurate prestack images of seismic data in overthrust areas containing both highly variable near‐surface velocities and rough topography. In tomo‐datuming, we downward continue shot records from the topography to a horizontal datum using velocities estimated from tomography. Turning‐ray tomography often provides a more accurate near‐surface velocity model than that from refraction statics. The main advantage of tomo‐datuming over tomo‐statics (tomography plus static corrections) or refraction statics is that instead of applying a vertical time‐shift to the data, tomo‐datuming propagates the recorded wavefield to the new datum. We find that tomo‐datuming better reconstructs diffractions and reflections, subsequently providing better images after migration. In the datuming process, we use a recursive finite‐difference (FD) scheme to extrapolate wavefield without applying the imaging condition, such that lateral velocity variations can be handled properly and approximations in traveltime calculations associated with the raypath distortions near the surface for migration are avoided. We follow the downward continuation step with a conventional Kirchhoff prestack depth migration. This results in better images than those migrated from the topography using the conventional Kirchhoff method with traveltime calculation in the complicated near surface. Since FD datuming is only applied to the shallow part of the section, its cost is much less than the whole volume FD migration. This is attractive because (1) prestack depth migration usually is used iteratively to build a velocity model, so both efficiency and accuracy are important factors to be considered; and (2) tomo‐datuming can improve the signal‐to‐noise (S/N) ratio of prestack gathers, leading to more accurate migration velocity analysis and better images after depth migration. Case studies with synthetic and field data examples show that tomo‐datuming is especially helpful when strong lateral velocity variations are present below the topography.


Geophysics ◽  
1992 ◽  
Vol 57 (11) ◽  
pp. 1435-1443 ◽  
Author(s):  
Craig Beasley ◽  
Walt Lynn

Seismic data acquired in areas with irregular topography are usually corrected to a flat datum before migration. A time‐honored technique for handling elevation changes is to time shift the data before application of migration. This simple time shift, or elevation‐static correction, cannot properly represent wide‐angle or dipping reflections as they would have been recorded at the datum. As a result, when elevation varies significantly, accuracy in event positioning may be compromised for migration and other wave‐equation processes, such as dip moveout processing (DMO). Traditionally, such over‐ and under‐migration artifacts have been dealt with by increasing or decreasing the migration velocity. However, simple adjustment of the migration velocity cannot undo the wave‐field distortions induced in seismic data acquired over varying elevations. More sophisticated and accurate solutions such as wave‐equation datuming are too computationally demanding for routine use. Here, we propose an efficient and accurate technique for doing migration from irregular surfaces using conventional migration algorithms. As in elevation‐static corrections, surface‐recorded data are time‐shifted to a horizontal datum; for our process, we choose to have that datum elevation lie at or above the highest elevation in the survey. After migration, the datum elevation can always be adjusted to any other level by means of a bulk time shift. In the migration step, the velocity is set to zero (or some very small value) in the layer between the surface and the datum; below the original surface, the interval velocity represents the best estimate of the subsurface geology. By adding a zero‐velocity layer, the migration algorithm is applied to the data from the flat datum and no lateral propagation is allowed until a nonzero velocity is encountered at the recording surface. Synthetic and field data examples demonstrate that use of the “zero‐velocity layer” significantly improves imaging accuracy relative to conventional migration from a flat datum. Moreover, the geologically derived migration‐velocity field need not be adjusted to compensate for shortcomings in the datum‐static procedure. The technique can be extended to prestack processes such as DMO, shot‐ and receiver‐gather downward extrapolation, and migration and thus suggests a unified approach to processing data from irregular surfaces.


2021 ◽  
Vol 13 (14) ◽  
pp. 2684
Author(s):  
Eldert Fokker ◽  
Elmer Ruigrok ◽  
Rhys Hawkins ◽  
Jeannot Trampert

Previous studies examining the relationship between the groundwater table and seismic velocities have been guided by empirical relationships only. Here, we develop a physics-based model relating fluctuations in groundwater table and pore pressure with seismic velocity variations through changes in effective stress. This model justifies the use of seismic velocity variations for monitoring of the pore pressure. Using a subset of the Groningen seismic network, near-surface velocity changes are estimated over a four-year period, using passive image interferometry. The same velocity changes are predicted by applying the newly derived theory to pressure-head recordings. It is demonstrated that the theory provides a close match of the observed seismic velocity changes.


Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. Q27-Q37
Author(s):  
Yang Shen ◽  
Jie Zhang

Refraction methods are often applied to model and image near-surface velocity structures. However, near-surface imaging is very challenging, and no single method can resolve all of the land seismic problems across the world. In addition, deep interfaces are difficult to image from land reflection data due to the associated low signal-to-noise ratio. Following previous research, we have developed a refraction wavefield migration method for imaging shallow and deep interfaces via interferometry. Our method includes two steps: converting refractions into virtual reflection gathers and then applying a prestack depth migration method to produce interface images from the virtual reflection gathers. With a regular recording offset of approximately 3 km, this approach produces an image of a shallow interface within the top 1 km. If the recording offset is very long, the refractions may follow a deep path, and the result may reveal a deep interface. We determine several factors that affect the imaging results using synthetics. We also apply the novel method to one data set with regular recording offsets and another with far offsets; both cases produce sharp images, which are further verified by conventional reflection imaging. This method can be applied as a promising imaging tool when handling practical cases involving data with excessively weak or missing reflections but available refractions.


Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 721-732 ◽  
Author(s):  
Lanlan Yan ◽  
Larry R. Lines

Seismic imaging of complex structures from the western Canadian Foothills can be achieved by applying the closely coupled processes of velocity analysis and depth migration. For the purposes of defining these structures in the Shaw Basing area of western Alberta, we performed a series of tests on both synthetic and real data to find optimum imaging procedures for handling large topographic relief, near‐surface velocity variations, and the complex structural geology of steeply dipping formations. To better understand the seismic processing problems, we constructed a typical foothills geological model that included thrust faults and duplex structures, computed the model responses, and then compared the performance of different migration algorithms, including the explicit finite difference (f-x) and Kirchhoff integral methods. When the correct velocity was used in the migration tests, the f-x method was the most effective in migration from topography. In cases where the velocity model was not assumed known, we determined a macrovelocity model by performing migration/velocity analysis by using smiles and frowns in common image gathers and by using depth‐focusing analysis. In applying depth imaging to the seismic survey from the Shaw Basing area, we found that imaging problems were caused partly by near‐surface velocity problems, which were not anticipated in the modeling study. Several comparisons of different migration approaches for these data indicated that prestack depth migration from topography provided the best imaging results when near‐surface velocity information was incorporated. Through iterative and interpretive migration/velocity analysis, we built a macrovelocity model for the final prestack depth migration.


2020 ◽  
Vol 8 (3) ◽  
pp. T651-T665
Author(s):  
Yalin Li ◽  
Xianhuai Zhu ◽  
Gengxin Peng ◽  
Liansheng Liu ◽  
Wensheng Duan

Seismic imaging in foothills areas is challenging because of the complexity of the near-surface and subsurface structures. Single seismic surveys often are not adequate in a foothill-exploration area, and multiple phases with different acquisition designs within the same block are required over time to get desired sampling in space and azimuths for optimizing noise attenuation, velocity estimation, and migration. This is partly because of economic concerns, and it is partly because technology is progressing over time, creating the need for unified criteria in processing workflows and parameters at different blocks in a study area. Each block is defined as a function of not only location but also the acquisition and processing phase. An innovative idea for complex foothills seismic imaging is presented to solve a matrix of blocks and tasks. For each task, such as near-surface velocity estimation and static corrections, signal processing, prestack time migration, velocity-model building, and prestack depth migration, one or two best service companies are selected to work on all blocks. We have implemented streamlined processing efficiently so that Task-1 to Task-n progressed with good coordination. Application of this innovative approach to a mega-project containing 16 3D surveys covering more than [Formula: see text] in the Kelasu foothills, northwestern China, has demonstrated that this innovative approach is a current best practice in complex foothills imaging. To date, this is the largest foothills imaging project in the world. The case study in Kelasu successfully has delivered near-surface velocity models using first arrivals picked up to 3500 m offset for static corrections and 9000 m offset for prestack depth migration from topography. Most importantly, the present megaproject is a merge of several 3D surveys, with the merge performed in a coordinated, systematic fashion in contrast to most land megaprojects. The benefits of this approach and the strategies used in processing data from the various subsurveys are significant. The main achievement from the case study is that the depth images, after the application of the near-surface velocity model estimated from the megasurveys, are more continuous and geologically plausible, leading to more accurate seismic interpretation.


Geophysics ◽  
1980 ◽  
Vol 45 (3) ◽  
pp. 361-375 ◽  
Author(s):  
D. R. Judson ◽  
J. Lin ◽  
P. S. Schultz ◽  
J. W. C. Sherwood

The conventional methods for migrating a seismic section, e.g., the finite‐difference method and the Kirchhoff summation method, are inadequate in the presence of significant lateral variations in velocity. For this type of velocity distribution, the basic migration output should be in true depth, although for practical purposes it may be preferable to display it with a nonlinear depth scale. A finite‐difference method has been implemented for obtaining migrated depth sections. The concept underlying this involves all the usual assumptions of a dip line and primary reflections only, with the seismic section considered as the surface measurement of an upcoming wave field which we process with downward continuation in small increments of depth, rather than the customary increments of traveltime. The specified velocity variation laterally along a thin layer results in transmission time changes which must be corrected by a small static time shift applied to each seismic trace. This additional operation within the migration algorithm can be difficult and expensive to implement and is the main reason for its prior omission. Results are given of depth migration applications to both synthetic and real seismic data.


Geophysics ◽  
2008 ◽  
Vol 73 (2) ◽  
pp. S27-S34 ◽  
Author(s):  
Tong W. Fei ◽  
Christopher L. Liner

When a subsurface is anisotropic, migration based on the assumption of isotropy will not produce accurate migration images. We develop a hybrid wave-equation migration algorithm for vertical transversely isotropic (VTI) media based on a one-way acoustic wave equation, using a combination of Fourier finite-difference (FFD) and finite-difference (FD) approaches. The hybrid method can suppress an additional solution that exists in the VTI acoustic wave equation, and it offers speed and other advantages over conventional FFD or FD methods alone. The algorithm is tested on a synthetic model involving log data from onshore eastern Saudi Arabia, including estimates of both intrinsic and layer-induced VTI parameters. Results indicate that VTI imaging in this region offers some improvement over isotropic imaging, primarily with respect to subtle structure and stratigraphy and to image continuity. These benefits probably will be overshadowed by perennial land seismic data issues such as near-surface distortions and multiples.


Sign in / Sign up

Export Citation Format

Share Document