IMAGING BENEATH COMPLEX STRUCTURE: A CASE HISTORY

1981 ◽  
Vol 21 (1) ◽  
pp. 112
Author(s):  
K. Lamer ◽  
B. Gibson ◽  
R. Chambers
Keyword(s):  
Velocity Structure ◽  
Deep Structure ◽  
Complex Structure ◽  
Velocity Model ◽  
Important Conclusion ◽  
Depth Migration ◽  
Interval Velocity ◽  
Velocity Models ◽  
Time Migration ◽  
Zero Offset

Migration is recognised as the essential step in converting seismic, data into a representation of the earth's subsurface structure. Ironically, conventional migration often fails where migration is needed most—when the data are recorded over complex structures. Processing field data shot in Central America, and synthetic data derived for that section, demonstrates that time migration actually degrades the image of the deep structure that lies below a complicated overburden.In the Central American example, velocities increase nearly two-fold across an arched and thrust-faulted interface. Wavefront distortion introduced by this feature gives rise to distorted reflections from depth. Even with interval velocity known perfectly, no velocity is proper for time migrating the data here; time migration is the wrong process because it does not honour Snell's Law. Depth migration of the stacked data, on the other hand, produces a reasonable image of the deeper section. The depth migration, however, leaves artifacts that could be attributed to problems that are common in structurally complicated areas: (1) departures of the stacked section from the ideal, a zero-offset section; (2) incorrect specification of velocities; and (3) loss of energy transmitted through the complex zoneFor such an inhomogeneous velocity structure, shortcomings in CDP stacking are directly related to highly non- hyperbolic moveout. As with migration velocity, no proper stacking velocity can be developed for these data, even from the known interval-velocity model. Proper treatment of nonzero-offset reflection data could be accomplished by depth migration before stacking. Simple ray-theoretical correction of the complex moveouts, however, can produce a stack that is similar to the desired zero-offset section.Overall, the choice of velocity model most strongly influences the results of depth migration. Processing the data with a range of plausible velocity models, however, leads to an important conclusion: although the velocities can never be known exactly, depth migration is essential for clarifying structure beneath complex overburden.

Incorporating geologic information into reflection tomography

Geophysics ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 533-546 ◽  
Author(s):  
Robert G. Clapp ◽  
Biondo L. Biondi ◽  
Jon F. Claerbout
Keyword(s):  
Velocity Structure ◽  
Null Space ◽  
Velocity Model ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Interval Velocity ◽  
Velocity Models ◽  
Regularized Inversion ◽  
Symmetric Operators ◽  
Reflection Tomography

In areas of complex geology, prestack depth migration is often necessary if we are to produce an accurate image of the subsurface. Prestack depth migration requires an accurate interval velocity model. With few exceptions, the subsurface velocities are not known beforehand and should be estimated. When the velocity structure is complex, with significant lateral variations, reflection‐tomography methods are often an effective tool for improving the velocity estimate. Unfortunately, reflection tomography often converges slowly, to a model that is geologically unreasonable, or it does not converge at all. The large null space of reflection‐tomography problems often forces us to add a sparse parameterization of the model and/or regularization criteria to the estimation. Standard tomography schemes tend to create isotropic features in velocity models that are inconsistent with geology. These isotropic features result, in large part, from using symmetric regularization operators or from choosing a poor model parameterization. If we replace the symmetric operators with nonstationary operators that tend to spread information along structural dips, the tomography will produce velocity models that are geologically more reasonable. In addition, by forming the operators in helical 1D space and performing polynomial division, we apply the inverse of these space‐varying anisotropic operators. The inverse operators can be used as a preconditioner to a standard tomography problem, thereby significantly improving the speed of convergence compared with the typical regularized inversion problem. Results from 2D synthetic and 2D field data are shown. In each case, the velocity obtained improves the focusing of the migrated image.

Depth migration by the Gaussian beam summation method

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S81-S93 ◽  
Author(s):  
Mikhail M. Popov ◽  
Nikolay M. Semtchenok ◽  
Peter M. Popov ◽  
Arie R. Verdel
Keyword(s):  
Wave Equation ◽  
Model Building ◽  
Velocity Structure ◽  
Velocity Model ◽  
Summation Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Depth Migration ◽  
Time Migration ◽  
Ray Methods

