Prestack depth migration of primary and surface-related multiple reflections: Part II — Identification and removal of residual multiples

Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. S71-S76 ◽  
Author(s):  
Remco Muijs ◽  
Johan O. Robertsson ◽  
Klaus Holliger

Depth imaging using primary and multiple reflections (DIPMR), as described in Part I of this study, allows subsurface information carried by multiple reflections to be utilized. In the presence of strong lateral heterogeneity, however, the migration results may be distorted by artifacts originating from reflections associated with layers above the imaging plane that are commonly referred to as crosstalk. We present an image enhancement procedure that allows such artifacts to be effectively suppressed by predicting the initial crosstalk in a second migration phase. This second migration uses reflections imaged at shallower depth levels and requires knowledge of the total and primary downgoing wavefield at the receiver level. The predicted crosstalk image it-self may assist the interpretational effort by indicating areas where artifacts may result in incorrect identification of geologic structures or may cause local distortions of amplitude informa-tion. Furthermore, a clean depth image can be obtained by adap-tively subtracting the predicted crosstalk-related noise from the original image. The proposed method is an extension of the DIPMR procedure outlined in Part I of this study. However, it is also applicable to seismic data imaged by using conventional mi gration techniques that are based exclusively on primary reflec-tions. In the latter case, the enhancement procedure allows the ef-fects of imperfect multiple removal during preprocessing to be revealed in the migrated sections. When applied to synthetic data simulated for a complex salt model, the proposed enhancement procedure proved to be valid and effective.

Geophysics ◽  
1993 ◽  
Vol 58 (6) ◽  
pp. 873-882 ◽  
Author(s):  
Roelof Jan Versteeg

To get a correct earth image from seismic data acquired over complex structures it is essential to use prestack depth migration. A necessary condition for obtaining a correct image is that the prestack depth migration is done with an accurate velocity model. In cases where we need to use prestack depth migration determination of such a model using conventional methods does not give satisfactory results. Thus, new iterative methods for velocity model determination have been developed. The convergence of these methods can be accelerated by defining constraints on the model in such a way that the method only looks for those components of the true earth velocity field that influence the migrated image. In order to determine these components, the sensitivity of the prestack depth migration result to the velocity model is examined using a complex synthetic data set (the Marmousi data set) for which the exact model is known. The images obtained with increasingly smoothed versions of the true model are compared, and it is shown that the minimal spatial wavelength that needs to be in the model to obtain an accurate depth image from the data set is of the order of 200 m. The model space that has to be examined to find an accurate velocity model from complex seismic data can thus be constrained. This will increase the speed and probability of convergence of iterative velocity model determination methods.


2005 ◽  
Vol 45 (1) ◽  
pp. 421
Author(s):  
P. Bocca ◽  
L. Fava ◽  
E. Stolf

3D pre-stack depth migration (PSDM) reprocessing was conducted in 2003 on a portion of the Onnia 3D seismic cube, located in exploration permit AC/P-21, Timor Sea.The main objective of the reprocessing was to obtain the best seismic depth image and the most realistic structural reconstruction of the sub-surface to mitigate the risk factors associated with trap definition (trap retention and trap efficiency). This represents one of the main challenges for oil exploration in the area.The 3D PSDM methodology was chosen as the most appropriate imaging tool to define the correct sub-surface geometry and fault imaging through the use of an appropriate velocity field. An integrated approach to building the final velocity model was adopted, with a substantial contribution from the regional geological model.Several examples are given to demonstrate that the 3D PSDM reprocessing significantly improved the seismic image and thus the confidence in the interpretation, contributing strongly to the definition of the exploration targets.The interpretation of the new seismic data has resulted in a new structural picture in which higher confidence in seismic imaging has improved fault correlation. This has enabled better structural definition at the Middle Jurassic Plover Formation level that has reduced the complexity of the large Vesta Prospect, in the centre of the Swan Graben to the northwest of East Swan–1. Improved understanding of the fault reactivation mechanism and the structural elements of the trap (trap integrity) were eventually incorporated in the prospect risking.In the Swan Graben 3D PSDM has proved to be a very powerful instrument capable of producing significant impact on the exploration even in an area with a complex geological setting and a fairly poor seismic data quality.


