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

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 ◽  
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 ◽  
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.


Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 1782-1791 ◽  
Author(s):  
M. Graziella Kirtland Grech ◽  
Don C. Lawton ◽  
Scott Cheadle

We have developed an anisotropic prestack depth migration code that can migrate either vertical seismic profile (VSP) or surface seismic data. We use this migration code in a new method for integrated VSP and surface seismic depth imaging. Instead of splicing the VSP image into the section derived from surface seismic data, we use the same migration algorithm and a single velocity model to migrate both data sets to a common output grid. We then scale and sum the two images to yield one integrated depth‐migrated section. After testing this method on synthetic surface seismic and VSP data, we applied it to field data from a 2D surface seismic line and a multioffset VSP from the Rocky Mountain Foothills of southern Alberta, Canada. Our results show that the resulting integrated image exhibits significant improvement over that obtained from (a) the migration of either data set alone or (b) the conventional splicing approach. The integrated image uses the broader frequency bandwidth of the VSP data to provide higher vertical resolution than the migration of the surface seismic data. The integrated image also shows enhanced structural detail, since no part of the surface seismic section is eliminated, and good event continuity through the use of a single migration–velocity model, obtained by an integrated interpretation of borehole and surface seismic data. This enhanced migrated image enabled us to perform a more robust interpretation with good well ties.


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 ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 533-546 ◽  
Author(s):  
Robert G. Clapp ◽  
Biondo L. Biondi ◽  
Jon F. Claerbout

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.


2014 ◽  
Vol 69 (6) ◽  
Author(s):  
Sudra Irawan ◽  
Sismanto Sismanto ◽  
Adang Sukmatiawan

Seismic data processing is one of the three stages in the seismic method that has an important role in the exploration of oil and gas. Without good data processing, it is impossible to get seismic image cross section for good interpretation. A research using seismic data processing was done to update the velocity model by horizon based tomography method in SBI Field, North West Java Basin. This method reduces error of seismic wave travel time through the analyzed horizon because the existence velocity of high lateral variation in research area. There are three parameters used to determine the accuracy of the resulting interval velocity model, namely, flat depth gathers, semblance residual moveout that coincides with the axis zero residual moveout, and the correspondence between image depth (horizon) with wells marker  (well seismic tie). Pre Stack Depth Migration (PSDM) form interval velocity model and updating using horizon-based tomography method gives better imaging of under-surfaced structure results than PSDM before using tomography. There are three faults found in the research area, two normal faults have southwest-northeast strike and the other has northwest-southeast strike. The thickness of reservoir in SBI field, North West Java Basin, is predicted between 71 to 175 meters and the hydrocarbon (oil) reserve is predicted about  with 22.6% porosity and 70.7% water saturation. 


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.


2010 ◽  
Vol 50 (2) ◽  
pp. 723
Author(s):  
Sergey Birdus ◽  
Erika Angerer ◽  
Iftikhar Abassi

Processing of multi and wide-azimuth seismic data faces some new challenges, and one of them is depth-velocity modelling and imaging with azimuthal velocity anisotropy. Analysis of multi-azimuth data very often reveals noticeable fluctuations in moveout between different acquisition directions. They can be caused by several factors: real azimuthal interval velocity anisotropy associated with quasi-vertical fractures or present day stress field within the sediments; short-wavelength velocity heterogeneities in the overburden; TTI (or VTI) anisotropy in the overburden; or, random distortions due to noise, multiples, irregularities in the acquisition geometry, etcetera. In order to build a velocity model for multi-azimuth pre-stack depth migration (MAZ PSDM) taking into account observed azimuthal anisotropy, we need to recognise, separate and estimate all the effects listed above during iterative depth-velocity modelling. Analysis of seismic data from a full azimuth 3D seismic land survey revealed the presence of strong spatially variable azimuthal velocity anisotropy that had to be taken into consideration. Using real data examples we discuss major steps in depth processing workflow that took such anisotropy into account: residual moveout estimation in azimuth sectors; separation of different effects causing apparent azimuthal anisotropy (see A–D above); iterative depth-velocity modelling with azimuthal anisotropy; and, subsequent MAZ anisotropic PSDM. The presented workflow solved problems with azimuthal anisotropy in our multi-azimuth dataset. Some of the lessons learned during this MAZ project are relevant to every standard narrow azimuth seismic survey recorded in complex geological settings.


2020 ◽  
Vol 82 (6) ◽  
pp. 29-37
Author(s):  
Sudra Irawan ◽  
Siti Noor Chayati ◽  
Sismanto Sismanto

The tomography method requires an excellent initial velocity model. On the horizon based tomography, it will correct the travel time error of seismic waves along the horizon which is analysed using input results from the analysis of residual depth moveout. In this study, a semblance residual moveout analysis will be conducted after the interval velocity model has applied to the SBI field seismic data (CDP Gathers and RMS velocity). Based on the imaging results generated by the PSDM running process, an aperture value of 550 for inline and 800 for crossline is selected. PSDM generated from the initial interval velocity model has an acoustic impedance value between 1000 kg/m2s to 14339.2 kg/m2s. The PSDM process, residual moveout analysis, and horizon-based tomography are carried out iteratively until the error in the interval velocity model approaches zero. In this study, five iterations were performed. The resulting residual moveout is increasingly oscillating around zero after the 5th iteration, which indicates that the error in the interval velocity model is getting smaller. There are two types of residual moveout, namely residuals moveout positively and residuals moveout negatively. Residual moveout positive indicates that the velocity used is too high, while the residual moveout negative indicates that the velocity used is too low. The identification of interval velocity model errors with analysis of residual moveout semblance is calculated from depth gathers. The semblance residual moveout analysis is used for the Pre Stack Depth Migration (PSDM) depth image analysis stage along with the marker (well data). .


Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1903-1912 ◽  
Author(s):  
Igor B. Morozov ◽  
Alan Levander

Wide‐aperture, prestack depth migration requires application of challenging and time‐consuming velocity analysis and depth focusing, collectively referred to here as depth focusing. We present an approach to depth focusing using (1) a detailed starting velocity model obtained by a 1‐D transformation of the first‐arrival times, followed iteratively by (2) interactive analysis of the common‐image gathers, (3) computation of coherency attributes of the wavefield in the depth domain, and (4) 2‐D traveltime tomography to update the background velocity model. We employ two interactive method of migration velocities refinement. In the first method (similar to the common‐midpoint velocity spectrum approach), residual velocity updates are picked directly from the common‐image gathers. In another method (analogous to the common velocity stacks), we pick the velocity updates from the areas of maximum coherency in depth sections that are migrated using rescaled traveltime maps. Both types of migration velocity picks, optionally combined with the first arrivals, are inputs for a 2‐D traveltime inversion scheme that uses either the infinite‐frequency or a finite‐bandwidth approximation. This flexible and versatile depth focusing approach is implemented for several prestack depth migration algorithms and illustrated on an application to a real, ultrashallow seismic data set. The technique resolves overburden velocity variations and facilitates reliable high‐resolution reflection imaging of a paleochannel that was the target of the study.


Sign in / Sign up

Export Citation Format

Share Document