Seismic depth migration aims to produce an image of seismic reflection interfaces. Ray methods are suitable for subsurface target-oriented imaging and are less costly compared to two-way wave-equation-based migration, but break down in cases when a complex velocity structure gives rise to the appearance of caustics. Ray methods also have difficulties in correctly handling the different branches of the wavefront that result from wave propagation through a caustic. On the other hand, migration methods based on the two-way wave equation, referred to as reverse-time migration, are known to be capable of dealing with these problems. However, they are very expensive, especially in the 3D case. It can be prohibitive if many iterations are needed, such as for velocity-model building. Our method relies on the calculation of the Green functions for the classical wave equation by per-forming a summation of Gaussian beams for the direct and back-propagated wavefields. The subsurface image is obtained by cal-culating the coherence between the direct and backpropagated wavefields. To a large extent, our method combines the advantages of the high computational speed of ray-based migration with the high accuracy of reverse-time wave-equation migration because it can overcome problems with caustics, handle all arrivals, yield good images of steep flanks, and is readily extendible to target-oriented implementation. We have demonstrated the quality of our method with several state-of-the-art benchmark subsurface models, which have velocity variations up to a high degree of complexity. Our algorithm is especially suited for efficient imaging of selected subsurface subdomains, which is a large advantage particularly for 3D imaging and velocity-model refinement applications such as subsalt velocity-model improvement. Because our method is also capable of providing highly accurate migration results in structurally complex subsurface settings, we have also included the concept of true-amplitude imaging in our migration technique.

Unified imaging of multichannel seismic data: Combining traveltime inversion and prestack depth migration

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE269-VE280 ◽  
Author(s):  
Priyank Jaiswal ◽  
Colin A. Zelt
Keyword(s):  
Seismic Data ◽  
Joint Inversion ◽  
Velocity Model ◽  
Depth Image ◽  
Prestack Depth Migration ◽  
Traveltime Inversion ◽  
Depth Migration ◽  
Interval Velocity ◽  
Wide Aperture ◽  
Zero Offset

Imaging 2D multichannel land seismic data can be accomplished effectively by a combination of traveltime inversion and prestack depth migration (PSDM), referred to as unified imaging. Unified imaging begins by inverting the direct-arrival times to estimate a velocity model that is used in static corrections and stacking velocity analysis. The interval velocity model (from stacking velocities) is used for PSDM. The stacked data and the PSDM image are interpreted for common horizons, and the corresponding wide-aperture reflections are identified in the shot gathers. Using the interval velocity model, the stack interpretations are inverted as zero-offset reflections to constrain the corresponding interfaces in depth; the interval velocity model remains stationary. We define a coefficient of congruence [Formula: see text] that measures the discrepancy between horizons from the PSDM image andtheir counterparts from the zero-offset inversion. A value of unity for [Formula: see text] implies that the interpreted and inverted horizons are consistent to within the interpretational uncertainties, and the unified imaging is said to have converged. For [Formula: see text] greater than unity, the interval velocity model and the horizon depths are updated by jointly inverting the direct arrivals with the zero-offset and wide-aperture reflections. The updated interval velocity model is used again for both PSDM and a zero-offset inversion. Interpretations of the new PSDM image are the updated horizon depths. The unified imaging is applied to seismic data from the Naga Thrust and Fold Belt in India. Wide-aperture and zero-offset data from three geologically significant horizons are used. Three runs of joint inversion and PSDM are required in a cyclic manner for [Formula: see text] to converge to unity. A joint interpretation of the final velocity model and depth image reveals the presence of a triangle zone that could be promising for exploration.

scholarly journals Method Comparison of Model Based Tomography and Grid Based Tomography to refine interval velocity

2016 ◽  
Vol 4 (01) ◽  
pp. 63
Author(s):  
Yuninggar Dwi Nugroho ◽  
Sudarmaji S
Keyword(s):  
Method Comparison ◽  
Velocity Model ◽  
Seismic Section ◽  
Initial Interval ◽  
Depth Migration ◽  
Interval Velocity ◽  
Model Based ◽  
Time Migration ◽  
Tomography Method ◽  
Grid Based

<span>The input data for pre stack time migration and pre stack depth migration is velocity model. <span>The exact velocity model can provide maximum result in seismic section. The best seismic <span>section can minimize possibility of errors during interpretation. Model based and grid based <span>tomography are used to refine the interval velocity model. The interval velocity will be used as <span>input in the pre stack depth migration. Initial interval velocity is obtained from RMS velocity<br /><span>using Dix formula. This velocity will be refined by global depth tomography method. The <span>global depth tomography method is divided into model based and grid based tomography. <span>Velocity analysis is performed along the horizon (depth model). Residual depth move out is <span>obtained from picking velocity. It is used as input in tomography method. The flat gather is <span>obtained at tenth iteration. The interval velocity that is obtained from tenth iteration has the <span>small errors. Tomography method can provide maximum result on velocity refinement. That is <span>shown by the result that the pre stack depth migration is much better than using initial interval <span>velocity. The pull up effect can be corrected by tomography method.</span></span></span></span></span></span></span></span></span></span></span></span><br /></span>