Geophysics ◽  
1967 ◽  
Vol 32 (2) ◽  
pp. 207-224 ◽  
Author(s):  
John D. Marr ◽  
Edward F. Zagst

The more recent developments in common‐depth‐point techniques to attenuate multiple reflections have resulted in an exploration capability comparable to the development of the seismic reflection method. The combination of new concepts in digital seismic data processing with CDP techniques is creating unforeseen exploration horizons with vastly improved seismic data. Major improvements in multiple reflection and reverberation attenuation are now attainable with appropriate CDP geometry and special CDP stacking procedures. Further major improvements are clearly evident in the very near future with the use of multichannel digital filtering‐stacking techniques and the application of deconvolution as the first step in seismic data processing. CDP techniques are briefly reviewed and evaluated with real and experimental data. Synthetic data are used to illustrate that all seismic reflection data should be deconvolved as the first processing step.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE269-VE280 ◽  
Author(s):  
Priyank Jaiswal ◽  
Colin A. Zelt

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.


Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 568-576 ◽  
Author(s):  
Young C. Kim ◽  
Worth B. Hurt, ◽  
Louis J. Maher ◽  
Patrick J. Starich

The transformation of surface seismic data into a subsurface image can be separated into two components—focusing and positioning. Focusing is associated with ensuring the data from different offsets are contributing constructively to the same event. Positioning involves the transformation of the focused events into a depth image consistent with a given velocity model. In prestack depth migration, both of these operations are achieved simultaneously; however, for 3-D data, the cost is significant. Prestack time migration is much more economical and focuses events well even in the presence of moderate velocity variations, but suffers from mispositioning problems. Hybrid migration is a cost‐effective depth‐imaging approach that uses prestack time migration for focusing; inverse migration for the removal of positioning errors; and poststack depth migration for proper positioning. When lateral velocity changes are moderate, the hybrid technique can generate a depth image that is consistent with a velocity field. For very complex structures that require prestack depth migration, the results of the hybrid technique can be used to create a starting velocity model, thereby reducing the number of iterations for velocity model building.


2017 ◽  
Vol 5 (1) ◽  
pp. T1-T9 ◽  
Author(s):  
Rui Zhang ◽  
Kui Zhang ◽  
Jude E. Alekhue

More and more seismic surveys produce 3D seismic images in the depth domain by using prestack depth migration methods, which can present a direct subsurface structure in the depth domain rather than in the time domain. This leads to the increasing need for applications of seismic inversion on the depth-imaged seismic data for reservoir characterization. To address this issue, we have developed a depth-domain seismic inversion method by using the compressed sensing technique with output of reflectivity and band-limited impedance without conversion to the time domain. The formulations of the seismic inversion in the depth domain are similar to time-domain methods, but they implement all the elements in depth domain, for example, a depth-domain seismic well tie. The developed method was first tested on synthetic data, showing great improvement of the resolution on inverted reflectivity. We later applied the method on a depth-migrated field data with well-log data validated, showing a great fit between them and also improved resolution on the inversion results, which demonstrates the feasibility and reliability of the proposed method on depth-domain seismic data.


Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1239-1246 ◽  
Author(s):  
Robert W. Vestrum ◽  
Don C. Lawton ◽  
Ron Schmid

Seismic anisotropy in dipping shales causes imaging and positioning problems for underlying structures. We developed an anisotropic depth‐migration approach for P-wave seismic data in transversely isotropic (TI) media with a tilted axis of symmetry normal to bedding. We added anisotropic and dip parameters to the depth‐imaging velocity model and used prestack depth‐migrated image gathers in a diagnostic manner to refine the anisotropic velocity model. The apparent position of structures below dipping anisotropic overburden changes considerably between isotropic and anisotropic migrations. The ray‐tracing algorithm used in a 2-D prestack Kirchhoff depth migration was modified to calculate traveltimes in the presence of TI media with a tilted symmetry axis. The resulting anisotropic depth‐migration algorithm was applied to physical‐model seismic data and field seismic data from the Canadian Rocky Mountain Thrust and Fold Belt. The anisotropic depth migrations offer significant improvements in positioning and reflector continuity over those obtained using isotropic algorithms.


Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. B87-B96 ◽  
Author(s):  
Ammanuel Fesseha Woldearegay ◽  
Priyank Jaiswal ◽  
Alexander R. Simms ◽  
Hanna Alexander ◽  
Leland C. Bement ◽  
...  

