Overthrust imaging with tomo‐datuming: A case study

Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 25-38 ◽  
Author(s):  
Xianhuai Zhu ◽  
Burke G. Angstman ◽  
David P. Sixta
Keyword(s):  
Velocity Model ◽  
Surface Velocity ◽  
Time Shift ◽  
Lateral Velocity ◽  
Prestack Depth Migration ◽  
Near Surface ◽  
Depth Migration ◽  
Static Corrections ◽  
Velocity Variations ◽  
Kirchhoff Method

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.

Novel strategies for complex foothills seismic imaging — Part 1: Mega-near-surface velocity estimation

2020 ◽  
Vol 8 (3) ◽  
pp. T651-T665
Author(s):  
Yalin Li ◽  
Xianhuai Zhu ◽  
Gengxin Peng ◽  
Liansheng Liu ◽  
Wensheng Duan
Keyword(s):  
Seismic Imaging ◽  
Velocity Model ◽  
Velocity Estimation ◽  
Innovative Approach ◽  
Surface Velocity ◽  
Near Surface ◽  
Depth Migration ◽  
Static Corrections ◽  
Over Time

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.

High‐resolution seismic traveltime tomography incorporating static corrections applied to a till‐covered bedrock environment

Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 1082-1090 ◽  
Author(s):  
Björn Bergman ◽  
Ari Tryggvason ◽  
Christopher Juhlin
Keyword(s):  
High Resolution ◽  
Velocity Model ◽  
Surface Velocity ◽  
Data Set ◽  
Near Surface ◽  
Seismic Processing ◽  
Simultaneous Inversion ◽  
Reflection Seismic ◽  
Static Corrections ◽  
Velocity Variations

A major obstacle in tomographic inversion is near‐surface velocity variations. Such shallow velocity variations need to be known and correctly accounted for to obtain images of deeper structures with high resolution and quality. Bedrock cover in many areas consists of unconsolidated sediments and glacial till. To handle the problems associated with this cover, we present a tomographic method that solves for the 3D velocity structure and receiver static corrections simultaneously. We test the method on first‐arrival picks from deep seismic reflection data acquired in the mid‐ late to 1980s in the Siljan Ring area, central Sweden. To use this data set successfully, one needs to handle a number of problems, including time‐varying, near‐surface velocities from data recorded in winter and summer, several sources and receivers within each inversion cell, varying thickness of the cover layer in each inversion cell, and complex 3D geology. Simultaneous inversion for static corrections and velocity produces a much better image than standard tomography without statics. The velocity model from the simultaneous inversion is superior to the velocity model produced using refraction statics obtained from standard reflection seismic processing prior to inversion. Best results using the simultaneous inversion are obtained when the initial top velocity layer is set to the near‐surface bedrock velocity rather than the velocity of the cover. The resulting static calculations may, in the future, be compared to refraction static corrections in standard reflection seismic processing. The preferred final model shows a good correlation with the mapped geology and the airborne magneticmap.

Recent applications of turning-ray tomography

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE243-VE254 ◽  
Author(s):  
Xianhuai Zhu ◽  
Paul Valasek ◽  
Baishali Roy ◽  
Simon Shaw ◽  
Jack Howell ◽  
...  
Keyword(s):  
Reservoir Characterization ◽  
Velocity Model ◽  
Surface Velocity ◽  
Seismic Survey ◽  
Structural Imaging ◽  
Prestack Depth Migration ◽  
Near Surface ◽  
Depth Migration ◽  
Time Migration ◽  
Prestack Time Migration

