Integrated prestack depth migration of vertical seismic profile and surface seismic data from the Rocky Mountain Foothills of southern Alberta, Canada

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 ◽  
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 ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. SI15-SI26 ◽  
Author(s):  
Xiang Xiao ◽  
Gerard T. Schuster

We have developed a novel vertical seismic profile (VSP) imaging method that requires only a local velocity model around the well. This method is denoted as a “local migration,” where the model-based forward-propagation and backward-propagation operators are computed using a local velocity model between the well and the target body. The velocity model for the complex overburden and salt body is not needed, and the source-side statics are automatically accounted for. In addition, kinematic and dynamic effects, including anisotropy, absorption, and other unknown/undefined rock effects outside of this local velocity model, are mostly accounted for. Numerical tests on an acoustic 2D Sigsbee VSP data set, an elastic salt model, and offset 2D VSP data in the Gulf of Mexico partly validate the effectiveness of this method by accurately imaging the sediments next to the vertical well.


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.


1984 ◽  
Vol 24 (1) ◽  
pp. 429
Author(s):  
F. Sandnes W. L. Nutt ◽  
S. G. Henry

The improvement of acquisition and processing techniques has made it possible to study seismic wavetrains in boreholes.With careful acquisition procedures and quantitative data processing, one can extract useful information on the propagation of seismic events through the earth, on generation of multiples and on the different reflections coming from horizons that may not all be accessible by surface seismic.An extensive borehole seismic survey was conducted in a well in Conoco's contract area 'Block B' in the South China Sea. Shots at 96 levels were recorded, and the resulting Vertical Seismic Profile (VSP) was carefully processed and analyzed together with the Synthetic Seismogram (Geogram*) and the Synthetic Vertical Seismic Profile (Synthetic VSP).In addition to the general interpretation of the VSP data, i.e. time calibration of surface seismic, fault identification, VSP trace inversion and VSP Direct Signal Analysis, the practical inclusion of VSP data in the reprocessing of surface seismic data was studied. Conclusions that can be drawn are that deconvolution of surface seismic data using VSP data must be carefully approached and that VSP can be successfully used to examine phase relationships in seismic data.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. WB219-WB224 ◽  
Author(s):  
Weiping Cao ◽  
Gerard T. Schuster

An antialiasing formula has been derived for interferometric redatuming of seismic data. More generally, this formula is valid for numerical implementation of the reciprocity equation of correlation type, which is used for redatuming, extrapolation, interpolation, and migration. The antialiasing condition can be, surprisingly, more tolerant of a coarser trace sampling compared to the standard antialiasing condition. Numerical results with synthetic vertical seismic profile (VSP) data show that interferometry artifacts are effectively reduced when the antialiasing condition is used as a constraint with interferometric redatuming.


GeoArabia ◽  
1999 ◽  
Vol 4 (3) ◽  
pp. 363-378
Author(s):  
Mohammed A. Badri ◽  
Taha M. Taha ◽  
Robert W. Wiley

ABSTRACT In 1995 oil was discovered in the pre-Miocene Matulla and Nubia Sandstones in the Ras El Ush field, Gulf of Suez, Egypt. The discovery was based on an aeromagnetic anomaly from a basement high. After drilling several delineation wells, based on a geological model, it became evident that the field is very complex as it is broken into tilted and rotated compartmental blocks by two perpendicular fault systems. Also the 2-D seismic data were of poor quality beneath the thick Miocene South Gharib Evaporite. Since part of the field lies below shallow-water, 3-D seismic was considered to be too costly. When a delineation well did not encounter the reservoir, due to an unanticipated fault, a 2-D walkaway Vertical Seismic Profile (VSP) was acquired. It clearly revealed the presence of a cross fault. The success of the 2-D VSP in imaging the fault led to the acquisition of the first Middle East 3-D VSP survey in the following well. A downhole, tri-axial, five geophone array tool was used to acquire the 3-D VSP. The 3-D volume of the final migrated VSP data provided the means for the reliable mapping of horizons beneath the South Gharib Evaporite. These maps improved the definition of the field and helped detect previously unrecognized prospective blocks. Four further successful delineation wells confirmed the 3-D VSP interpretation.


Geophysics ◽  
1993 ◽  
Vol 58 (11) ◽  
pp. 1676-1688
Author(s):  
Ronald C. Hinds ◽  
Neil L. Anderson ◽  
Richard Kuzmiski

On the basis of conventional surface seismic data, the 13–15–63–25W5M exploratory well was drilled into a low‐relief Leduc Formation reef (Devonian Woodbend Group) in the Simonette area, west‐central Alberta, Canada. The well was expected to intersect the crest of the reef and encounter about 50–60 m of pay; unfortunately it was drilled into a flank position and abandoned. The decision to abandon the well, as opposed to whipstocking in the direction of the reef crest, was made after the acquisition and interpretive processing of both near( and far‐offset (252 and 524 m, respectively) vertical seismic profile (VSP) data, and after the reanalysis of existing surface seismic data. The near‐ and far‐offset VSPs were run and interpreted while the drill rig remained on‐site, with the immediate objectives of: (1) determining an accurate tie between the surface seismic data and the subsurface geology; and (2) mapping relief along the top of the reef over a distance of 150 m from the 13–15 well location in the direction of the adjacent productive 16–16 well (with a view to whipstocking). These surveys proved to be cost‐effective in that the operators were able to determine that the crest of the reef was out of the target area, and that whipstocking was not a viable alternative. The use of VSP surveys in this situation allowed the operators to avoid the costs associated with whipstocking, and to feel confident with their decision to abandon the well.


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 ◽  
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 ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE13-VE23 ◽  
Author(s):  
Frank Adler ◽  
Reda Baina ◽  
Mohamed Amine Soudani ◽  
Pierre Cardon ◽  
Jean-Baptiste Richard

Velocity-model estimation with seismic reflection tomography is a nonlinear inverse problem. We present a new method for solving the nonlinear tomographic inverse problem using 3D prestack-depth-migrated reflections as the input data, i.e., our method requires that prestack depth migration (PSDM) be performed before tomography. The method is applicable to any type of seismic data acquisition that permits seismic imaging with Kirchhoff PSDM. A fundamental concept of the method is that we dissociate the possibly incorrect initial migration velocity model from the tomographic velocity model. We take the initial migration velocity model and the residual moveout in the associated PSDM common-image gathers as the reference data. This allows us to consider the migrated depth of the initial PSDM as the invariant observation for the tomographic inverse problem. We can therefore formulate the inverse problem within the general framework of inverse theory as a nonlinear least-squares data fitting between observed and modeled migrated depth. The modeled migrated depth is calculated by ray tracing in the tomographic model, followed by a finite-offset map migration in the initial migration model. The inverse problem is solved iteratively with a Gauss-Newton algorithm. We applied the method to a North Sea data set to build an anisotropic layer velocity model.


Sign in / Sign up

Export Citation Format

Share Document