Depth imaging in ultrashallow ([Formula: see text]) environments presents twofold challenge: (1) coda available for depth migration is very limited; and (2) conventional time processing with limited coda generally fails to estimate reliable velocity models for depth migration. We studied the combining of first-arrival traveltime inversion and prestack depth migration (PSDM) for depth imaging of ultrashallow paleochannel stratigraphy associated with the Bull Creek drainage system, Oklahoma. Restricted by a limited number of geophones (24) we acquired data for inversion and migration through two coincident profiles. The first profile for inversion has a wider survey-aperture (115-m maximum shot-receiver spacing) and consequently sparse CMP spacing (2.5 m), whereas the second profile for PSDM has denser CMP spacing (1 m) and consequently a narrower survey aperture (46-m maximum shot-receiver spacing). We also found that the velocity model from traveltime inversion of the wider-aperture data set is more preferable for depth-migration than the velocity model from time processing of the denser data set. The preferred depth image showed three episodes of incision whose chronological order is resolved through radio-carbon dating of terrace sediments. Results suggested that even with limited geophones, depth imaging of ultrashallow targets can be achieved by combining first-arrival traveltime inversion and PSDM through coincident wide- and narrow-aperture acquisitions.


Geophysics ◽  
2008 ◽  
Vol 73 (1) ◽  
pp. V1-V9 ◽  
Author(s):  
Chun-Feng Li ◽  
Christopher Liner

Although the passage of singularity information from acoustic impedance to seismic traces is now well understood, it remains unanswered how routine seismic processing, mode conversions, and multiple reflections can affect the singularity analysis of surface seismic data. We make theoretical investigations on the transition of singularity behaviors from acoustic impedances to surface seismic data. We also perform numerical, wavelet-based singularity analysis on an elastic synthetic data set that is processed through routine seismic processing steps (such as stacking and migration) and that contains mode conversions, multiple reflections, and other wave-equation effects. Theoretically, seismic traces can be approximated as proportional to a smoothed version of the [Formula: see text] derivative of acoustic impedance,where [Formula: see text] is the vanishing moment of the seismic wavelet. This theoretical approach forms the basis of linking singularity exponents (Hölder exponents) in acoustic impedance with those computable from seismic data. By using wavelet-based multiscale analysis with complex Morlet wavelets, we can estimate singularity strengths and localities in subsurface impedance directly from surface seismic data. Our results indicate that rich singularity information in acoustic impedance variations can be preserved by surface seismic data despite data-acquisition and processing activities. We also show that high-resolution detection of singularities from real surface seismic data can be achieved with a proper choice of the scale of the mother wavelet in the wavelet transform. Singularity detection from surface seismic data thus can play a key role in stratigraphic analysis and acoustic impedance inversion.


Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S73-S79
Author(s):  
Ørjan Pedersen ◽  
Sverre Brandsberg-Dahl ◽  
Bjørn Ursin

One-way wavefield extrapolation methods are used routinely in 3D depth migration algorithms for seismic data. Due to their efficient computer implementations, such one-way methods have become increasingly popular and a wide variety of methods have been introduced. In salt provinces, the migration algorithms must be able to handle large velocity contrasts because the velocities in salt are generally much higher than in the surrounding sediments. This can be a challenge for one-way wavefield extrapolation methods. We present a depth migration method using one-way propagators within lateral windows for handling the large velocity contrasts associated with salt-sediment interfaces. Using adaptive windowing, we can handle large perturbations locally in a similar manner as the beamlet propagator, thus limiting the impact of the errors on the global wavefield. We demonstrate the performance of our method by applying it to synthetic data from the 2D SEG/EAGE [Formula: see text] salt model and an offshore real data example.


Sign in / Sign up

Export Citation Format

Share Document