scholarly journals Zero-offset sections with a deblurring filter in the time domain

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. S239-S249
Author(s):  
Shihang Feng ◽  
Oz Yilmaz ◽  
Yuqing Chen ◽  
Gerard T. Schuster
Keyword(s):  
Time Domain ◽  
Constant Velocity ◽  
Field Data ◽  
Velocity Model ◽  
Velocity Models ◽  
Time Migration ◽  
Prestack Time Migration ◽  
Zero Offset ◽  
Image Volume ◽  
The Time Domain

The conventional common-midpoint stack is not equivalent to the zero-offset section due to the existence of velocity uncertainty. To obtain a zero-offset reflection section that preserves most reflections and diffractions, we have developed a velocity-independent workflow for reconstructing a high-quality zero-offset reflection section from prestack data with a deblurring filter. This workflow constructs a migration image volume by prestack time migration using a series of constant-velocity models. A deblurring filter for each constant-velocity model is applied to each time-migration image to get a deblurred image volume. To preserve all events in the image volume, each deblurred image panel is demigrated and then summed over the velocity axis. Compared with the workflow without a deblurring filter, the composite zero-offset reflection section has higher resolution and fewer migration artifacts. We evaluate applications of our method to synthetic and field data to validate its effectiveness.

scholarly journals Identification of Reservoir Thickness and Estimation of Hydrocarbon Reservation used the Pre-Stack Depth Migration (PSDM) in Cikung Area

2017 ◽  
Vol 1 (01) ◽  
pp. 5-10
Author(s):  
Sudra Irawan ◽  
Hanah Khoirunnisa
Keyword(s):  
Water Saturation ◽  
Complex Structure ◽  
Velocity Model ◽  
Horizontal Velocity ◽  
Geological Data ◽  
Depth Migration ◽  
Interval Velocity ◽  
Hydrocarbon Deposits ◽  
Hydrocarbon Reserves

Information from the geological data shows that the Cikung field has the complex structure which is characterized by carbonate reef. Therefore, it is required in the domain depth migration (PSDM) using interval velocity model is sensitive to variations in vertical and horizontal velocity to identify the thickness of the reservoir and estimation of hydrocarbon reserves. Stages include making the depth structure maps, map making the thickness of the reservoir, and the calculation of hydrocarbon reserves. The results showed that the thickness of the reservoir in Cikung Field, estimated to range from 71 meters to 175 meters with the prospect reservoir is at a depth of 1216 meters to 1247 meters from the surface. Hydrocarbon deposits (oil) is estimated at 1,134 × 106 𝑆𝑇𝐵 or 1,311 × 108 kilo litre by the porosity of 22.6% and 70.7% water saturation.

Kinematic interpretation in the prestack depth‐migrated domain

Geophysics ◽  
1997 ◽  
Vol 62 (4) ◽  
pp. 1226-1237 ◽  
Author(s):  
Irina Apostoiu‐Marin ◽  
Andreas Ehinger
Keyword(s):  
Signal To Noise Ratio ◽  
Velocity Model ◽  
Point Of View ◽  
Prestack Depth Migration ◽  
Depth Migration ◽  
Velocity Models ◽  
Time Domain Data ◽  
And Migration ◽  
Point To Point ◽  
Homogeneous Media

Prestack depth migration can be used in the velocity model estimation process if one succeeds in interpreting depth events obtained with erroneous velocity models. The interpretational difficulty arises from the fact that migration with erroneous velocity does not yield the geologically correct reflector geometries and that individual migrated images suffer from poor signal‐to‐noise ratio. Moreover, migrated events may be of considerable complexity and thus hard to identify. In this paper, we examine the influence of wrong velocity models on the output of prestack depth migration in the case of straight reflector and point diffractor data in homogeneous media. To avoid obscuring migration results by artifacts (“smiles”), we use a geometrical technique for modeling and migration yielding a point‐to‐point map from time‐domain data to depth‐domain data. We discover that strong deformation of migrated events may occur even in situations of simple structures and small velocity errors. From a kinematical point of view, we compare the results of common‐shot and common‐offset migration. and we find that common‐offset migration with erroneous velocity models yields less severe image distortion than common‐shot migration. However, for any kind of migration, it is important to use the entire cube of migrated data to consistently interpret in the prestack depth‐migrated domain.

Migration velocity analysis based on focusing diffractions generated using the CRS wavefront attributes: Application to real land data

Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima
Keyword(s):  
Time Domain ◽  
Real Data ◽  
Normal Incidence ◽  
Velocity Model ◽  
Migration Velocity ◽  
Velocity Analysis ◽  
Velocity Models ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
The Time Domain

