Sensitivity of prestack depth migration to the velocity model

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.


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


Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 594-603 ◽  
Author(s):  
Robert J. Ferguson ◽  
Gary F. Margrave

A new depth migration method suitable for heterogeneous media is presented. The well‐known phase shift plus interpolation (PSPI) method and the recently introduced nonstationary phase‐shift (NSPS) method are combined into a single symmetric operator with improved accuracy and stability and with similar computational effort. For prestack depth migration, the symmetric operator is used in a recursive wavefield extrapolation to compute incident and reflected wavefields at any desired depth, and the ratio of the incident and reflected wavefields at a particular depth is used to estimate seismic reflectivity. When the velocity model is made piecewise constant laterally, the symmetric extrapolation operator can be computed efficiently using ordinary phase‐shift extrapolation for a series of reference velocities and appropriate spatial windowing. Migration of the Marmousi synthetic data set by symmetric nonstationary phase shift (SNPS) provides an image that compares favorably with an image of the zero‐offset reflectivity derived from the Marmousi velocity model.


2019 ◽  
Vol 38 (4) ◽  
pp. 268-273
Author(s):  
Maheswara Phani ◽  
Sushobhan Dutta ◽  
Kondal Reddy ◽  
Sreedurga Somasundaram

Raageshwari Deep Gas (RDG) Field is situated in the southern part of the Barmer Basin in Rajasthan, India, at a depth of 3000 m. With both clastic and volcanic lithologies, the main reservoirs are tight, and hydraulic fracturing is required to enhance productivity, especially to improve permeability through interaction of induced fractures with natural fractures. Therefore, optimal development of the RDG Field reservoirs requires characterization of faults and natural fractures. To address this challenge, a wide-azimuth 3D seismic data set over the RDG Field was processed to sharply define faults and capture anisotropy related to open natural fractures. Anisotropy was indicated by the characteristic sinusoidal nature of gather reflection events processed using conventional tilted transverse imaging (TTI); accordingly, we used orthorhombic imaging to correct for these, to quantify fracture-related anisotropy, and to yield a more correct subsurface image. During prestack depth migration (PSDM) processing of the RDG data, TTI and orthorhombic velocity modeling was done with azimuthal sectoring of reflection arrivals. The resultant PSDM data using this velocity model show substantial improvement in image quality and vertical resolution at the reservoir level compared to vintage seismic data. The improved data quality enabled analysis of specialized seismic attributes like curvature and thinned fault likelihood for more reliable characterization of faults and fractures. These attributes delineate the location and distribution of probable fracture networks within the volcanic reservoirs. Interpreted subtle faults, associated with fracture zones, were validated with microseismic, production, and image log data.


2017 ◽  
Vol 5 (3) ◽  
pp. SJ81-SJ90 ◽  
Author(s):  
Kainan Wang ◽  
Jesse Lomask ◽  
Felix Segovia

Well-log-to-seismic tying is a key step in many interpretation workflows for oil and gas exploration. Synthetic seismic traces from the wells are often manually tied to seismic data; this process can be very time consuming and, in some cases, inaccurate. Automatic methods, such as dynamic time warping (DTW), can match synthetic traces to seismic data. Although these methods are extremely fast, they tend to create interval velocities that are not geologically realistic. We have described the modification of DTW to create a blocked dynamic warping (BDW) method. BDW generates an automatic, optimal well tie that honors geologically consistent velocity constraints. Consequently, it results in updated velocities that are more realistic than other methods. BDW constrains the updated velocity to be constant or linearly variable inside each geologic layer. With an optimal correlation between synthetic seismograms and surface seismic data, this algorithm returns an automatically updated time-depth curve and an updated interval velocity model that still retains the original geologic velocity boundaries. In other words, the algorithm finds the optimal solution for tying the synthetic to the seismic data while restricting the interval velocity changes to coincide with the initial input blocking. We have determined the application of the BDW technique on a synthetic data example and field data set.


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


Sign in / Sign up

Export Citation Format

Share Document