Estimation of Near-Surface Velocity for 3-D Complicated Topography and Seismic Tomographic Static Corrections

2001 ◽  
Vol 44 (2) ◽  
pp. 268-274 ◽  
Author(s):  
Yi-Ke LIU ◽  
Xu CHANG ◽  
Hui WANG ◽  
Fu-Zhong LI
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.


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 ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. S117-S129 ◽  
Author(s):  
Jianfeng Zhang ◽  
Jincheng Xu ◽  
Hao Zhang

We have developed a modified 3D prestack time migration (PSTM) scheme that can handle rugged topography as well as high near-surface velocities in land seismic imaging. The proposed topography PSTM can be applied to seismic data recorded on a 3D irregular surface without static corrections. Two effective velocity parameters were found to describe wave propagation through inhomogeneous media above and below a chosen datum. As a result, wave propagation phenomena in the complex near surface, such as near-vertical incidences through a weathering layer and raypaths bending away from vertical in the presence of high near-surface velocities, are correctly considered. The two effective velocity parameters can be estimated by flattening events in imaging gathers. Hence, it is not necessary to have detailed knowledge of the near-surface velocity model and velocity field below the datum when applying topography PSTM. We integrated residual static corrections into topography PSTM. This eliminated the distortions along the events better than conventional residual static corrections, which are usually applied before migration. The computational cost of the topography PSTM was only slightly higher than that of conventional PSTM due to the use of a table-driven algorithm. Three-dimensional synthetic and field data sets were used to test the proposed topography PSTM. High-quality imaging results were obtained.


Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 1082-1090 ◽  
Author(s):  
Björn Bergman ◽  
Ari Tryggvason ◽  
Christopher Juhlin

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.


Geophysics ◽  
1990 ◽  
Vol 55 (6) ◽  
pp. 757-760 ◽  
Author(s):  
G. A. McMechan ◽  
H. W. Chen

Static effects due to surface topography and near‐surface velocity variations may be accurately compensated for, in an implicit way, during prestack reverse‐time migration of common‐source gathers, obviating the need for explicit static corrections. Receiver statics are incorporated by extrapolating the observed data from the actual recorder positions; source statics are incorporated by computing the excitation‐time imaging conditions from the actual source positions.


Geophysics ◽  
2007 ◽  
Vol 72 (4) ◽  
pp. U55-U66 ◽  
Author(s):  
Robbert van Vossen ◽  
Jeannot Trampert

Near-surface wavefield perturbations can be very complex and completely mask the target reflections. Despite this complexity, conventional methods rely on parameterizations characterized by simple time and amplitude anomalies to compensate for these perturbations. Determining and compensating for time shifts is generally referred to as (residual) static corrections, whereas surface-consistent deconvolution techniques deal with amplitude anomalies. We present an approach that uses the full waveform to parameterize near-surface perturbations. Therefore, we refer to this method as waveform statics. Important differences from conventional static corrections are that this approach allows time shifts to vary with frequency and takes amplitude variations directly into account. Furthermore, the procedure is fully automated and does not rely on near-surface velocity information. The waveform static corrections are obtained usingblind channel identification and applied to the recordings using multichannel deconvolution. As a result, the method implicitly incorporates array forming. The developed method is validated on synthetic data and applied to part of a field data set acquired in an area with significant near-surface heterogeneity. The source and receiver responses obtained are strongly correlated to the near-surface conditions and show changes, both in phase and frequency content, along the spread. The application of the waveform statics demonstrates that they not only correct for near-surface wavefield perturbations, but also strongly reduce coherent noise. This results in substantial improvements, both in trace-to-trace coherency and in depth resolution. In addition, the procedure delineates reflection events that are difficult to detect prior to our proposed correction. Based on these results, we conclude that complex near-surface perturbations can be successfully dealt with using the multichannel, full-waveform, static-correction procedure.


Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. U23-U34
Author(s):  
Raul Cova ◽  
David Henley ◽  
Kristopher A. Innanen

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.


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCB35-WCB46 ◽  
Author(s):  
Hua-wei Zhou ◽  
Peiming Li ◽  
Zhihui Yan ◽  
Hui Liu

Although first-arrival tomography provides an effective way to estimate near-surface velocities and static corrections, the undulation of velocity interfaces such as the base of the weathered zone may not be easily determined by this method. The main reason is that first arrivals are insensitive to small geometric changes in velocity interfaces because their raypaths tend to traverse along those interfaces. To improve the solution of interface geometry, we developed a deformable layer tomostatics method that approximates the near-surface velocity field as several layers of constant velocity and variable thickness that can be inverted for the geometry of the velocity interfaces. We use a multiscale model parameterization in the inversion for interface geometry. Synthetic and field data tests showed that the method can determine the interface geometry. Constraining the depth range of the basal boundary of the weathered zone increases the convergence rate of the iterative inversion process. Tests on field data showed greater reflection coherency in a stacked section based on constrained static corrections than in one from unconstrained static corrections. The method yielded a better match with statics computed from sand-dune curves than does a match obtained by using two commercial grid tomography packages.


Sign in / Sign up

Export Citation Format

Share Document