In land seismic data processing, the prestack time migration (PSTM) image remains the standard imaging output, but a reliable migrated image of the subsurface depends on the accuracy of the migration velocity model. We have adopted two new algorithms for time-domain migration velocity analysis based on wavefield attributes of the common-reflection-surface (CRS) stack method. These attributes, extracted from multicoverage data, were successfully applied to build the velocity model in the depth domain through tomographic inversion of the normal-incidence-point (NIP) wave. However, there is no practical and reliable method for determining an accurate and geologically consistent time-migration velocity model from these CRS attributes. We introduce an interactive method to determine the migration velocity model in the time domain based on the application of NIP wave attributes and the CRS stacking operator for diffractions, to generate synthetic diffractions on the reflection events of the zero-offset (ZO) CRS stacked section. In the ZO data with diffractions, the poststack time migration (post-STM) is applied with a set of constant velocities, and the migration velocities are then selected through a focusing analysis of the simulated diffractions. We also introduce an algorithm to automatically calculate the migration velocity model from the CRS attributes picked for the main reflection events in the ZO data. We determine the precision of our diffraction focusing velocity analysis and the automatic velocity calculation algorithms using two synthetic models. We also applied them to real 2D land data with low quality and low fold to estimate the time-domain migration velocity model. The velocity models obtained through our methods were validated by applying them in the Kirchhoff PSTM of real data, in which the velocity model from the diffraction focusing analysis provided significant improvements in the quality of the migrated image compared to the legacy image and to the migrated image obtained using the automatically calculated velocity model.

Redatuming through a salt canopy and target-oriented salt-flank imaging

Geophysics ◽  
2008 ◽  
Vol 73 (3) ◽  
pp. S63-S71 ◽  
Author(s):  
Rongrong Lu ◽  
Mark Willis ◽  
Xander Campman ◽  
Jonathan Ajo-Franklin ◽  
M. Nafi Toksöz
Keyword(s):  
Velocity Gradient ◽  
Model Building ◽  
Velocity Structure ◽  
Velocity Model ◽  
Seismic Profile ◽  
Velocity Estimation ◽  
Reverse Time ◽  
Estimation Model ◽  
Depth Migration ◽  
Prestack Migration

We describe a new shortcut strategy for imaging the sediments and salt edge around a salt flank through an overburden salt canopy. We tested its performance and capabilities on 2D synthetic acoustic seismic data from a Gulf of Mexico style model. We first redatumed surface shots, using seismic interferometry, from a walkaway vertical seismic profile survey as if the source and receiver pairs had been located in the borehole at the positions of the receivers. This process creates effective downhole shot gathers by completely moving surface shots through the salt canopy, without any knowledge of overburden velocity structure. After redatuming, we can apply multiple passes of prestack migration from the reference datum of the bore-hole. In our example, first-pass migration, using only a simple vertical velocity gradient model, reveals the outline of the salt edge. A second pass of reverse-time, prestack depth migration using full two-way wave equation was performed with an updated velocity model that consisted of the velocity gradient and salt dome. The second-pass migration brings out dipping sediments abutting the salt flank because these reflectors were illuminated by energy that bounced off the salt flank, forming prismatic reflections. In this target-oriented strategy, the computationally fast redatuming process eliminates the need for the traditional complex process of velocity estimation, model building, and iterative depth migration to remove effects of the salt canopy and surrounding overburden. This might allow this strategy to be used in the field in near real time.

Migration moveout analysis and depth focusing

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 91-100 ◽  
Author(s):  
Claude F. Lafond ◽  
Alan R. Levander
Keyword(s):  
Heterogeneous Media ◽  
A Priori ◽  
Complex Structure ◽  
Tomographic Reconstruction ◽  
Synthetic Data ◽  
Real Data ◽  
Velocity Model ◽  
Velocity Analysis ◽  
Data Set ◽  
Velocity Models

Prestack depth migration still suffers from the problems associated with building appropriate velocity models. The two main after‐migration, before‐stack velocity analysis techniques currently used, depth focusing and residual moveout correction, have found good use in many applications but have also shown their limitations in the case of very complex structures. To address this issue, we have extended the residual moveout analysis technique to the general case of heterogeneous velocity fields and steep dips, while keeping the algorithm robust enough to be of practical use on real data. Our method is not based on analytic expressions for the moveouts and requires no a priori knowledge of the model, but instead uses geometrical ray tracing in heterogeneous media, layer‐stripping migration, and local wavefront analysis to compute residual velocity corrections. These corrections are back projected into the velocity model along raypaths in a way that is similar to tomographic reconstruction. While this approach is more general than existing migration velocity analysis implementations, it is also much more computer intensive and is best used locally around a particularly complex structure. We demonstrate the technique using synthetic data from a model with strong velocity gradients and then apply it to a marine data set to improve the positioning of a major fault.

Sign in / Sign up

Export Citation Format

Share Document