Prestack processing of land data with complex topography

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1875-1886 ◽  
Author(s):  
Sara Rajasekaran ◽  
George A. McMechan

A new wave‐equation–based prestack seismic processing system is proposed. This system has only two essential elements; velocity analysis and depth migration. This approach applies truly surface‐consistent statics corrections, regardless of the amount of elevation, change or of near‐surface velocity variation. It uses tomography for estimating the details of shallow velocities and a finite‐difference solution of the two‐way wave‐equation both for computation of image times and for data extrapolation in migration. A field data set that violates most of the assumptions in conventional common midpoint (CMP) processing, because of severe elevation changes and near‐surface velocity variations, is successfully processed. The final depth section reveals a complicated fold‐thrust geometry that was not visible after CMP processing.

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.


1997 ◽  
Vol 34 (8) ◽  
pp. 1167-1184 ◽  
Author(s):  
S. Winardhi ◽  
R. F. Mereu

The 1992 Lithoprobe Abitibi–Grenville Seismic Refraction Experiment was conducted using four profiles across the Grenville and Superior provinces of the southeastern Canadian Shield. Delay-time analysis and tomographic inversion of the data set demonstrate significant lateral and vertical variations in crustal velocities from one terrane to another, with the largest velocity values occurring underneath the Central Gneiss and the Central Metasedimentary belts south of the Grenville Front. The Grenville Front Tectonic Zone is imaged as a southeast-dipping region of anomalous velocity gradients extending to the Moho. The velocity-anomaly maps suggest an Archean crust may extend, horizontally, 140 km beneath the northern Grenville Province. Near-surface velocity anomalies correlate well with the known geology. The most prominent of these is the Sudbury Structure, which is well mapped as a low-velocity basinal structure. The tomography images also suggest underthrusting of the Pontiac and Quetico subprovinces beneath the Abitibi Greenstone Belt. Wide-angle PmP signals, indicate that the Moho varies from a sharp discontinuity south of the Grenville Front to a rather diffuse and flat boundary under the Abitibi Greenstone Belt north of the Grenville Front. A significant crustal thinning near the Grenville Front may indicate post-Grenvillian rebound and (or) the extensional structure of the Ottawa–Bonnechere graben. Crustal thickening resulting from continental collision may explain the tomographic images showing the Moho is 4–5 km deeper south of the Grenville Front.


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 ◽  
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 39 (5) ◽  
pp. 324-331
Author(s):  
Gary Murphy ◽  
Vanessa Brown ◽  
Denes Vigh

As part of a wide-reaching full-waveform inversion (FWI) research program, FWI is applied to an onshore seismic data set collected in the Delaware Basin, west Texas. FWI is routinely applied on typical marine data sets with high signal-to-noise ratio (S/N), relatively good low-frequency content, and reasonably long offsets. Land seismic data sets, in comparison, present significant challenges for FWI due to low S/N, a dearth of low frequencies, and limited offsets. Recent advancements in FWI overcome limitations due to poor S/N and low frequencies making land FWI feasible to use to update the shallow velocities. The chosen area has contrasting and variable near-surface conditions providing an excellent test data set on which to demonstrate the workflow and its challenges. An acoustic FWI workflow is used to update the near-surface velocity model in order to improve the deeper image and simultaneously help highlight potential shallow drilling hazards.


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.


2018 ◽  
Vol 64 (245) ◽  
pp. 450-464 ◽  
Author(s):  
COLLEEN A. MORTIMER ◽  
MARTIN SHARP ◽  
WESLEY VAN WYCHEN

ABSTRACTRepeat airborne laser altimetry measurements show widespread thinning (surface lowering) of glaciers in Canada's Queen Elizabeth Islands since 1995. Thinning rates averaged for 50 m elevation bins, were more than three times higher during the period 2005/06 to 2012/14 pentad than during the previous two pentads. Strongly negative thickness change (dh/dt) anomalies from 2005/06 to 2012/14, relative to the 1995–2012/14 mean, suggest that most of the measured thinning occurred during the most recent 5–6 year period when mean summer land surface temperatures (LSTs) were anomalously high and the mean summer black-sky shortwave broadband albedos (BSA) were anomalously low, relative to the 2000/01–15/16 period, and upper-air (700 hPa) and near surface (2 m) air temperatures were between 0.8°C and 1.5°C higher than 1995–2012 mean. Comparisons of dh/dt with mean summer LST and BSA measurements from the Moderate Resolution Imaging Spectroradiometer and with surface longitudinal strain rates computed from surface velocity fields derived from RADARSAT 1/2 and Landat-7 ETM + data suggest that surface elevation changes were driven mainly by changes in climate. An exception to this occurs along many fast-flowing outlet glaciers where ice dynamics appear also to have played an important role in surface elevation changes.


Geophysics ◽  
2021 ◽  
pp. 1-45
Author(s):  
Guofeng Liu ◽  
Xiaohong Meng ◽  
Johanes Gedo Sea

Seismic reflection is a proven and effective method commonly used during the exploration of deep mineral deposits in Fujian, China. In seismic data processing, rugged depth migration based on wave-equation migration can play a key role in handling surface fluctuations and complex underground structures. Because wave-equation migration in the shot domain cannot output offset-domain common-image gathers in a straightforward way, the use of traditional tools for updating the velocity model and improving image quality can be quite challenging. To overcome this problem, we employed the attribute migration method. This worked by sorting the migrated stack results for every single-shot gather into the offset gathers. The value of the offset that corresponded to each image point was obtained from the ratio of the original migration results to the offset-modulated shot-data migration results. A Gaussian function was proposed to map every image point to a certain range of offsets. This helped improve the signal-to-noise ratio, which was especially important in handing low quality seismic data obtained during mineral exploration. Residual velocity analysis was applied to these gathers to update the velocity model and improve image quality. The offset-domain common-image gathers were also used directly for real mineral exploration seismic data with rugged depth migration. After several iterations of migration and updating the velocity, the proposed procedure achieved an image quality better than the one obtained with the initial velocity model. The results can help with the interpretation of thrust faults and deep deposit exploration.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. B153-B165 ◽  
Author(s):  
Kai Yang ◽  
Hong-Ming Zheng ◽  
Li Wang ◽  
Yu-Zhu Liu ◽  
Fan Jiang ◽  
...  

An integrated wave-equation datuming scheme improves the imaging quality of seismic data from overthrust areas. It can be regarded as integrated because upward-layer replacement is included. In this scheme, data are downward continued to a nonplanar datum (such as the base of the weathering layer), followed by upward continuation from the nonplanar datum to a final planar datum using a one-way extrapolator. When compared with a Kirchhoff integral, this method can deal better with the strong lateral velocity variation within the near surface. After a test on synthetic data, the scheme is applied successfully to real 2D overthrust data acquired in the Qi-Lian foothills, western China. Compared with results using static corrections, integrated wave-equation datuming results lead to better reconstruction of the diffractions and reflections, more reliable migration-velocity analyses, and stronger stack and final depth images.


Sign in / Sign up

Export Citation Format

Share Document