Recent applications of 2D and 3D turning-ray tomography show that near-surface velocities are important for structural imaging and reservoir characterization. For structural imaging, we used turning-ray tomography to estimate the near-surface velocities for static corrections followed by prestack time migration and the near-surface velocities for prestack depth migration. Two-dimensional acoustic finite-difference modeling illustrates that wave-equation prestack depth migration is very sensitive to the near-surface velocities. Field data demonstrate that turning-ray tomography followed by prestack time migration helps to produce superior images in complex geologic settings. When the near-surface velocity model is integrated into a background velocity model for prestack depth migration, we find that wave propagation is very sensitive to the velocities immediately below the topography. For shallow-reservoir characterization, we developed and applied azimuthal turning-ray tomography to investigate observed apparent azimuthal-traveltime variations, using a wide-azimuth land seismic survey from a heavy-oil field at Surmont, Canada. We found that the apparent azimuthal velocity variations are not necessarily related to azimuthal anisotropy, or horizontal transverse isotropy (HTI), induced by the stress field or fractures. Near-surface heterogeneity and the acquisition footprint also could result in apparent azimuthal variations.

Tomographic velocity analysis and wave-equation depth migration in an overthrust terrain: A case study from the Tuha Basin, China

2018 ◽  
Vol 6 (1) ◽  
pp. T1-T13
Author(s):  
Bin Lyu ◽  
Qin Su ◽  
Kurt J. Marfurt
Keyword(s):  
Wave Equation ◽  
Data Quality ◽  
Model Building ◽  
Velocity Model ◽  
Lateral Velocity ◽  
Near Surface ◽  
Depth Migration ◽  
Time Migration ◽  
Head Waves ◽  
Velocity Variations

Although the structures associated with overthrust terrains form important targets in many basins, accurately imaging remains challenging. Steep dips and strong lateral velocity variations associated with these complex structures require prestack depth migration instead of simpler time migration. The associated rough topography, coupled with older, more indurated, and thus high-velocity rocks near or outcropping at the surface often lead to seismic data that suffer from severe statics problems, strong head waves, and backscattered energy from the shallow section, giving rise to a low signal-to-noise ratio that increases the difficulties in building an accurate velocity model for subsequent depth migration. We applied a multidomain cascaded noise attenuation workflow to suppress much of the linear noise. Strong lateral velocity variations occur not only at depth but near the surface as well, distorting the reflections and degrading all deeper images. Conventional elevation corrections followed by refraction statics methods fail in these areas due to poor data quality and the absence of a continuous refracting surface. Although a seismically derived tomographic solution provides an improved image, constraining the solution to the near-surface depth-domain interval velocities measured along the surface outcrop data provides further improvement. Although a one-way wave-equation migration algorithm accounts for the strong lateral velocity variations and complicated structures at depth, modifying the algorithm to account for lateral variation in illumination caused by the irregular topography significantly improves the image, preserving the subsurface amplitude variations. We believe that our step-by-step workflow of addressing the data quality, velocity model building, and seismic imaging developed for the Tuha Basin of China can be applied to other overthrust plays in other parts of the world.

A strategy for computing seismic attributes on depth-migrated data using time-shift gathers

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. O37-O46
Author(s):  
Mark Roberts
Keyword(s):  
Interpolation Method ◽  
Seismic Attributes ◽  
Depth Image ◽  
Time Shift ◽  
Lateral Velocity ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Velocity Variations ◽  
Made In

Assumptions made during postprocessing can be as important as those made during migration. Prestack depth migration (PSDM) is often used due to its ability to handle lateral velocity variations and dipping events. However, most postprocessing flows still use a simple 1D “depth-to-time” vertical stretch, violating the very assumptions that led us to use PSDM. Postprocessing workflows based on time-shift depth image gathers allow for postprocessing flows that make no further approximations other than those made in migration. For cases in which time-shift gathers are not computed during migration, they can be approximated by use of the “exploding-reflector” model or through an orthogonal shift and interpolation method.

Seismic imaging and velocity analysis for an Alberta Foothills seismic survey

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 721-732 ◽  
Author(s):  
Lanlan Yan ◽  
Larry R. Lines
Keyword(s):  
Seismic Imaging ◽  
Migration Velocity ◽  
Surface Velocity ◽  
Seismic Survey ◽  
Velocity Analysis ◽  
Prestack Depth Migration ◽  
Near Surface ◽  
Depth Migration ◽  
Migration Velocity Analysis ◽  
Imaging Results

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.

Response by Arthur E. Barnes to the reply by the authors

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1947-1947 ◽  
Author(s):  
Arthur E. Barnes
Keyword(s):  
Lateral Velocity ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Time Migration ◽  
Pulse Distortion ◽  
Velocity Variations ◽  
Zero Offset ◽  
Complex Geology

I appreciate the thoughtful and thorough response given by Tygel et al. They point out that even for a single dipping reflector imaged by a single non‐zero offset raypath, pulse distortion caused by “standard processing” (NM0 correction‐CMP sort‐stack‐time migration) and pulse distortion caused by prestack depth migration are not really the same, because the reflecting point is mispositioned in standard processing. Within a CMP gather, this mispositioning increases with offset, giving rise to “CMP smear.” CMP smear degrades the stack, introducing additional pulse distortion. Where i‐t is significant, and where lateral velocity variations or reflection curvature are large, such as for complex geology, the pulse distortion of standard processing can differ greatly from that of prestack depth migration.

Where is the zero‐velocity layer?

Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 266-269 ◽  
Author(s):  
Samuel H. Gray
Keyword(s):  
Finite Difference ◽  
Seismic Velocity ◽  
Surface Velocity ◽  
Time Shift ◽  
Near Surface ◽  
Depth Migration ◽  
Static Correction ◽  
Zero Velocity ◽  
Conventional Procedure ◽  
Velocity Layer

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.

Computing near-surface velocity models for S-wave static corrections using raypath interferometry

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. U23-U34
Author(s):  
Raul Cova ◽  
David Henley ◽  
Kristopher A. Innanen
Keyword(s):  
Wave Velocity ◽  
Velocity Model ◽  
P Wave ◽  
Surface Velocity ◽  
S Wave ◽  
Converted Wave ◽  
Near Surface ◽  
Velocity Models ◽  
Consistent Solution ◽  
Static Corrections

A near-surface velocity model is one of the typical products generated when computing static corrections, particularly in the processing of PP data. Critically refracted waves are the input usually needed for this process. In addition, for the converted PS mode, S-wave near-surface corrections must be applied at the receiver locations. In this case, however, critically refracted S-waves are difficult to identify when using P-wave energy sources. We use the [Formula: see text]-[Formula: see text] representation of the converted-wave data to capture the intercept-time differences between receiver locations. These [Formula: see text]-differences are then used in the inversion of a near-surface S-wave velocity model. Our processing workflow provides not only a set of raypath-dependent S-wave static corrections, but also a velocity model that is based on those corrections. Our computed near-surface S-wave velocity model can be used for building migration velocity models or to initialize elastic full-waveform inversions. Our tests on synthetic and field data provided superior results to those obtained by using a surface-consistent solution.

Scattering effect on shallow gas-obscured zone imaging in Bohai PL19-3 area

Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. B43-B53 ◽  
Author(s):  
Xianhuai Zhu ◽  
Kirk Wallace ◽  
Qingrong Zhu ◽  
Robert Hofer
Keyword(s):  
Velocity Model ◽  
Data Depth ◽  
Surface Velocity ◽  
Bohai Bay ◽  
Shallow Gas ◽  
Near Surface ◽  
Depth Migration ◽  
Precise Knowledge ◽  
Streamer Data ◽  
Viscoelastic Modeling

Seismic imaging is challenging in the Bohai Bay PL19-3 area, offshore China. Bohai Bay Field is seismically obscured, making well penetrations the only reliable source of data for subsurface interpretation. With the help of turning-ray tomography, we are able to obtain a reliable near-surface velocity model approximately down to 700 m below sea level, using the first arrivals picked from streamer data. Depth migration using velocities estimated from turning-ray tomography has improved shallow structures and fault definition. However, reservoir level structures from 800 to 1500 m are still poorly imaged. A viscoelastic modeling study with assigned variable Q and shallow velocity profiles, with and without shallow gas-induced scatterers, demonstrates that scattering is the primary controlling phenomenon causing imaging difficulty within the obscured zone. Due to scattering, imaging tests at the target level were unsuccessful even with precise knowledge of velocity.

Sign in / Sign up

Export Citation Format

Share Document